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Executive  Summary 


This  report  presents  the  results  of  a  study  performed  by  the  Center  for  Transportation 
Research  (CTR)  at  Virginia  Polytechnic  Institute  and  State  University  concerning  the 
development  of  rapid  runway  turnoffs  to  be  used  in  existing  and  future  airport  scenarios  This 
study  was  conducted  for  the  Federal  Aviation  Administration  (FAA)  Office  of  Airport  Research 
and  Development  (ARD-200)  to  assess  the  impact  of  optimal  runway  turnoff  locations  and  their 
corresponding  geometries  in  runway  occupancy  time  and  ultimately  in  runway  capacity.  The 
report  emphasizes  in  the  development  of  a  combined  simulation  and  optimization 
methodology  to  ascertain  the  impact  of  runway  turnoff  placement  in  the  weighted  average 
runway  occupancy  time  (WAROT)  for  realistic  aircraft  populations  operating  from  a  single 
runway.  The  methodology  developed  was  extended  into  a  user-friendly  computer  program 
called  REDIM  -  Runway  Exit  Design  Interactive  Model  -  to  estimate  the  WAROT  performance 
index  for  user  defined  aircraft  populations  under  various  airport  atmospheric  conditions.  The 
accompanying  computer  software  developed  in  this  research  requires  a  minimum  of  512 
Kbytes  of  memory  and  an  ordinary  IBM  or  compatible  computer. 

The  results  of  this  Investigation  can  be  summarized  as  follows: 

1)  Evidence  suggests  that  existing  runway  turnoff  geometries  are  not  being  used  near  their 
design  speeds.  The  small  existing  database  in  high-speed  runway  exit  use  indicates  that 
commercial  aircraft  are  regularly  using  these  exits  10-15  knots  below  their  design  speeds. 
The  major  consequence  of  this  being  higher  runway  occupancy  times  as  aircraft  have  to 
spend  more  time  decelerating  on  the  runway. 
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2)  Significant  reductions  in  runway  occupancy  time  (ROT)  were  observed  (i.e.,  from  7  to  18 
%)  with  the  use  of  three,  four  and  five  optimally  located  runway  turnoffs.  The  magnitude  of 
these  ROT  reductions  was,  as  expected,  a  strong  function  of  the  exit  speeds  used. 

3)  For  all  practical  purposes  a  gobd  approximation  of  a  fully  variable  variable  turnoff  geometry 
can  be  obtained  with  two  radii  of  curvature  ending  at  a  generally  small  exit  angle.  The 
differences  observed  between  the  approximation  and  the  actual  fully  variable  geometries  are 
insignificant  for  small  exit  angles  and  should  not  compromise  safety  and  comfort  factors  while 
maneuvering  on  the  turnoff. 

4)  It  seems  feasible  to  design  rapid  runway  turnoff  designs  for  entry  speeds  of  up  to  35  m/s. 
(78  MPH)  for  aircraft  TERP  categories  C  and  D.  For  these  designs  a  smaller  exit  angle  would 
be  required  to  increase  the  deceleration  distance  available  along  the  tangent  turnoff 
geometry.  This  will  add  confidence  to  the  pilot  before  arriving  to  the  nearest  taxiway  junction. 
An  alternative  to  this  reduction  in  exit  angle  is  to  increase  the  current  lateral  separation 
minima  between  a  runway  and  parallel  taxiway  centerline  distance  from  183  m.  (600  ft.)  to  228 
m.  (750  ft ) 

5)  It  is  suggested  that  further  investigations  be  made  on  the  aircraft  stability  and  control 
implications  at  speeds  above  35  m/s.  to  determine  safety  boundaries  while  executing  turns 
at  moderate  ground  turn  rates. 

B)  Very  low  exit  angles  (i  e.,  <  15  degrees)  should  be  avoided  as  they  contribute  to  large 
turnoff  and  runway  occupancy  times  due  to  long  exposure  of  the  aircraft  wing  or  tailplane  tip 
within  the  runway  bounds.  Optimal  exit  angle  geometries  seem  to  be  in  the  17-30  degree 
range  depending  on  the  specific  aircraft  and  exit  speed.  The  lower  exit  angles  might  have 
better  acceptance  among  pilots  as  they  would  provide  added  distances  to  decelerate  before 
a  taxiway  is  reached.  However,  they  would  involve  slightly  higher  costs  and  land  use. 

7)  A  computer  simulation/optimization  model  to  assess  the  runway  optimal  location  and 
geometries  of  runway  turnoffs  was  successfully  implemented  in  a  personal  computer.  The 
program  requires  a  definition  of  the  aircraft  population  and  environmental  conditions  of  the 
airport  facility  to  optimize  the  weighted  average  runway  occupancy  time  (WAROT).  The 
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program  executes  a  typical  linear  network  optimal  assignment  problem  including  the 
simulation  of  the  aircraft  kinematic  trajectories.  The  program  is  suitable  for  use  at  existing 
facilities  in  process  of  upgrading  as  well  a«  for  entirely  new  runways  planned  for  the  future. 


8)  Preliminary  results  at  the  runway  capacity  level  indicate  that  rapid  runway  designs  could 
increase  runway  hourly  capacity  by  up  to  10-12%  if  further  reductions  are  accomplished  in  the 
in-trail  separation  distances  allowed  for  successive  arrivals  planned  for  future  airport 
environments  (i.e.,  2/4/5  nautical  mile  rule). 
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1.0  Introduction 


The  subject  of  airport  congestion  and  delays  has  received  a  great  deal  of  attention  in 
recent  years  due  to  the  rapid  growth  of  air  transportation  services  coupled  with  a  relatively 
stagnant  airport  infrastructure.  Current  statistics  indicate  that  approximately  two  billion 
dollars  are  paid  by  air  travelers  due  to  system  imposed  delays  in  the  United  States  alone 
[FAA,  1987].  These  delays  are  likely  to  increase  as  air  travel  demand  builds  up  from  416 
billion  passenger  revenue  mile  flown  in  1988  to  an  estimated  750  billion  passenger  revenue 
mile  by  the  end  of  the  century  or  equivalent  to  an  average  annual  growth  of  roughly  6  percent 
[Aviation  Week  &  Space  Technology,  1989],  The  problem  is  further  aggravated  when  one 
considers  that  the  current  air  transportation  system  has  been  operating  in  a  pseudo-stagnant 
mode  with  almost  the  same  infrastructure  in  terms  of  airport  facilities  since  the  early 
seventies.  The  Federal  Aviation  Administration  (FAA)  has  estimated  that  eleven  major 
airports  now  experience  severe  chronic  operational  delays  -more  than  20,000  hours  of  system 
imposed  delays  per  year  -  as  a  result  of  traffic  congestion  [FAA,  1989],  According  to  the  FAA 
this  number  will  increase  to  thirty  two  by  the  year  1996  and  possibly  fifty  by  the  end  of  the 
century  [FAA,  1988],  One  fifth  of  these  airport  facilities  will  experience  more  than  50,000  hours 
of  system  imposed  delays  according  to  the  same  study.  Delays  at  these  key  airports  are  not 
simply  local  problems;  the  effects  ripple  outward  to  other  airports  with  flights  connecting  to 
these  hubs  and  ultimately  to  the  entire  air  transportation  network  [Transportation  Research 
Board.  1988], 

To  illustrate  the  magnitude  of  the  capacity  restrictions  at  a  typical  major  airport  facility 
consider  the  Los  Angeles  International  Airport  (LAX).  According  to  the  City  of  Los  Angeles 
Department  of  Airports,  Facilities  Planning  Bureau,  the  design  capacities  of  LAX  are  147  and 
128  aircraft  operations  per  hour  under  Visual  and  Instrument  Flight  Rules  (VFR  and  IFR), 
respectively  [Los  Angeles  City  Department  of  Airports,  1988].  Due  to  strict  noise  abatement 
procedures  this  capacity  is  reduced  to  114  operations  per  hour  for  both  VFR  and  IFR 
conditions.  Furthermore,  during  the  hours  of  midnight  to  6:30  A.M.  the  terminal  airport  capacity 
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is  restricted  to  32  operations  per  hour.  This  last  point  has  tremendous  operational  effects  on 
this  facility  due  to  its  strategic  location  for  transpacific  flights.  Similar  capacity  restrictions 
have  been  applied  to  Chicago  O'Hare  International  Airport  prompted  by  a  targe  increase  in 
the  number  of  operational  errors  at  the  terminal  control  center  [Aviation  Week  and  Space 
Technology,  1988). 

The  capacity  of  an  airport  facility  is  dictated  by  the  critical  capacity  of  the  following  four 
components:  1)  the  airspace,  2)  the  runway,  3)  the  taxiway,  and  4)  the  apron-gate  component 
(i.e.,  assuming  a  well  designed  ground  access  system).  Although  the  interrelations  between 
these  four  components  could  be  significant  for  certain  airport  configurations  it  has  been 
customary  to  study  in  detail  each  component  independently  and  then  select  the  most 
restrictive  one  as  that  defining  the  capacity  of  the  facility.  Obviously,  capacity  is  also  affected 
by  many  external  factors  such  as  meteorological  conditions,  airfield  configuration,  aircraft 
characteristics,  and  air  traffic  control  system  performance. 

The  challenge  faced  by  today's  system  developers  is  how  to  increase  system  capacity 
without  violating  present  operational  norms  and  degrading  system  safety.  Research  and 
development  programs  at  FAA  and  NASA  are  addressing  several  issues  of  the  airfield 
problems,  foremost  of  which  are  the  improvement  of  operational  use  of  runways,  provision 
of  efficient  flow  control,  spacing  and  management  of  aircraft  in  the  terminal  airspace, 
upgrading  of  the  computer/communication  technology  usage,  and  the  resolution  of  the  effects 
of  wake  vortex  and  aircraft  noise. 

Improving  the  operational  use  of  runways,  the  reduction  of  runway  occupancy  times  using 
high-speed  exits  is  one  of  the  research  activities  carried  by  the  FAA  and  NASA.  The  efficiency 
of  the  runway  component  is  dictated  primarily  by  the  runway  occupancy  time  (ROT)  and  its 
variability  from  aircraft  to  aircraft.  ROT  is  the  time  an  aircraft  spends  on  the  runway  or  its 
vicinity  until  a  new  arrival  or  departure  can  be  processed.  Table  1.1  illustrates  typical  values 
of  ROT  and  its  variablity  for  several  aircraft  classes  using  current  and  future  technologies 
according  to  a  recent  study  conducted  at  the  MITRE  Corporation  [Barrer  and  Dielh,  1988).  The 
study  quantified  the  potential  increases  to  runway  capacity  resulting  from  improvements  to  the 
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Table  1.1  Arrival  Runway  Occupancy  Times  (Adapted  from  Barter  &  Dlelh,  1988). 

TERP  Category 

ROT  Mean  Val^e  (seconds) 

Present  Future  (1996) 

A 

43 

35 

B 

45 

37 

C 

46 

40 

D 

50 

45 

Air  Traffic  Control  System  (ATC)  performance  parameters  (i.e.,  reducing  in-trail  landing 
separations,  better  planned  runway  exits,  improved  ground-based  radar  surveillance 
capability,  etc.)  and  concluded  that  gains  of  up  to  20%  in  the  capacity  of  a  single  runway  are 
possible  if  these  control  actions  were  to  be  implemented.  Other  studies  support  similar  gains 
if  advanced  systems  are  used  [Lebron,  1987;  Simpson  et  al,  1988).  On  the  airspace  component 
a  critical  parameter  directly  related  to  capacity  is  the  arrival  in-trial  separation  (AIT).  AIT 
values  are  dictated  by  safety  criteria  to  avoid  the  effects  of  the  wake  vortex  generated  by  a 
leading  aircraft  during  landing.  Current  and  future  arrival  in-trail  separation  values  taken  from 
the  same  study  are  shown  in  Table  1.2.  Values  of  equivalent  arrival  in-trail  separation  have 
been  added  to  show  that  even  with  the  projected  fixed  distance  separation  standards  large 
in-trail  separation  times  remain  tor  some  operational  scenarios  (i.e.,  a  small  aircraft  operating 
behind  a  heavy  heavy  jet)  Moreover,  one  must  recognize  that  challenges  remain  to  be  solved 
in  terms  of  the  technology  available  to  support  some  of  these  control  actions.  For  example,  the 
wake  vortex  separation  criteria  envisioned  in  the  MITRE  study  would  necessitate  of  a 
completely  revised  scheme  in  how  aircraft  control  systems  are  designed  or  through 
sophisticated  approach  sequences  which  might  not  only  include  curved  paths  but  also 
real-time  advisories  of  the  progress  of  preceeding  landing  aircraft.  The  first  point  is  probably 
the  most  debatable  since  the  vortex  wake  generated  by  a  heavy  aircraft  requires  a  significant 
amount  of  time  to  decay  to  harmless  levels  and  the  reduction  of  the  current  separation  criteria 
could  result  in  substantial  increases  in  vortex  wake  penetration  [Rossow  and  Tinglin,  1988j. 
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The  second  approach  seems  more  promising  with  the  introduction  of  the  microwave  landing 
system  (MLS)  coupled  with  area  navigation  procedures  (RNAV)  and  its  well  publicized 
multipath  approach  capabilities  although  more  flight  testing  is  necessary  to  validate  potential 
MLS-RNAV  procedures  for  heavy  aircraft  operations  in  the  "near"  terminal  control  area 
[Branstetter  et  al,  1988).  The  automation  of  the  future  Air  Traffic  Control  (ATC)  activities  in  the 
terminal  control  aerea  seem  to  promise  a  small  reduction  in  the  interarrival  time  and  runway 
occupancy  time  buffers  currently  applied  under  manual  control  ATC  conditions  [Swedish  , 
1979].  If  these  reductions  were  to  take  place  in  the  near  future,  the  runway  subsystem  could 
in  fact  become  the  "bottleneck"  of  the  airport  system  and  thus  airport  capacity  could  then  be 
dictated  by  the  runway  occupancy  times  achieved  during  realistic  operations  Under  the 
simple  assumption  of  considering  an  average  ROT  value  in  the  future,  if  the  2/4/5  n.m. 
distance  in-trail  separation  rule  were  implemented  in  the  future  then  the  average  runway 
occupancy  time  (ROT)  necessary  to  accomodate  this  level  of  arrival  rates  would  be  39 
seconds  with  a  1.4%  intervention  rate  [Swedish,  1979],  According  to  observations  made  by 
Koenig  (Koenig,  19/8]  at  six  major  airports  this  average  is  below  the  42-51  second  average 
ROT  experienced  by  short/medium  size  transport-type  aircraft  and  the  51-58  second  average 
ROT  observed  in  heavy  jet  VFR  operations.  It  is  interesting  to  observe  that  these  averages 
have  not  significantly  changed  for  more  than  a  decade  suggesting  that  the  current  runway 
subsystem  needs  to  be  modified  if  future  gains  in  average  ROT  are  to  be  achieved  It  seems 
desirable  to  investigate  high-speed  turnoff  alternatives  that  could  reduce  the  present  average 
ROT  values  by  4  to  10  seconds  (i.e.,  depending  upon  the  operating  aircraft  mix)  to  balance  the 
expected  future  interarrival  separation  standards  with  the  average  runway  occupancy  times 
and  thus  improving  runway  capacity  under  mixed  operations. 

Parallel  research  studies  are  currently  being  conducted  to  address  the  air  space  and 
runway  capacity  issues  in  order  to  improve  the  level  of  service  of  existing  and  future  facilities. 
It  is  important  to  emphasize  that  a  reduction  in  the  in-trail  separation  rules  will  have  to  be 
followed  by  a  corresponding  reduction  in  ROT  times  if  significant  improvements  in  the  airport 
capacity  are  to  be  achieved.  The  use  of  highly  advanced  avionics  coupled  with  the  installation 
of  modern  navigation  facilities  (i.e.,  microwave  landing  system,  four-dimensional  navigation, 
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etc.)  and  the  strategic  location  of  turnoffs  seem  to  promise  some  reductions  in  the  in-t rail 
separation,  ROT  and  their  corresponding  variabilities  to  warrant  changes  in  runway  capacity. 

Among  the  several  alternatives  to  improve  airport  capacity  one  being  addressed  in  this 
research  is  the  use  of  high-speed  turnoffs  in  airport  runways  to  reduce  the  aircraft  runway 
occopancy  times.  It  should  be  clear  that  although  some  of  the  major  airport  facilities  have 
adopted  the  30-degree  high-speed  FAA  turnoff  standard  and  newly  proposed  spiral  designs  it 
seems  possible  to  improve  the  level  of  service  of  a  runway  facility  (i.e.,  decreasing  the  runway 
occupancy  time  or  increasing  the  runway  acceptance  rate)  by  tailoring  the  turnoff  geometry 
and  location  to  an  existing  or  forecasted  population  of  aircraft.  That  is  the  population  of  aircraft 
operating  in  a  particular  airfield  should  dictate  not  orlly  the  location  but  also  the  characteristic 
turnoff  geometries  associated  with  that  runway. 


1.1  Previous  Turnoff  Geometry  Research 


Past  studies  by  HoronjefT  et  al  [Horonjeff  et  al,  1958, 1959,  and  1960]  recognized  the  critical 
relationship  between  turnoff  location  and  turnoff  geometry  and  developed  a  mathematical 
model  to  locate  exit  taxiways  for  a  limited  number  of  scenarios  (i.e.,  two  exit  taxiway  speeds 
and  a  reduced  aircraft  population).  The  results  of  this  model  concluded  that  the  optimum 
location  of  runway  turnoffs  is  quite  sensitive  to  aircraft  population,  number  of  exits,  and  exit 
speeds.  The  same  model  used  external  atmospheric  corrections  to  modify  the  baseline  results 
due  to  meteorological  and  geographical  conditions.  However,  only  two  exit  speeds  (i.e.,  40  and 
60  mph)  and  a  limited  number  of  aircraft  populations  were  investigated  thus  making  the  model 
of  limited  use.  Furthermore,  since  the  aircraft  populations  used  comprised  "old'  aircraft  by 
current  standards  the  results  need  revision.  The  pioneering  effort  of  the  Horonjeff  team, 
however,  generated  a  good  amount  of  information  regarding  the  cornering  capabilities  of 
aircraft  and  also  obtained  data  on  several  lighting  schemes  to  help  pilots  negotiate  these 
turnoffs  under  adverse  weather  conditions.  The  Horonjeff  team  performed  extensive 
experiments  to  find  the  acceptable  turning  radius  at  a  given  exit  speed.  The  results  suggested 
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Table  1.2  Arrival  In-Trall  Separation  Criteria. 


Lead 

Aircraft 

Small 

Current  Values 
(nautical  miles) 
(seconds) 

Large 

Future  Goal  (1996) 

(nautical  miles) 

(seconds) 

Heavy 

Small 

Large 

Heavy 

Small 

2.5 

2.5 

2.5 

2.0 

2.0 

2.0 

82 

64 

60 

65 

51 

48 

Large 

4.0 

2.5 

2.5 

3.0 

2.0 

2.0 

131 

64 

60 

98 

51 

48 

Heavy 

6.0 

5.0 

4.0 

5.0 

40 

3.0 

196 

129 

96 

163 

103 

72 

Assumed  speeds  for  three  aircraft  classes: 
Small  -  1 10  knots 
Large  -  140  knots 
Heavy  -  150  knots 


two  centered  curves  for  the  turnoff  geometry  with  specifications  shown  in  Table  1.3  and  Figure 
1.1  IHoronjeff  et.al.,  1958], 

In  1970,  FAA  (1970,  AC  150/5335-1A)  made  standards  of  high  speed  exits,  angled  exits  with 
30°  and  45  0  .  The  new  FAA  standard  employed  Horonjeff's  suggested  ending  radius  of  1800 
feet  which  equates  to  an  exit  speed  of  60  MPH.  Another  standard  adopted  for  smaller  aircraft 
was  a  45  degree  angle  exit  with  a  radius  of  curvature  of  800  ft.  which  could  well  serve  aircraft 
exiting  at  40  MPH  according  to  Horonjeff's  findings.  One  important  missing  item  from  these 
high  speed  geometries  was  the  easement  or  transition  curve  necessary  to  provide  passenger 
comfort  while  executing  the  initial  portion  of  the  turnoff  maneuver.  Interestingly  enough  the 
International  Civil  Aviation  Organization  (ICAO)  adopted  the  easement  curves  as  well  as  the 
second  radius  of  curvature  suggested  by  the  1958  UC  Berkeley  team  (ICAO,  1977). 

Schoen  et.  al.  [Schoen  et  al.,  1985]  investigated  the  turnoff  trajectory  of  high  speed  taxiing 
aircraft  in  an  isolated  basis.  The  resulting  ape  of  the  aircraft  turnoff  was  a  variable  curvature 
geometry  with  a  continuously  decreasing  radius  of  curvature.  The  end  result  of  this  research 
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Table  1.3  Turnoff  Geometry  Specification  Table. 


Velocity 

(MPH) 

ul> 

<fl) 

ith 

40 

1724 

189 

821 

54 

2936 

236 

1282 

64 

3138 

283 

1846 

A’,  radius  ofrr, 
/ ,  length  ot  et 
/»’  /odw  f  of  ft 


?0  V)  40  M)  U)  f  0 

0 peed,  in  miles  per  hot/' 


Figure  1.1  Turnoff  Geometry  Specification  Graph  [Horonjeff  et.  al.,  1958]. 
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was  a  computer  program  to  calculate  the  (x,  y)  coordinates  of  the  geometry,  considcnng  exit 
speed  and  aircraft  turning  ability.  The  findings  of  this  research  suggested  that  aircraft 
rotational  inertia  played  an  important  factor  in  dictating  the  initial  trajectory  of  the  turnoff 
maneuver.  This  research  also  showed  that  ROT  values  of  30  seconds  are  possible  at  the 
expense  of  large  turning  radius  and  extremely  high  exit  speeds  (e  g.,  110  MPH  for  a  Boeing 
747).  Very  high-speed  turnoff  results  should,  however,  be  treated  cautiosly  since  at  such  high 
speeds  the  controllability  of  aircraft  on  the  ground  could  become  a  serious  operational 
deterrent. 

The  most  recent  research  of  turnoff  geometry  was  conducted  by  Aviation  Department 
staffs  of  Dade  County,  Florida  (  Carr  et.  al.  .  1980,  Witteveen,  1987,  and  Haury,  1987  ).  They 
tested  various  types  of  geometry,  lighting,  and  marking  in  an  L1011  flight  simulator.  Figure 
1.4  is  a  proposed  wide  entrance  geometry  which  shows  its  wider  throat  than  that  of  a  standard 
FAA  high-speed  exit  geometry.  Figure  1.5  shows  the  difference  of  the  standard  geometry  and 
suggested  geometry.  The  shaded  area  is  the  standard  geometry.  The  new  geometry  was 
implemented  in  Miami  International  Airport,  Baltimore-Washington  International  Airport,  and 
Orlando  International  Airport,  and  is  expected  to  be  constructed  in  Cincinnati  International 
Airport  and  the  new  Denver  International  Airport  (Witteveen,  1987  and  Haury,  1987],  The  wide 
entrance  throat  of  this  geometry  is  appealing  in  situations  where  lateral  spacing  restrictions 
between  the  runway  and  the  nearest  parallel  taxiway  are  severe  (i.e.,  less  than  600  ft  ). 
However,  the  ending  radius  of  curvature  of  only  800  ft.  might  be  a  limiting  factor  in  the 
operational  capabilities  of  this  exit  to  handle  large  aircraft  above  50  knots  in  a  routine  basis. 
The  FAA  is  currently  engaged  in  evaluating  this  geometry  in  the  Boeing  727-100  simulator  and 
aircraft. 


1.2  Mathematical  Models  for  Optimal  Turnoff  Locations 


The  earliest  effort  to  make  a  model  for  the  optimal  runway  exit  locating  problem  is  also 
found  in  the  pioneering  work  of  Robert  Horonjeff  (Horonjeff  et.  al.,  1959,  1960).  The  objective 
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of  their  model  was  to  find  exit  locations  that  maximize  the  landing  acceptance  rate  of  a  runway 
in  the  saturated  situation,  assuming  the  aircraft  arrival  pattern  followed  a  fixed  time  or  fixed 
distance  separation  criterion.  The  saturated  situation  means  that  aircraft  try  to  land 
continuously  with  a  separation  rule.  The  acceptance  rate  Is  determined  by  : 


E(AC) 


1  1 
E(6)  *  (1+q) 


{1.1} 


where, 

E(AC)  =  the  expected  acceptance  rate  {  aircraft/hr ), 

E(6)  =  the  expected  interarrival  time  (  hr/aircraft ), 
q  =  weighted  average  of 'wave  off'  prooabilities. 

Wave  off  occurs  if  the  previous  aircraft  remains  on  the  runway  when  the  next  aircraft 
reaches  the  runway  threshold.  The  expected  acceptance  rate  can  be  maximized  by 
minimizing  the  weighted  average  of  wave  off  probabilities  which  expressed  by  : 


n 

Q  =  £p/  •  9/  {12} 

i 

where, 

p,  =  the  proportion  of  aircraft  type  i, 

q,  =  the  wave  off  probability  of  aircraft  type  i, 
n  =  the  number  of  aircraft  types. 

p,  should  be  provided  to  the  model,  and  q,  is  calculated  by  : 
q,  =  Pr  {T,  >  6}  {13} 

where, 

T,  =  runway  occupancy  time  of  aircraft  type  i. 

Since  T,  is  a  function  of  exit  locations  {  0,,  Dt,  ,  Dm),  q  can  be  expressed  as  a  function 
of  exit  locations,  and  thus,  by  calculus,  the  optimal  exit  locations,  which  minimize  q,  are  found 
The  equation  for  T,  involves  bivariate  random  variables,  (dt,  t,)  which  are  the  mean  dish  nee 
and  time  for  aircraft  of  type  i  to  decelerate  to  the  predetermined  exit  speed,  respectively.  The 
differential  equations  for  optimal  exit  locations  are  not  simple,  and  can  not  be  solved 
analytically.  Hence,  finding  the  optimal  locations  requires  a  numerical  computation  algorithm 
which  consumes  a  lot  of  computation  time. 
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The  joint  distributions  of"(d,,  t,)  of  every  aircraft  type  are  another  input  data  (or  the  model 
The  values  of  di  and  f/  vary  according  to  the  operational  factors  such  as  the  design  exit  speed 
and  the  landing  weight,  and  environmental  factors  such  as  runway  surlaco  conditions,  altitude; 
of  airport,  and  temperature,  even  though  we  consider  only  one  aircraft  type.  Hence  the 
parameters  of  the  joint  distribution  of  d,  and  f,  should  be  estimated  again  if  an  influencing 
factor  is  changed  thus  posing  a  great  computational  and  labor  intensive  challenge. 

In  1974,  Daellenbach  [1974]  developed  a  dynamic  programming  model  which  is  equivalent 
to  the  Horonjeff's  approach  with  some  extensions.  Horonjeff's  model  imposes  a  strict 
assumption  on  the  aircraft  arrival  pattern.  Daellenbach  released  the  assumption,  and 
permitted  a  generalized  arrival  pattern.  He  showed  his  model  to  be  more  efficient 
computationally  and  more  flexible  for  modelling  than  Horonjeff's  model.  Daellenbach's  model, 
however,  also  requires  the  joint  distributions  of  (dit  t!)  as  input.  The  data  for  estimating  the 
parameters  of  the  joint  distributions  are  difficult  to  collect  and  almost  impossible  when  the 
influencing  factors  vary. 

In  the  same  year,  Joline  [1974]  developed  another  dynamic  programming  model  to  find 
the  optimal  number  of  exits  and  their  locations  with  respect  to  the  combined  objective  function 
of  ROT  and  exit  construction  cost.  He  incorporated  the  ROT  gain  and  the  exit  construction  cost 
into  an  objective  function  by  equating  1  second  gain  in  ROT  with  $100,000  in  construction  cost 
While  Horonjeff's  model  and  Daellenbach's  model  require  the  joint  distributions  of  (d,,  t,)  for 
each  aircraft  type,  Joline's  model  needs  only  an  univariate  distribution  of  'ideal  exit  location' 
for  a  mixed  aircraft  population.  Joline  classified  aircraft  into  three  categories  based  on  the 
aircraft  size,  and  found  the  distributions  of  ideal  exit  locations  for  these  three  aircraft  classes 
based  on  the  observations  of  aircraft  landing  operations  in  Chicago  O'Hare  Airport.  The  ideal 
exit  location  distribution  for  entire  aircraft  population  is  found  by  combining  the  three 
distributions  according  to  the  proportions  of  the  three  aircraft  classes.  As  mentioned  earlier, 
there  are  several  factors  influencing  the  aircraft  landing  distance  such  as  the  design  exit 
speed,  landing  weight,  etc..  Joline's  model,  like  the  previous  models,  makes  the  effects  of 
these  influencing  factors  hard  to  incorporate. 


Introduction 


12 


ROT  consists  of  the  time  from  the  runway  threshold  to  the  exit  location  and  the  time  from 
the  beginning  of  the  turn  to  the  clearance  of  the  runway.  The  second  term  obviously  varies 
according  to  the  turnoff  trajectory  and  its  magnitude  varies  according  to  the  desired  turnoff 
exit  angle.  None  of  three  models  above,  however,  takes  into  account  the  ROT  variation  due 
to  the  change  of  the  turnoff  trajectory.  The  turnoff  trajectory  also  varies  according  to  design 
exit  speed,  aircraft  turning  ability,  runway  surface  conditions,  etc..  Therefore  an  attempt  is 
made  in  this  report  to  bridge  the  gap  between  practical  and  theoretical  models. 

The  three  models  above  implicitly  assume  that  an  aircraft  type  can  use  more  than  one 
exit  for  turnoff  with  different  ROT  and  exiting  probability.  If  we  want  to  decide  separation  times 
between  the  landing  aircraft  based  on  the  ROT  of  the  aircraft,  it  is  desirable  to  assign  an 
aircraft  to  an  exit  with  high  exiting  probability,  say  99%.  Thereby  the  variation  of  ROT  of  an 
aircraft  can  be  reduced.  This  situation  is  expected  to  occur  in  the  near  future  with  an 
improvements  to  aircraft  traffic  control  (  ATC  )  systems  and  better  crew  situational  awareness 
provided  by  enhanced  on-board  ground  navigation  avionic  systems.  The  purpose  of  this  new 
research  effort  is  to  develop  a  more  general  model  that  will  be  able  to  predict  that  location 
of  runway  turnoffs  in  a  myriad  of  scenarios  are  where  changes  to  the  aircraft  population  and 
airfield  conditions  can  be  easily  defined  by  the  user  thus  making  it  possible  to  be  used  under 
practical  airport  conditions. 

1.3  Approach  and  Scope  of  the  UCTRIVPI  Research 

The  University  Center  for  Transportation  Research  (UCTR)  at  Virginia  Tech  University  (VPI 
&  SU)  has  been  asked  by  the  FAA  and  NASA  Langley  to  investigate  the  feasibility  of 
implementing  high-speed  turnoffs  at  major  airports  to  reduce  the  runway  occupancy  times. 
The  goal  of  this  research  project  is  to  develop  a  user-friendly  computer  simulation  model  to 
estimate  the  optimal  location  of  high-speed  turnoffs  at  an  airport  facility  in  order  to  reduce  the 
runway  occupancy  time  (ROT).  The  mode,  .^corporates  environmental  factors  such  as  airfield 
elevation,  runway  configuration,  weather  conditions,  etc.  and  operational  factors  such  as 


aircraft  mix  (i.e.,  terminal  operation  aircraft  types)  and  aircraft  piloting  technique  (i.e.,  aircraft 
touchdown  speed  dispersion  values)  to  determine  the  potential  location  of  high-speed  turnoffs. 
Factors  that  in  past  studies  either  have  been  neglected  or  treated  as  externalities  to  the 
models. 

A  problem  scenario  addressed  in  this  research  can  be  better  understood  with  the  help 
of  Fig.  1.6  which  depicts  a  histogram  of  landing  distances  for  a  typical  airport  facility.  It  is 
observed  that  a  grouping  phenomena  of  the  potential  location  of  the  runway  turnoffs  takes 
place  at  discrete  distances  from  the  active  runway  threshold  as  each  aircraft  type  has  unique 
landing  rollout  performance  characteristics  coupled  with  some  inherent  variability  due  to 
dissimilar  piloting  techniques.  The  point  to  be  stressed  here  is  that  by  carefully  locating  the 
runway  turnoffs  one  could,  at  least  in  principle,  reduce  the  ROT  time  for  a  given  aircraft 
population  and  airfield  conditions  below  a  desired  level  (e.g.,  40  seconds).  Fig.  1.7  shows  an 
example  of  projected  turnoff  locations  for  a  hypothetical  runway. 

A  second  argument  to  this  already  complex  problem  is  the  fact  that  due  to  the  unique 
landing  characteristics  of  each  aircraft  it  is  also  possible  to  determine  the  most  efficient  turnoff 
geometry  for  the  same  aircraft  population.  Fig.  1.8  depicts  five  different  turnoff  geometries  for 
dissimilar  aircraft  representative  of  the  existing  Terminal  Instrument  Procedure  categories 
(TERPS)  using  a  probabilistic  model  developed  by  the  McDonnell  Douglas  Corporation  for 
NASA  (Schoen  et.  al.,  1985],  These  results  consider  the  aircraft  landing  gear  friction 
characteristics  to  be  the  only  source  of  aircraft  directional  control  on  the  ground,  an 
approximation  well  suited  for  medium  to  slow  landing  rollout  speeds  (e  g.,  30-60  knots)  where 
the  aerodynamic  efficiency  of  the  aircraft  primary  control  surfaces  is  negligible.  However,  the 
model  is  very  restrictive  as  it  only  optimizes  the  geometry  of  a  turnoff  for  a  single  aircraft  with 
no  consideration  for  specific  environmental  nor  operational  factors  of  the  airfield.  Again,  the 
issue  of  considering  the  aircraft  mix  as  part  of  the  solution  of  the  problem  arises  when  one 
considers  that  the  selection  of  a  critical  aircraft  might  not  warrant  the  overall  'best”  alternative 
if  that  critical  aircraft  seldom  operates  the  facility.  In  other  words,  if  an  optimal  location  is 
sought  for  a  large  population  of  aircraft  it  might  be  desirable  to  penalize  the  critical  aircraft 
in  terms  of  their  ROT  parameter  if  the  rest  of  the  population  is  large  and  can  be 
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Figure  1.6  Hypothetical  Aircraft  Landing  Distance  Distribution. 


accommodated  in  a  particular  geometry  that  reduces  the  overall  ROT  (i.e.,  average  ROT)  of 
the  facility  in  question.  This  clearly  demonstrates  that  an  optimal  solution  must  heavily  depend 
upon  the  aircraft  population  mix  operating  in  the  airport  environment. 

From  this  last  point  an  interesting  question  emerges  regarding  the  applicability  of  an 
"optimal"  design  for  current  and  future  conditions.  One  might  say  that  if  a  runway  is  at  the 
planning  stage,  the  planner  should  attempt  to  forecast  the  future  aircraft  population  of  its 
environment  with  the  aid  of  airlines  and  aviation  authorities,  while  if  only  improvements  are 
sought  (i.e.,  the  runway  is  already  in  place)  the  use  of  a  current  aircraft  population  constitutes 
a  better  choice.  The  bottom  line  seems  to  be  that  the  model  should  be  flexible  enough  to  allow 
the  analyst  to  execute  environmental  and  operational  changes  with  a  minimum  of  effort. 
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The  scope  of  the  model  should  be  viewed  not  as  an  isolated  effort  to  address  all  the 
problems  regarding  runway  operations,  but  as  a  novel  approach  to  solve  some  of  the  imposed 
demands  generated  with  the  growth  of  air  traffic  operations  and  the  need  for  using  the  airport 
infrastructure  available  more  efficiently.  The  effort  presented  here  should  be  integrated  and 
coordinated  with  future  and  present  complementary  efforts.  The  main  emphasis  from  the 
user's  point  of  view  will  be  in  the  general  aviation  area  but  the  model  avoids  preferences 
towards  a  selected  group  of  aircraft  and  can  be  expanded  to  suit  any  operating  aircraft  mix. 

The  model,  once  developed,  should  be  extended  and  integrated  to  the  airport  capacity 
analysis  level.  At  this  stage  the  interactions  between  runway,  taxiway  ,  and  apron  should  be 
studied  to  perceive  potential  problems  of  integrating  variable  geometry  turnoffs  in  the  airfield 
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Figure  1.8  Sample  Turnoff  Geometries  for  Different  Aircraft  Categories. 


environment.  Following  the  implications  of  these  results  a  cost  effectiveness  analysis  would 
have  to  be  performed  in  order  to  perceive  the  gains  in  ROT  time  translated  into  benefits  to  the 
user,  and  operators.  Finally,  the  results  will  have  to  be  tested  for  pilot's  acceptance  in  terms 
of  procedures  and  workload.  These  studies  could  be  carried  out  in  flight  simulators  and  the 
aircraft  themselves. 


1.4  Model  Overview 


The  Runway  Exit  Design  Interactive  Model  (REDIM)  developed  in  this  research  effort 
incorporates  several  specific  airfield  variables  that  affect  the  landing  performance  of  the 
aircraft  as  well  as  other  important  operational  constraints  (e  g.,  aircraft  mix)  that  have  a  direct 
impact  on  the  selection  of  the  turnoff  location  and  their  geometry.  A  simplified  flow  diagram 
of  the  model  proposed  is  shown  in  Fig.  1.9.  It  can  be  seen  that  five  modules  comprise  the 
program:  1)  an  input  module  ,  2)  a  dynamic  simulation  module  to  estimate  the  ROT  times 
for  individual  aircraft,  3)  a  selection/optimization  algorithm  to  determine  candidate  turnoff 
locations  to  comply  with  a  desired  reliability  threshold  value  and  4)  an  aircraft  data  module 
containing  relevant  aircraft  performance  and  geometric  parameters  (also  named  master  file 
in  this  research  and  5)  an  output  module  to  shown  graphically  and  in  tabular  form  the 
suggested  runway  turnoff  configuration  and  display  some  measures  of  effectiveness  desired 
by  the  analyst.  Fig.  1.10  shows  a  typical  sequence  of  events  occuring  during  the  flare  transition 
and  landing  rollout  that  will  be  modeled  explicitly  in  the  dynamic  portion  of  this  program.  The 
program  contains  a  library  of  geometric  and  operational  aircraft  characteristics  to  allow  the 
analyst  to  choose  from  a  wide  selection  of  operational  airport  scenarios.  Obviously,  the  user 
is  also  capable  of  editing  his/her  own  aircraft  data  if  desired  through  simple  steps  in  the  input 
module. 

The  program  considers  three  broad  types  of  analyses:  1)  evaluation  of  an  existing  runway, 
2)  redesign  of  an  existing  runway  and  3)  design  of  a  new  runway  facility.  In  the  evaluation 
mode  REDIM  estimates  several  measures  of  effectiveness  indicative  of  the  operational 
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capabilities  of  an  existing  runway  facility.  In  this  mode  the  user  inputs  the  number,  type  and 
location  of  existing  turnoffs  as  well  as  the  relevant  aircraft  population  data  and  the  model 
predicts  the  average  runway  occupancy  time  (WAROT),  the  particular  exit(s)  that  an  aircraft 
can  take  ,  and  the  probability  of  each  aircraft  i  taking  the  assigned  exit(s).  Another  potential 
use  of  this  mode  is  to  serve  as  a  benchmark  to  perform  valid  comparisons  between  different 
runway  proposals.  This  way  the  analyst  can  perceive  the  operational  gains  of  various 
modification  alternatives. 


The  second  mode  of  operation  deals  with  the  redesign  of  a  runway  facility.  In  this  scenario 
it  is  expected  that  the  user  might  want  to  explore  the  possibility  of  adding  new  high-speed 
turnoffs  to  an  existing  facility  and  examine  their  impact  in  the  operational  efficiency  of  the 
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facility.  Inputs  In  this  mode  are  the  number  and  type  of  existing  turnoffs,  their  locations,  the 
number  of  new  turnoffs  to  be  constructed  and  a  reliability  parameter.  The  outputs  are  the 
location  and  geometry  of  each  new  turnoff,  the  weighted  average  runway  occupancy  time,  and 
an  aircraft  assignment  table  containing  individual  runway  occupancy  times  for  the  desired 
reliability  factor  specified  by  the  user. 

In  the  third  mode  of  operation  REDIM  estimates  the  optimal  location  of  runway  turnoffs 
and  their  corresponding  geometries.  An  assignment  table  is  given  to  the  user  indicating  the 
turnofffs)  associated  with  each  aircraft  and  their  individual  runway  occupancy  times.  The 
weighted  average  runway  occupancy  time  is  also  estimated  as  a  global  runway  operational 
parameter  and  sensitivity  studies  can  be  easily  be  conducted  by  changing  the  number  of 
turnoffs  allocated  to  a  specific  runway.  Inputs  by  the  user  in  this  mode  are  the  number  of  exits 
to  be  constructed  and  the  desired  exit  reliability  parameter.  More  detailed  descriptions  of 
these  three  modes  of  operation  will  be  given  in  the  remaining  chapters  of  this  report. 

REDIM  blends  the  principles  of  continuous  simulation  with  those  of  mathematical 
optimization  to  find  the  best  turnoff  locations  and  corresponding  turnoff  geometries  for  a 
myriad  of  possibilities.  The  program  was  designed  to  be  interactive  and  a  great  effort  was 
made  to  reduce  the  number  of  inputs  expected  from  the  user.  A  large  aircraft  data  base  is 
included  to  simplify  the  analyst  input  task  but  flexibility  is  also  built-in  to  allow  future  aircraft 
additions.  The  overall  effort  was  to  make  the  program  interactive  and  easy  to  use  Interactive 
input  and  output  menus  are  easy  to  follow  providing  the  user  with  graphic  results  on  the 
screen  or  a  regular  line  printer. 
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Figure  1.10  Aircraft  Landing  Events  Modeled  in  the  Analysis. 


2.0  Input  Module 


The  Input  Module  comprises  a  series  of  Interactive  screens  that  allow  the  user  to  input 
and  edit  data  necessary  for  the  analysis  portion  of  the  program  (i.e.,  Dynamic  and 
Optimization  Modules).  This  module  is  controlled  by  menus  or  key-stroke  commands  such  as 
"Esc”  key.  A  flow  chart  depicting  the  sequence  of  events  comprising  the  Input  Module  is  shown 
in  Fig.  2.1.  As  can  be  observed,  the  "Main  Menu"  placed  at  the  top  level  of  the  flow  chart  has 
seven  modes:  1)  "Start  for  a  New  Problem,"  2)  "Edit  Data,"  3)  "Begin  Analysis,"  4)  "Edit  Master 
File,"  5)  "Go  To  Output  Module,"  6)  "Help"  and  7)  "Quit."  The  details  of  these  modes  are  given 
in  Section  2.4  and  shown  schematically  in  Fig.  2.1. 

Input  Data  is  classified  into  six  broad  categories:  1)  analysis  type  and  related  data,  2) 
aircraft  mix  and  characteristics  data,  3)  airport  operational  data,  4)  airport  environmental  data, 
5)  runway  gradients,  and  6)  weather  and  exit  speeds.  All  of  these  are  necessary  for  the 
analysis,  and  should  be  saved  in  a  'working  data  Tile'  specified  by  the  user  with  an  arbitrary 
name.  For  the  convenience  of  the  user,  all  the  aircraft  characteristics  are  kept  in  a  master 
data  file  named  "MACT.DAT"  and  are  transferred  to  a  working  data  file  automatically  if 
necessary. 


2.1  Data  Classification 


There  are  three  kinds  of  data  for  analysis:  1)  input  data,  2)  constant  data,  and  3)  calculated 
data.  Among  these  kinds  of  data,  constant  data  and  calculated  data  are  determined  in  the 
Dynamic  Module  detailed  in  the  next  Chapter.  Input  data  is  provided  by  the  user  via  the  Input 
Module  and  its  user-friendly  screens.  The  input  data  is  classified  into  six  categories  as 
mentioned  previously.  The  following  paragraphs  define  the  categories  in  more  detail. 

Analysis  Type  and  Related  Data 
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Figure  2.1  Input  Module  Functional  Flow  Chart. 


The  program  provides  Ihe  user  with  ihroe  choices  for  the  lype  ol  analysis  lo  Im  pcilormed 
For  each  type  of  analysis,  (lieie  are  some  accompanying  data  which  vary  depending  on  the 
user's  choice  A  morn  detailed  explanation  of  this  is  provided  m  Section  2  2  devoted  lo  Input 
/  Output  relationships. 

Aircraft  Mix  and  Characteristics  Data 

In  this  category,  the  percentages  of  the  aircraft  comprising  the  airport  population  mix  and 
aircraft  geometric  characteristics  used  in  the  program  are  included. 

Airport  Operational  Data 

In  this  category,  the  free  roll  time  between  the  touchdown  and  the  beginning  of  braking,  the 
free  roll  time  between  the  end  of  braking  and  the  beginning  of  turn  off,  taxiing  speed,  and  their 
standard  deviations  are  included.  A  safety  factor  for  the  impending  skidding  condition  is  also 
pail  of  this  screen. 

Airport  Environmental  Data 

The  following  parameters  are  included  in  this  category:  wind  speed,  wind  direction,  airport 
elevation,  airport  temperature,  runway  orientation,  visual  range,  runway  width,  and  distance 
to  the  nearest  taxiway 

Runway  Gradient 

In  this  category,  runway  length,  and  the  effective  gradient  for  every  one  tenth  of  runway  are 
included 

Weather  and  Exit  Speeds 

The  relative  frequency  of  dry  and  wet  runway  surface  conditions  are  included  in  this  category. 
The  desired  exit  speeds  of  each  aircraft  category  on  each  surface  condition  are  also  included 
here 
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2.2  InputlOutput  Relationship 


As  stated  earlier,  the  user  can  select  one  of  three  types  of  analysis  :  t)  evaluation  of  an 
existing  runway  system,  2)  improvement  of  an  existing  runway  system,  and  3)  design  of  a  new 
runway  system. 

The  'evaluation  option  requires  1)  the  reliability  parameter,  2)  the  number  of  existing 
exits,  and  3)  the  types  of  existing  exits  as  input  data.  This  produces  as  output  1)  the  aircraft 
distribution  to  the  existing  exits  with  the  corresponding  ROT  and  cumulative  exiting 
probabilities  and  2)  the  weighted  average  ROT  as  results  of  analysis.  When  the  average  ROT 
is  calculated  for  all  the  aircraft  population,  only  one  ROT  is  considered  for  each  aircraft.  For 
example,  if  aircraft  I  is  able  to  take  exit  k  and  k  +  1  whose  corresponding  ROT's  and 
cumulative  probabilities  are  tk,  f*  + ,,  and  p*p*  +  1,  respectively.  If  p*  is  less  than  the  reliability 
specified  by  the  user,  and  is  greater  than  the  desired  reliability,  then  only  f. . .  is 
considered  as  the  representative  ROT  for  aircraft  i. 

The  'improvement'  option  assumes  that  a  few  exits  would  be  add^d  to  an  existing  runway. 
This  analysis  requires  1)  the  reliability  parameter,  2)  the  number  of  existing  exits,  3)  the  types 
of  existing  exits,  and  4)  the  number  of  new  exits  which  will  be  constructed.  The  results  are 
1)  optimal  exit  locations,  2)  aircraft  assignment  to  the  existing  and  new  exits,  3)  the  weighted 
average  ROT  which  is  minimized  by  the  optimal  exit  locations,  and  4)  turnoff  geometries  of  the 
exits. 

The  'design'  option  assumes  a  hypothetical  situation  with  no  exits  on  the  runway.  The 
reliability  parameter  and  the  number  of  new  exits  are  inputs  for  this  type  of  analysis  The 
results  are  similar  in  nature  to  those  of  the  'improvement'  option  with  the  only  difference  being 
is  that  this  new  option  takes  into  account  the  new  exits  for  aircraft  assignment,  while  the 
'improvement'  option  considers  both  new  and  existing  exits.  I/O  relationships  for  each 
analysis  type  are  shown  in  Figs  2  2  to  2  4. 
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Design  a  New  Runway  Facility  (Design  Mode) 


Figure  2.4  I/O  Relationship  for  "Design'  Analysis. 


2.3  Data  Input  Method 


In  the  Input  Module,  there  are  three  different  input  methods  used:  1)  menu  input,  2)  line 
input,  and  3)  table  input.  Menu  input  arises  when  the  user  selects  his  choice  among  the  list 
displayed  on  the  screen  using  the  arrow  keys  and  enter  key.  The  flow  in  the  module  is 
controlled  by  the  menu  input  method.  The  main  menu,  edit  menu  for  working  data  file,  edit 
menu  for  master  data  file,  selection  of  a  analysis  type,  etc.  are  the  examples  of  the  menu  input 
method.  Line  input  occurs  when  the  user  puts  a  numerical  value  like  runway  length  or  a  string 
datum  like  a  data  file  name  at  the  position  specified  on  a  screen.  The  user  inputs  file  names 
(  data  and/or  output  file),  the  number  of  exits,  the  locations  of  exits,  etc.  using  this  method. 
Table  input  is  similar  to  line  input.  However,  table  input  is  used  in  order  to  get  several 
numerical  data  on  the  same  screen,  while  line  input  is  used  in  order  to  get  one  numerical  or 
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string  datum  on  a  line.  By  the  table  input  method  the  user  inputs  aircraft  mix  data,  aircraft 
characteristics  data,  airport  environmental  data,  etc. 


2A  Procedures  in  Input  Module 


The  first  screen  presented  to  the  user  after  the  title  screen  is  the  main  menu  The  main 
menu  screen  gives  the  user  a  list  of  seven  choices:  1)  Start  a  New  Problem,  2)  Edit  Data,  3) 
Begin  Analysis,  4)  Edit  Master  File,  5)  Go  to  Output  Module,  6)  Help,  and  7)  Quit  as  shown  in 
the  Fig.  2.5.  The  functions  of  choices  3),  5),  and  7)  are  used  to  transfer  the  flow  control  to  other 
program  modules.  The  purpose  of  choice  6)  is  to  give  a  brief  explanation  of  the  screen  choices 
to  the  user.  Choices  1,2  and  4  invoke  the  procedures  belonging  to  the  Input  Module. 


2.4.1  Starting  a  New  Problem 


In  this  mode,  all  the  data  necessary  for  the  analysis  should  be  provided  by  the  user  Once 
the  user  enters  this  mode  a  complete  set  of  values  is  expected  before  completing  the  entire 
input  process.  The  first  set  of  data  which  the  user  needs  to  specify  is  'Analysis  Type  and 
Related  Data"  {  refer  to  the  Section  2.1  ).  The  type  of  analysis  is  selected  through  the  menu 
input  method.  Fig.  2.6  shows  the  screen  for  analysis  type  selection,  which  is  followed  by  the 
related  data  screen  which  might  vary  depending  upon  the  analysis  type  selected.  Fig.  2.7  is 
an  example  of  such  a  screen.  The  second  set  of  the  data  is  "Aircraft  Mix  and  Characteristics." 
The  aircraft  mix  sc  -een  shows  the  names  of  the  aircraft  whose  characteristics  are  included  in 
master  data  file  in  table  form  as  shown  in  Fig.  2.8.  The  user  inputs  the  percentages  of  the 
aircraft  which  comprise  the  aircraft  population,  expected  to  operate  at  the  runway  facility. 

Following  the  aircraft  mix  screen,  several  aircraft  characteristics  screens  for  each  aircraft 
selected  in  the  mix  screen  are  displayed  as  shown  in  Fig  2.9.  All  the  values  shown  in  this 
screen  are  transferred  from  the  master  data  file.  If  the  user  does  not  want  to  change  the 
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Analysis 


Press  escape  (Esc)  key  to  retire  to  rtftlK  HEW. 


Use  arrow  (t  or  t)  key  to  change  yoir  choice. 
Press  enter  (  *J  )  key  If  the  choice  Is  correct. 
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Figure  2.6  Analysis  Type  Menu. 
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values,  he/she  would  press  'Esc'  key  to  proceed  to  the  next  step  which  might  be  another 
aircraft  characteristics  screen  or  airport  environmental  data  screen. 

The  screens  for  'Airport  Environmental  Data"  and  the  "Airport  Operational  Data'  follow 
the  'Aircraft  Mix  and  Characteristics  Data'  screen.  These  screens  having  similar  table  formats 
are  shown  in  Figs  2.10  and  2.11.  The  next  screen  deals  with  "Runway  Gradients."  At  this  point, 
the  user  specifies  the  runway  length  with  the  line  input  method,  and  inputs  the  gradients  for 
every  one  tenth  section  of  the  runway  using  the  table  input  method  as  depicted  in  Fig.  2.12. 
The  final  screen  in  "Start  a  New  Problem'  mode  is  designed  for  "Weather  and  Exit  Speeds." 
On  this  screen,  the  relative  frequencies  of  occurence  for  dry  and  wet  conditions  are  specified. 
Also,  the  desired  exit  speeds  of  every  aircraft  category  under  each  weather  condition  are 
selected  by  the  table  input  method,  as  shown  in  Fig.  2.13. 


Figure  2.7  Data  Related  to  the  Type  of  Analysis. 
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2.4.2  Editing  Aircraft  Data 


This  portion  of  the  program  allows  the  user  to  modify  existing  data  files.  If  user  selects 
*2)  Edit  Data*  mode  at  the  main  menu,  the  edit  menu,  which  shows  six  groups  of  data,  is 
displayed.  In  the  "Edit  Data*  mode  the  user  can  select  the  group  of  data  which  he/she  wants 
to  change,  while  in  the  "Start  a  New  Problem*  mode  the  user  should  input  all  the  data 
sequentially.  The  details  for  editing  data  are  the  same  as  in  the  'Start  a  New  Problem"  mode. 
Fig.  2.14  shows  the  edit  menu. 

2.4.3  Editing  the  Master  File 


While  the  function  of  'Edit  Data"  mode  is  editing  the  working  data  file,  the  function  of  'Edit 
Master  File"  is  editing  the  master  data  file  which  keeps  the  aircraft  names  and  their  geometric 
characteristics.  If  'Edit  Master  File*  mode  is  selected,  the  edit  menu  for  master  data  file 
appears.  In  this  menu,  there  are  two  choices:  1)  "Add  a  New  Aircraft'  and  2)  "Change  some 
Specific  Data."  If  the  user  chooses  the  first,  he/she  has  to  select  one  out  of  five  aircraft 
categories  (TERPS  A-E)  and  input  the  new  aircraft  name.  Then  a  screen  for  editing  aircraft 
characteristics  appears.  If  the  user  opts  for  the  second  choice,  he/she  has  to  select  one 
aircraft  category  and  one  aircraft  name  included  in  the  category  selected.  Then  a  screen  for 
editing  aircraft  characteristics  appears.  Figs.  2.15  and  2.16  show  the  edit  menu  for  master 
data  file  and  the  screen  for  adding  a  new  aircraft,  respectively. 
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Input  Module 


Input  Module 


Figure  2.14  Edit  Menu  Screen. 


f  resi  escape  (Esc)  key  to  return  to  WIN  Will. 

Use  arrow  (t  or  I)  key  to  change  your  choice. 
Press  enter  (  *J  )  key  If  the  choice  Is  correct. 


2)Change  soar  Specific  lata 


Figure  2.15  Edit  Menu  Screen  for  Master  File. 


Input  Module 


3.0  Dynamic  Module 


The  Dynamic  Module  comprises  several  computational  subroutines  used  to  evaluate  the 
aircraft  performance  during  the  landing  flare,  the  runway  ground  roll,  and  the  turnoff 
maneuver.  For  modeling  purposes,  the  aircraft  flight  and  ground  paths  have  been  divided  into 
five  distinct  segments  as  shown  in  Fig.  1.10:  1)  an  air  distance,  2)  a  free  roll  distance  between 
the  air  and  braking  segments,  3)  the  braking  segment,  4)  a  second  free  roll  distance  between 
the  braking  and  turnoff  segments,  and  5)  the  turnoff  distance  until  clearing  the  runway  for  the 
next  arrival.  This  segmentation  eases  the  computational  tasks  in  the  model  and  also 
simplifies  the  inclusion  of  operational  policies  arising  at  specific  points  in  the  landing  phase. 
The  following  paragraphs  describe  in  detail  the  basic  assumptions  made  in  the  dynamic 
analysis  of  this  model.  It  should  be  understood  that  the  purpose  of  the  dynamic  model  is  the 
identification  of  potential  candidate  solutions  in  search  for  optimal  turnoff  locations  This  is 
accomplished  through  a  pre-screening  process  of  primary  and  secondary  candidates  as  'will 
be  explained  later  in  Chapter  5  of  this  report. 


3.1  Air  Distance 


The  air  distance  is  estimated  Indirectly  from  the  basic  aircraft  geometric  and  performance 
characteristics  contained  in  the  program  Master  File.  The  Master  File  contains  geometric  and 
basic  performance  characteristics  for  nearly  one  hundred  aircraft  in  current  use  These 
characteristics  are  used  to  estimate  the  approach  (VAP)  and  touchdown  speeds  (VTD)  for  each 
aircraft  selected  by  the  user.  Once  the  approach  speed  (or  reference  speed)  is  known,  an 
estimation  of  the  air  distance  can  be  made  assuming  a  circular  arc  flare  maneuver  flown  at 
constant  load  factor  to  transition  from  a  constant  rate  of  descent  flown  at  constant  descent 
flight  angle  y  on  final  approach  to  a  flat  flight  path  tangent  to  the  the  runway.  An  analytical 
expression  for  the  air  distance  can  be  found  by  equating  the  changes  in  kinetic  and  potential 
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DYNAMIC  MODULE  FLOW  DIAGRAM 


Exit  Path  Computations 

X  Coordinate  (XPAIII) 

Centripetal  Mu  Coeltaenl  (YMUC) 

Y -Coordinate  (YPATH) 

Inertial  Mu  Coollicienl  (YMUI) 

Instantaneous  Patti  Spoed  (VPATH) 

Scrubing  Mu  Coellicient  (YMUSC) 

Instantaneous  Yaw  Angle  (SAI) 

Total  Mu  Cooflicient  (YMUMAX) 

Time  in  the  Turn  (TTUBN) 

Instantaneous  Lateral  Acceleration  (NLAT) 

Instantaneous  Radius  ol  Turn  (RAPATH) 

Instantaneous  Rate  ol  Change  ol  RAPATH  (RDOT) 

Instantaneous  JerV  Acceleration  (JERK) 

implies  that  the  corresponding  run  oplion  has  been  selected  (i.e  ,  applicable  only  when  an  existing  90  deg 
angle  turnoll  has  been  selected) 


Figure  3.1  Dynamic  Module  Functional  Flow  Diagram. 


energy  of  the  aircraft  near  the  ground  with  the  product  of  a  retarding  force  F,  and  the  air 
distance  S,  as  shown  in  Eq.  3.1  [Nicolai,  1976;  Torenbeek,  1981;  Roskam,  1986] 


Sair  — 


2 

H thres _ Hare  V 

y  2  g  (nflare  -  1) 


(3.1} 


where,  Vn**  is  the  flare  speed  (taken  as  95%  of  the  approach  speed),  y  is  the  effective  descent 
flight  path,  the  threshold  crossing  altitude.  For  preliminary  analyses  the  flare  load  factor 
has  been  set  conservatively  to  1.15  g's  and  y  to  3  degrees  to  simulate  a  regular  ILS  approach 
flight  path.  Currently  the  dispersions  in  the  air  distance  are  set  internally  to  fixed  values  that 
depend  upon  the  aircraft  category  being  analyzed.  The  aircraft  categories  used  in  this 
research  are  consistent  with  those  implemented  in  the  FAA  Terminal  Operating  Procedures 
(TERPS)  and  defined  in  Table  3.1.  The  underlying  assumption  in  this  respect  is  that  slower 
aircraft  will  usually  experience  smaller  touchdown  dispersions  than  those  of  faster  aircraft  in 
absolute  distance  terms  (this  is  not  in  contradiction  to  the  fact  that  transport-type  pilots  might 
be  more  accurate  in  terms  of  touchdown  point  standard  deviations).  Actual  measurements 
of  lateral  and  longitudinal  landing  dispersions  for  transport-type  aircraft  made  by  Hosang 
[Hosang.  1975]  suggest  that  for  manual  control  landings  the  average  touchdown  dispersion 
(i.e.,  standard  deviation)  is  about  171  meters  (560  feet).  Although  little  data  is  available  in 
actual  instrument  meteorological  conditions  (IMC)  it  has  been  found  during  a  heads-up 
certification  display  that  reduced  touchdown  dispersions  prevail  under  this  circumstances 
[Desmond,  1986]. 


The  advantage  in  estimating  air  distances  relying  on  information  pertaining  to  each 
aircraft  is  two-fold;  1)  frees  the  analyst  from  relying  on  field  data  for  a  particular  aircraft  that 
in  most  cases  is  not  available  or  which  could  be  implemented  at  a  later  stage  for  calibration 
of  the  model,  2)  Introduces  more  realistic  variabilities  in  the  touchdown  locations  for  the  entire 
landing  aircraft  population  instead  of  assigning  a  fixed  touchdown  location  to  an  entire  aircraft 
category  population.  The  method  is  also  sensitive  to  specific  airfield  scenarios  since  more 
parameters  have  been  accounted  for.  For  example,  short  takeoff  and  landing  aircraft  can  be 
assigned  independently  different  values  for  the  flare  load  factor  and  descent  flight  angle  as 
they  occur  in  practice  thus  affecting  accordingly  the  air  distance  values  estimated  internally. 
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Table  3.1 

Aircraft  Approach  Category  Classification  (FAA,  1988). 

Category 

Landing  Speed  (1.3  Vst.n) 

A 

less  than  91  Knots 

B 

From  91  to  120  Knots  inclusive 

C 

From  121  to  140  Knots  inclusive 

D 

From  141  to  165  Knots  inclusive 

E 

166  Knots  and  higher 

The  time  consumed  in  the  air  phase  (7,/,)  Is  a  function  of  the  touchdown  location  (S,/f),  the 
approach  speed  ( V *),  and  the  touchdown  speed  (Vw).  Assuming  a  normal  distribution  for  the 
aircraft  touchdown  location,  T,„  and  its  corresponding  variance,  a},u,  are  given  as  follows: 


Tgrr  — 
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where,  V «,  and  Vtd  are  the  approach  and  touchdown  speeds,  respectively. 


3.2  Free  Roll  Distances 


Free  roll  distances  arise  in  the  aircraft  landing  operation  at  two  different  times:  1)  between 
the  air  distance  and  the  braking  stage,  7W1  and  2)  between  the  braking  segment  and  the 
turnoff  maneuver, 7W2.  The  firsl  free  roll  distance  tries  to  simulate  an  inherent  human  delay 
in  applying  aircraft  braking  mechanisms  such  as  thrust  reversers,  spoilers,  or  normal  wheel 

t 

braking.  A  conservative  average  value  of  three  seconds  has  been  allocated  for  this  transition 
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stage  with  a  typical  standard  deviation  of  one  second.  The  corresponding  free  roll  distance 
Sm,  and  its  variance  <r{mi  are  as  follows: 


Sfri  -  Tfr  1 


{3.4} 


2 


,/  2  2 
=  <*r, 


Ml 


{3.5} 


Note  that  any  reductions  in  aircraft  speed  during  wheel  spin-up  have  been  neglected  for  the 
sake  of  simplicity. 


The  second  transition  segment  tries  to  mimic  a  delay  time  arising  from  the  proper 
suppression  of  braking  action  and  a  recognition  time  of  the  turnoff  geometry  prior  to  exiting 
the  runway.  Under  all  visibility  conditions  there  is  a  delay  time  associated  with  the 
recognition  of  a  high-speed  turnoff  and  the  decision  of  whether  or  not  the  current  aircraft  state 
(i.e.,  speed,  braking  status,  etc.)  is  appropriate  to  negotiate  the  turn.  The  user  has  complete 
freedom  to  specify  this  delay  time,  TFR2,  in  the  Input  Menu  as  detailed  in  Chapter  3.  A  nominal 
value  of  2  seconds  has  been  allocated  for  this  parameter  as  a  representative  value  under 
clear  air  and  unlimited  visibility  conditions  (CAVU).  However,  the  analyst  could  increase  this 
value  accordingly  to  simulate  low  visibility  scenarios.  The  end  result  being  a  correspondingly 
higher  value  for  the  total  runway  occupancy  time  (ROT). 


A  complementing  assumption  in  this  analysis  is  that  free  roll  time,  TFRi  and  its  variance 
a\F„ j  are  known.  Then,  the  second  free  roll  distance  SfRl  and  its  variance  olFRt  are, 

SFR2  -  Vbrake„  TFR2  —  Vexlt  TFR2  (3  6) 


2 


—  WbrakaJ 


{3.7} 


where,  V*,.*.*  is  the  last  braking  speed  integrated  in  the  braking  stage  and  V is  the  aircraft 
exit  speed  selected  by  the  user. 
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3.3  Braking  Distance 


Under  normal  landing  conditions,  the  braking  segment  constitutes  the  largest  component 
of  the  Runway  Occupancy  Time  (ROT).  As  such,  it  becomes  necessary  to  estimate  with  some 
accuracy  the  braking  distance  if  one  is  to  have  some  confidence  in  the  total  distances  covered 
by  the  aircraft  on  the  ground.  The  previous  requirement  also  stems  from  the  incorporation 
of  airport  specific  variables  such  as  local  runway  slope  and  its  effect  on  aircraft  deceleration 
characteristics.  The  problem  seems  to  be  complicated  by  the  fact  that  many  aircraft 
parameters  necessary  to  determine  the  forces  and  moments  acting  on  the  aircraft  as  it  brakes 
are  not  only  time  dependent  (e.g.,  thrust  reverse  forces,  braking  forces,  parasitic  drag 
contributions,  etc.)  but  also  aircraft  specific  in  most  instances  (e.g.,  small  reciprocating  aircraft 
generally  do  not  have  thrust  reverse  capability  whereas  turbofan  and  large 
turbopropeller-driven  aircraft  do).  The  dilemma  is  then  to  use  a  model  that  will  provide  an 
accurate  answer  without  going  into  the  sophistication  and  computing  expense  of  a 
higher-order  model  (i.e.,  6-DOF  model). 

The  braking  algorithm  used  in  the  model  integrates  the  local  deceleration  of  the  aircraft, 
aKf,  as  it  travels  along  the  runway  (Eq.  3.8).  The  local  deceleration  is  estimated  from  the 
runway  Initial  conditions  specified  by  the  user  in  the  Input  portion  of  the  program.  At  the  same 
time  a  deceleration  multiplier,  cfrw,,  is  computed  throughout  the  integration  process  in  order 
to  correct  the  nominal  aircraft  deceleration  due  to  local  variations  of  runway  slope.  This 
simplistic  model  then  treats  the  aircraft  as  a  second  order  point  mass  model  whose  resultant 
deceleration  is  integrated  forward  in  time  to  obtain  the  velocity/distance  aircraft  state.  The 
assumption  of  a  constant  uncorrected  deceleration  rate  is  justifiable  if  one  realizes  that  in 
modern  aircraft  the  deceleration  rate  is  controlled  by  an  antiskid  system.  The  question  is  how 
can  we  estimate  the  deceleration  rate  for  either  each  aircraft  or  for  the  entire  aircraft 
population  from  the  known  runway  conditions?  As  the  reader  recalls  there  are  two  different 
scenarios,  dry  and  wet,  defined  in  the  runway  friction  characteristics  as  part  of  the  Input 
Module.  The  correlation  of  actual  aircraft  data  (Janes's,  1988;  Aviation  Week  &  and  Space 
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Technology,  1988;  Business  and  Commercial  Aviation,  1989]  is  done  backwards  to  estimate 
the  necessary  friction  coefficient  and  its  corresponding  deceleration  rate  necessary  to  match 
the  data  published  for  some  Known  conditions  such  as  those  corresponding  to  the  aircraft 
maximum  allowable  landing  mass  (MALW)  and  dry  pavement  conditions.  The  wet  condition 
braking  analysis  is  performed  with  the  introduction  of  a  deceleration  degradation  multiplier, 
c/rf,  into  the  baseline  deceleration  equation  (Eq.  3.9).  The  multiplier  has  been  derived  using 
NASA  and  ICAO  empirical  data  (Yager  and  White,  1981;  ICAO,  1966], 

A  second  correction  multiplier  is  also  introduced  in  this  analysis  to  modify  the 
instantaneous  deceleration  due  to  variations  in  the  local  runway  slope.  The  evaluation  of  this 
multiplier  has  been  done  outside  the  current  REDIM  Model  using  complete  drag/thrust  data 
for  a  Boeing  727-200  and  for  simplification  purposes  it  is  assumed  to  be  constant  for  all  the 
aircraft  population.  Equation  3.8  illustrates  the  approximation  made  of  the  braking  distance, 
Santo- 

n 

S brake  ~  ^  dt  V brake,  f3-8} 

I  =  1 

Furthermore,  decomposing  Vbr*tl  as  a  function  of  the  instantaneous  aircraft  deceleration  ( aKf ) 
and  the  deceleration  correction  factors  for  runway  friction  (cfrf)  and  runway  slope  ( cf, we 
obtain, 


n 

$  brake  =  (dt)  ( n )  (Vjcf)  +  (dt)  ^  (n  +  1  /)  cfrf.  cfrws .  agcf  {3.9} 

i=i 

where,  n  is  the  number  of  iterations  computed  in  the  simulation  of  the  braking  process  and 
whose  numerical  value  is  determined  by  the  integration  step  size,  dt,  and  Wo  is  the  touchdown 
speed.  The  computation  of  the  variance  of  Sbrtkt  denoted  by  is  estimated  as  follows, 


2 


{3.10} 
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where,  a1.^,  is  the  variance  of  the  deceleration  rate  (another  user  input)  and  dt  is  the 
simulation  step  size.  Note  that  the  time  consumed  in  the  braking  process  and  its  variance  are 
estimated  according  to  Eq.  3.11. 


^brake 
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{3.11} 


3.4  Turnoff  Algorithm 


The  turnoff  algorithm  integrates  the  aircraft  path  throughout  the  exit  maneuver.  The  exit 
maneuver  is  initiated  when  the  aircraft  reaches  the  user-defined  exit  speed  and  finalizes  with 
the  complete  clearance  of  the  runway  by  the  landing  aircraft  as  shown  in  Fig.  3.2.  In  order  to 
simplify  the  number  of  inputs  to  the  model  it  is  assumed  that  the  aircraft  wingtip  point  controls 
the  time  to  clear  the  runway.  This  is  generally  true  for  all  aircraft  exiting  at  high  speed. 
Exceptions  to  this  rule  are  small  aircraft  and  Short  Takeoff  and  Landing  Aircraft  STOL  (i.e., 
requiring  abnormally  large  tailplane  wingspans)  exiting  at  low  speed  (e  g.,  less  than  fifteen 
meters  per  second).  However,  since  the  objective  of  this  research  is  the  investigation  of 
high-speed  turnoffs  these  exceptions  would  seldom  occur  and  therefore  the  prediction  of  the 
clearing  point  can  be  done  adequately  with  a  single  aircraft  control  point. 

The  characteristic  motion  of  an  aircraft  turning  at  speeds  where  insignificant  aerodynamic 
control  can  be  exerted  by  conventional  primary  aerodynamic  surfaces  is  simplified  to  the 
forces  acting  on  the  nose  landing  gear.  An  algorithm  developed  by  Schoen  et  al  [Schoen  et 
al,  1985]  and  used  in  a  previous  NASA  research  effort  on  this  topic  considers  three  side  force 
contributions  acting  on  the  aircraft  nose  landing:  1)  the  centripetal  force,  2)  the  aircraft  inertia, 
and  3)  the  tire  scrubbing  resistance  to  the  turn.  Mathematically,  the  nondimensional 
contributions  to  the  nose  gear  are, 

+  Ac  +  A*C  {3-12} 
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where,  jifkW  is  tho  nose  gear  tire  skidding  friction  coefficient,  /i<„  is  the  aircraft  inertia 
contribution  term  to  ‘.he  nose  gear  side  load.  /ie  is  the  centripetal  acceleration  contribution, 
and  n,c  is  the  tire  scrubing  resistance.  These  contributions  are  calculated  In  Eqs.  3  13  to  3.15. 


where,  /„  is  aircraft  moment  of  inertia  about  the  vertical  axis,  in  Kg-m-m,  a  is  the  angular 
acceleration  (rad/sec,)  of  the  aircraft  fuselage  as  it  executes  the  turning  maneuver,  wb  is  the 
aircraft  wheelbase  (meters),  Im  is  the  aircraft  mass  supported  by  the  main  gear  (in  percent), 
g  is  the  gravitational  constant  (m/sec-sec),  m  Is  the  total  aircraft  mass  (Kg.).  V  is  the 
instantaneous  speed  (m./sec.)  of  the  nose  gear,  and  R  is  the  instantaneous  radius  of  the  curve 
(m.).  Further  breakdown  of  the  angular  acceleration  yields  for  Eq.  3.13  the  following. 


(V  R-V  R 


**'**  m  g  wb  Iml  100  (1  -  /m/100) 


{3.16} 


where,  R  represents  the  rate  of  change  of  the  turning  radius  of  curvature,  V  if  the 
Instantaneous  velocity  rate  of  change  of  the  nose  gear,  and  R  and  V  are  the  state  variables 
of  our  system.  A  further  simplification  can  be  introduced  if  the  term  V  R  is  neglected  on  the 
grounds  of  very  small  values  for  the  deceleration  rate  through  the  turnoff  maneuver.  Fact  that 
has  been  found  true  In  the  empirical  studies  of  Horonjeff  and  Hosang  [Horonjeff  et.  al.,  1958; 
Hosang,  1975].  The  new  expression  becomes, 
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fa[-VW] 

>a  R2  m  g  wb  /m/100  (1  -  /m/100) 


{3.17} 


solving  for  the  rate  variable,  R  and  integrating  over  time  it  is  possible  to  estimate  the  state 
variables  of  the  motion, 

n,  R2 

R  =  m  g  wb  /m/100  (1  -  /m/100)  {3.18} 
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where,  X  and  Y  are  the  position  coordinates  of  the  vehicle  as  it  progresses  into  the  turn  and 
^  is  the  heading  angle  that  the  nose  gear  makes  with  a  global  axis  system  centered  about 
initial  position  of  the  turnoff  path  (Fig.  3.2).  Eqn.  3.18  is  further  restricted  by  passenger  comfort 
factor  limitations  as  will  be  explained  in  Section  3.5  of  this  Chapter. 


The  neglect  of  aerodynamic  effects  in  this  analysis  might  be  realistic  for  up  to  about  two 
thirds  of  the  landing  speed  (V/*,)  as  this  is  known  to  be  the  threshold  for  significant 
aerodynamic  control  for  conventional  aircraft  [Miller,  1967).  Even  with  this  restriction,  the 
evaluation  of  runway  turnoffs  can  be  accomplished  for  a  large  range  of  aircraft  speed  values 
ranging  from  10  to  45  m/sec.  (22.3-100.4  MPH)  for  transport-type  category  aircraft.  As  will 
become  evident  during  the  discussion  of  results  of  this  research  turnoff  designs  above  45 
m/sec.  (100.4  MPH)  are  probably  unlikely  due  to  lateral  space  limitations  following  the  turn. 
This  topic  is  currently  being  investigated  in  a  continuous  research  carried  by  the  Center  for 
Transportation  Research  with  the  FAA  and  NASA. 
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Lateral  Range  (m.)  o  Lift-to-Weight  Ratio 


Speed  (m/sec) 

ire  3.3  Comparative  Lift-to-Welght  Ratios  for  High-Speed  Taxiing  Aircraft  (Low  «). 


■  Jerk  Limitation 
•  No  Lift 

■  Lift  Included 


Downrange  (meters) 

Figure  3.4  Aircraft  Turnoff  Tracks  with  and  without  Lifting  Forces  (Boeing  727-200). 
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The  inclusion  of  the  lifting  forces  acting  on  the  aircraft  at  high  speed  can  be  added  to  Eqn.  3.18 
replacing  the  mass  term  by  an  equivalent  force  that  accounts  for  the  potentially  large  lifting 
forces  experienced  at  high  speeds. 


R  = 


r2 
Izz  V 


(m-.SpV2S  CL)  g  wb  /m/100  (1  -  /m/100) 


{3.22} 


where,  p  is  the  air  density,  CL  is  the  aircraft  lift  coefficient  in  ground  effect  and  the  landing  flap 
configuration,  S  is  the  wing  area  and  V  is  the  aircraft  speed.  Fig.  3.3  depicts  the  lift-to-weight 
ratio  for  three  types  of  aircraft.  The  significance  of  this  plot  is  that  it  allows  us  to  establish 
desired  exit  speed  boundaries  for  each  TERP  aircraft  category.  It  can  be  seen  from  this  plot 
that  for  a  typical  single-engine  aircraft  the  L/W  ratio  below  15  m/sec.  can  be  practically 
ignored.  Similarly,  for  medium  sized  and  heavy  transport  aircraft  the  lift  effects  become 
noticeable  after  25  m/sec.  (56  MPH)  but  for  all  practical  purposes  the  turnoff  trajectories  are 


Coefficient  of  Skidding  Friction 


Figure  3.5  Skid  Friction  Variations  with  Tire  Pressure  and  Speed. 
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not  significantly  altered  but  after  45  m./sec.  or  more  as  depicted  in  Fig.  3.4  The  reader  is 
warned  that  these  results  were  derived  under  the  assumption  that  the  lift  contribution  acts 
near  the  aircraft  center  of  gravity  and  therefore  no  significant  pitching  moment  is  induced  to 
load  the  nose  gear.  In  practice  pilots  can  modify  the  nose  gear  load  distribution  by  deflection 
of  the  aircraft  elevator  surfaces  once  these  become  effective  (about  two  thirds  of  the  liftoff 
speed)  thus  making  more  difficult  the  task  of  arriving  at  a  unique  conclusion. 

Equally  important  is  the  fact  that  Fig.  3.3  depicts  the  lift  generated  by  an  aircraft  at  small 
angles  of  attack  and  flaps  down  (i.e.,  rolling  on  a  highspeed  turnoff).  However,  the  potential 
lifting  force  capable  of  being  generated,  is  much  larger  if  the  angle  of  attack  is  increased 

r  i 

through  the  use  of  elevator  power.  This  is  a  fact  of  fundamental  importance  if  one  is  to  restrict 
aircraft  from  turning  at  excessive  speeds  where  the  aircraft  might  actually  be  flying!  Taking 
as  example  a  light  single-engine  aircraft  such  as  the  Cessna  172  it  is  seen  from  Fig.  3.3  that 
at  25  m./sec.  the  lift-to-weight  ratio  is  about  .42  for  low  angles  of  attack  (i.e.,  2-3  degrees)  in 
the  landing  configuration  (i.e.,  flaps  fully  down).  However,  this  happens  to  be  the  stalling 
speed  for  this  aircraft  in  the  landing  configuration  (for  an  angle  of  attack  of  about  16  degrees) 
and  therefore  it  would  be  unreasonable  to  turn  this  aircraft  at  such  high  speeds  REDIM 
currently  has  a  high  speed  threshold  limiter  to  overcome  this  complication  by  announcing  ‘he 
potential  input  error  and  suggesting  upper  and  lower  bounds  for  the  aircraft  exit  speeds  A 
more  sophisticated  model  including  the  aerodynamic  terms  in  the  aircraft  equations  of  motion 
would  necessitate  the  knowledge  of  several  important  aerodynamic  derivatives,  geometric 
parameters,  and  reference  areas  for  every  aircraft  considered  imposing  a  large  computational 
burden  for  the  optimization  procedure  and  necessitating  a  more  complete  input  data  set  from 
the  user.  This  approach  was  then  ruled  out  due  to  the  complexity  of  the  input  data  needed. 

The  aforementioned  algorithm  has  been  modified  in  order  to  account  for  the  large 
variations  in  skidding  friction  coefficients  observed  for  a  large  aircraft  population.  It  is  well 
documented  in  the  literature  that  the  skidding  friction  coefficient  is  a  function  of  aircraft  tire 
pressure  and  aircraft  speed,  among  other  variables  [NASA  TN  4418,  1966;  Wong,  1978],  A 
summary  of  this  functional  relationship  ir.  'lopjcted  graphically  in  Fig  3.5  where  four  aircraft 
tire  pressures  are  represented  in  this  figure  and  they  correlate  well  with  the  four  different 
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aircraft  categories  modeled  in  this  research.  The  upper  curve  corresponds  to  a  tiro  pressure 
of  50  PSI  (pounds  per  square  inch)  and  is  representative  of  the  characteristics  of  TERP  A 
category  aircraft.  Similarly,  tho  lower  curve  represents  a  200  PSI  tire  pressure  typical  of 
current  heavy  aircraft  (i.e.,  Boeing  747,  DC-10,  L-1011,  etc  ). 

Complementing  this  algorithm  a  small  forward  deceleration  can  be  introduced  in  order  to 
account  for  the  sma'I  speed  losses  expected  while  turning.  The  rolling  friction  opposing  the 
motion  of  the  aircraft  on  the  ground  introduces  a  deceleration  rate  proportional  to  the  product 
of  g  and  fnn  where  this  last  term  is  the  coefficient  of  rolling  friction.  For  the  sake  of  simplicity 
fni,  is  taken  constant  with  speed  although  it  is  known  to  vary  with  tire  speed  as  well.  A  typical 
value  of  .03  is  used  for  fnU  for  the  base  model. 

a roll- 0  I roil  {3.23} 

An  Euler  first-order  integrating  scheme  is  used  to  solve  numerically  the  aircraft  equations 
of  motion  through  the  turnoff  maneuver.  The  time  spent  on  the  turn,  Tp,th  is  considered  to  be 
deterministic  in  nature.  A  baseline  step  size  of  one  hundredth  of  second  was  found  to  offer 
accurate  results  within  the  desired  computational  time  limitations  for  the  program.  The 
accuracy  of  the  method  is  evident  from  Fig.  3.6  where  the  first-order  solution  is  compared  with 
an  equivalent  Runge-Kutta  fourth-order  scheme  (notice  that  both  curves  overlap)  These 
solutions  were  obtained  using  the  geometric  and  performance  parameters  of  a  Boeing  727-200 
and  as  can  be  seen  from  the  turnoff  paths  generated  the  results  are  within  one  half  of  a 
percent  of  each  other  (i.e.,  less  than  half  a  meter  difference  between  both  solutions  at  the  end 
of  a  high-speed  turnoif).  Another  justification  for  the  Euler  algorithm  was  the  desired  accuracy 
in  stopping  the  simulation  as  closely  as  possible  to  the  runway  clearance  point  (Fig.  3.2).  With 
the  current  step  size  it  is  possible  to  ascertain  the  turnoff  time  (TOT)  and  the  lateral  range 
distance  within  very  small  windows,  .01  seconds  or  .15  meters,  for  an  aircraft  traveling  at  30 
m/sec.  (67.2  MPH)  and  reaching  the  runway  clearance  point  with  up  to  30  degrees  of  total 
heading  change. 
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Boeing  727-200 
35  nV»«c.  Exit  Speed 
Wet  Runway  Condition 


Euler 

RK-4  Method 


Downrange  (meters) 


Figure  3.6  Comparison  of  Track  Simulation  Results  Using  two  Integrating  Methods. 


The  aircraft  position  coordinates  in  the  turn  (Xp.w,  and  (Yp,«,),  the  aircraft  speed  ( Vp,th ),  and 
the  aircraft  instantaneous  heading  {^)  constitute  state  variables  through  the  turnoff  maneuver 
(Eqns.  3.20-22)  These  states  are  integrated  forward  in  time  to  assess  the  instantaneous 
turning  radius  (Rp,o,)  and  ultimately  estimate  the  position  changes  experienced  by  an  aircraft 
as  a  high-speed  turnoff  is  negotiated.  The  aircraft  is  considered  to  have  cleared  the  runway 
when  its  right  wingtip  has  traveled  the  lateral  distance  necessary  to  cross  the  runway  edge 
imaginary  plane.  Once  the  turnoff  path  and  times  are  estimated  it  is  possible  to  ascertain  the 
time  from  threshold  crossing  to  the  end  of  the  turnoff  maneuver.  Since  some  of  the  distances 
and  times  involved  in  the  process  are  random  variables  the  net  effect  is  that  runway 
occupancy  time  (ROT)  and  the  total  distance  to  initiate  the  turnoff  are  both  probabilistic  in 
nature. 


Stum  ~  $air  +  $FR  1  +  $  brake  +  SFR2 


{3.24} 


Dynamic  Module 


52 


{3.25} 


2  _  2  ,  2  ,2  ,2 

^tum  -  as*  +  +  *S6fM.  +  °SfRt 

hum  —  Tatr  4-  TFRA  4-  Tbrake  4-  TFR 2 

2  2  2  2  2 
=  +  "T",  +ffw#  +  ffrW2 

WOT  =  rtorr  4-  {3.28} 

4or*«L.  {3.29} 

where,  Smm  is  the  distance  from  the  threshold  to  the  initiation  of  the  turnoff  (i.e.,  exit  location 
distance),  a5,„m  is  the  variance  of  this  previous  parameter,  TWm  and  olwrn  are  the  time 
consumed  from  threshold  to  the  initiation  of  the  turnoff  and  its  corresponding  variance.  TcM. 
is  the  time  in  the  turnoff  and  ROT  is  the  total  runway  occupancy  time  for  a  single  aircraft  with 
variance  ohor. 


{3.20} 
13  27} 


3.5  Comfort  Factor  Considerations 


At  this  point  it  is  important  to  introduce  and  discuss  the  vehicle  limitations  due  to  the 
passenger  comfort  factor.  The  measures  of  effectiveness  used  to  estimate  passenger  comfort 
in  a  turning  vehicle  have  traditionally  been  the  normal  acceleration,  a™  and  the  two  vectorial 
components  of  the  jerk  usually  defined  in  the  literature  as  the  normal  and  tangential  jerk,  Jn 
and  Ju  respectively.  There  seems  to  be  little  information  in  the  literature  regarding  the  human 
comfort  'thresholds'  to  lateral  accelerations  and  jerks.  Most  of  the  data  seems  to  have 
concentrated  around  motions  in  the  plane  of  symmetry  of  transportation  vehicles  (i.e.,  pitch 
rate  and  vertical  accelerations).  Data  from  the  railroad  industry  seems  to  offer  the  only 
tractable  guidelines  for  both  lateral  acceleration  and  jerk  (Hulbert,  1979  and  Wright,  1989). 
Average  accepted  limitations  for  normal  jerk  oscillates  between  0.055  to  0.065  g's  (i.e.,  0.54-64 
m/sec-sec.)  whereas  that  for  lateral  acceleration  is  about  0.12-0.15  g's  (i.e.,  1.18-1.47 
m/sec-sec).  It  is  however  important  to  understand  that  these  represent  train  threshold  values 
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which  might  be  overly  conservative  when  applied  to  aircraft  passengers  Aircrall  pilots  can 
adjust  laterally  the  trajectory  of  the  vehicle  on  the  ground  and  consequently  have  a  large 
influence  in  the  lateral  acceleration  and  jerk  perceived  by  a  passenger  Also,  the  suspension 
system  found  In  typical  aircraft  is  better  suited  to  absorb  rolling  and  yawing  motions  than 
those  found  in  trains  and  this  should  favor  the  ground  riding  qualities  of  air  vehicles.  Finally, 
aircraft  seats  are  usually  designed  with  more  lateral  and  longitudinal  restraints  than  those 
found  in  trains  thus  hiding  the  passenger  perception  of  uncomfortable  motions.  These  facts 
need  further  research  to  assess  their  validity  and  should  be  interpreted  as  tentative. 


Schoen  et.  al.  [Schoen  et.  al.,  1985]  used  a  combined  performance  measure  to  relate  the 
upper  limits  of  the  normal  acceleration  and  the  jerk.  The  suggested  relationship  is  given  in 
Eq.  3.30  and  shown  graphically  in  Fig.  3.7.  It  must  be  pointed  oui  that  this  comfort  factor 
modeling  has  been  suggested  in  the  literature  but  has  not  been  correlated  with  experimental 
data  validating  the  results.  It  seems  possible  that  this  method  might  be  too  restrictive  when 
applied  to  aircraft  ground  motion  due  to  the  larger  lateral  restraint  mechanisms  offered  by 
aircraft  seats  when  compared  with  their  train  counterparts.  Also  the  reader  should  recall  that 
normal  acceleration  and  normal  jerk  are  related  according  to  the  functional  form  shown  in 
Eq n.  3.31. 


^max 
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4, 

‘Anax 


<  1.00 


{3.30} 


an 


Jn  I 

ve 


{3.31} 


where,  a„  is  the  normal  acceleration,  V,  is  the  entry  exit  speed  (i.e.,  assumed  constant 
throughout  the  turnoff  for  the  transition  curve  practical  analysis),  am„  is  the  maximum 
permissible  normal  acceleration,  J„  is  the  normal  jerk  and  Jm„  is  the  maximum  tolerated 
normal  jerk  value. 


The  kinematic  equations  of  motion  defining  a  turning  vehicle  through  a  transition  spiral 
are  shown  in  Eqns,  3.32  and  3.33  (in  cartesian  coordinates).  The  approximation  represents  a 
truncation  of  the  Taylor  series  expansion  up  to  the  third  term.  Note  that  an  iterative  solution 
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Figure  3.7  Possible  Passenger  Comfort  Indifference  Curves. 


to  find  x  and  y  can  easily  be  implemented  if  a  first  estimate  is  made  with  the  first  order  term 
of  Eqns  3.32  and  3.33  without  a  great  loss  in  accuracy  [Anderson,  1979], 
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{3.32} 


{3.33} 


where,  x  and  y  are  the  coordinates  of  the  turnoff,  J„  is  the  normal  jerk,  /  is  the  curve  length 
and  V,  is  the  aircraft  speed  at  the  entrance  of  the  turnoff. 


Figure  3.8  compares  two  high-speed  turnoff  tracks  for  two  transport-type  aircraft  and  two 
different  speeds  with  their  corresponding  jerk-limited  turnoff  tracks.  The  threshold  jerk  value 
used  for  this  figure  was  0.55  ml  sec5  which  has  been  accepted  by  many  researchers  in  the 
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Lateral  Range  (m.)  Lateral  Range  (m.) 


50%  Skid  Friction  Safety  Factor 
125 -t  30  mJtte.  Exit  Speed 


a  B- 727- 200 
♦  B- 747-400 
■  Jerk- Limited 


a  B-747  400 
•  B  727-200 
a  Jerk-Limited 


Downrange  (m.) 

Figure  3.8  Comparison  of  High-Speed  Turnoff  Tracks  for  Various  Aircraft. 
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train  industry,  it  should  be  pointed  out  that  both  paths  are  very  close  to  each  other  in  the  first 
one  hundred  meters  into  the  turnoff  maneuver.  The  jerk-limited  track  quickly  ceases  to  be  of 
concern  once  this  point  is  reached.  A  magnification  of  these  results  of  these  results  during  the 
first  one  hundred  meters  reveals  that  the  jerk-limited  transition  curve  slightly  dominates  the 
turnoff  geometry  and  thus  an  enlargement  of  the  local  radius  of  curvature  might  be  desirable. 
It  should  be  pointed  out  .  however,  that  the  differences  in  both  curves  are  so  small  that  the 
resulting  practical  geometry  remains  almost  unaffected.  The  conclusion  of  this  discussion 
seems  to  suggest  that  further  investigation  needs  to  be  carried  out  in  actual  vehicles  in  order 
to  determine  if  these  differences  are  really  important  from  a  practical  point  of  view.  Also,  it  is 
suggested  that  more  research  needs  to  be  carried  out  in  the  area  of  aircraft  ground  ride 
qualities  and  human  comfort  factors  to  lateral  accelerations.  The  data  and  models  reviewed 
in  our  research  are  inconclusive  in  terms  of  well  defined  thresholds  of  comfort.  It  is  possible 
that  the  normal  acceleration  and  jerk  tolerances  accepted  by  train  passengers  could  in  fact 
be  overly  conservative  when  applied  to  aircraft  ground  scenarios.  It  is  time  to  fill  this  void  in 
aerospace  research. 


The  equations  of  motion  of  the  point  mass  vehicle  can  then  be  effectively  modified  to 
account  for  the  comfort  limitations.  The  main  effect  of  restricting  the  turnoff  geometry  to 
comply  with  a  minimum  jerk-limitation  threshold  is  to  restrict  the  rate  of  change  in  the  turnoff 
curvature.  For  a  constant  speed  transition  spiral  the  first  order  differential  equation  defining 
the  rate  of  change  of  the  heading  angle  with  respect  to  time  is, 


dJ/_  _  dip  dl 
dt  ~  dl  dt 


d<P  djp_  _v 

dt  dl  R 


{3.34} 


{3.35} 


where,  \p  is  the  instantaneous  heading  angle,  /  is  turnoff  curve  length,  V  is  the  aircraft  turnoff 
speed,  R  is  the  instantaneous  turning  radius  and  f  is  the  time  in  the  turnoff  maneuver. 
Differentiating  Eqn.  3.34  with  respect  to  time  and  knowing  the  value  of  the  second  derivative 
of  ip  with  respect  to  time  and  neglecting  (or  a  moment  the  deceleration  term  the  limiting 
expression  of  R  limited  by  jerk  constraints  is  shown  in  Eq.  3  36 
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Equation  3.37  can  then  be  used  to  restrict  the  value  of  R  for  any  limiting  values  of  Jn  and  a„ 
(i.e.,  for  values  previously  defined  as  amtx  and  J^.  It  is  interesting  to  note  that  although 
previous  research  has  tried  to  justify  a  combined  passenger  comfort  performance  index  in 
terms  of  these  two  parameters  they  arise  naturally  in  the  geometric  definition  of  a  spiral  curve 
and  are  in  fact  related  as  seen  in  Eqn.  3.37.  Consequently,  a  linear  combination  of  the  two 
parameters  is  very  unlikely  to  dictate  passenger  a  comfort  indifference  curve.  In  these 
author's  opinion  the  indifference  curve  could  probably  resemble  a  quarter  of  an  ellipse  rather 
than  being  linear  although  further  research  will  ultimately  dictate  these  thresholds.  Careful 
examination  of  Eqns.  3.34  through  3.37  indicates  that  jerk  limitations  dominate  the  initial 
portion  of  a  constant  speed  transition  spiral  whereas  the  upper  value  of  normal  acceleration 
takes  precedence  as  the  turnoff  geometry  progresses  in  time  as  the  result  of  a  decreasing 
turnoff  curvature  over  time. 


3.6  Exit  Assignment  Algorithm 


The  exit  assignment  algorithm,  as  its  name  implies  estimates  ’he  probable  exit  that  an 
aircraft  would  take  under  a  given  set  of  operating  conditions.  The  algorithm  is  used  in  all  of 
REDIM  run  options  and  assigns  either  existing  or  potential  turnoffs  to  every  aircraft  according 
to  their  landing  performance  characteristics.  In  the  design  and  redesign  running  modes  the 
algorithm  assigns  all  exits  downrange  of  the  primary  candidate  generated  by  the  /  th  aircraft 
in  question.  This  is  necessary  as  every  other  exit  downrange  constitutes  a  potential  optimal 
solution  for  the  optimization  routine.  This  will  become  more  evident  in  Section  4  3  2  of  this 
report  where  the  dynamic  programming  formulalion  is  explained  in  detail. 


Dynamic  Module 


58 


3.7  Airport  Environmental  Variables 


It  was  said  in  Chapter  2  that  the  airport  environmental  variables  were  defined  in  a  single 
screen  bearing  a  similar  name.  The  environmental  characteristics  of  interest  are:  1)  wind 
speed  (WSPEED),  2)  wind  direction  (WDIR),  3)  airport  elevation  (AIRELV),  4)  airport 
temperature  (AIRTEMP),  5)  runway  orientation  (RUNOR),  6)  runway  visual  range  (RVR),  7) 
runway  width  (RUNWID),  and  8)  distance  to  nearest  taxiway  (DISTT). 

The  wind  vector  is  used  in  conjunction  with  the  runway  orientation  to  estimate  the 
longitudinal  and  lateral  wind  components  affecting  aircraft  operations.  The  longitudinal  wind 
component  affects  the  landing  speeds  of  the  aircraft  population  and  as  such  has  a  direct 
impact  in  the  runway  occupancy  time  and  turnoff  locations.  Regarding  the  use  of  a  single  wind 
vector  as  input  to  the  model,  the  the  user  is  urged  to  execute  the  baseline  program  under  the 
average  prevailing  wind  conditions  at  the  airport  facility  just  as  he/she  would  do  under  the 
average  prevailing  temperature. 

Temperature  and  airfield  elevation  have  a  direct  impact  in  the  performance  of  the  aircraft 
in  the  air  and  on  the  ground.  Changes  to  the  aircraft  equivalent  airspeed  (EAS)  due  to 
temperature  and  fie.  ‘  elevation  can  have  large  impact  in  the  ROT  and  the  turnoff  location 
parameters  as  will  be  seen  in  Chapter  6  of  this  report.  The  model  converts  equivalent  speeds 
(EAS)  to  true  air  speeds  (TAS)  to  estimate  the  stalling  {Vst,n)  and  approach  speeds  and 
ultimately  predict  the  aircraft  landing  roll  performance.  The  runway  width  and  runway 
distance  to  nearest  taxiway  are  included  in  this  set  of  parameters  in  order  to  estimate  the  time 
sp#*nt  on  the  turnoff  maneuver  by  each  aircraft. 
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3.8  Aircraft  Characteristics 


The  aircraft  characteristics  used  in  the  model  are  shown  in  Table  3.2.  These  are 
necessary  to  estimate  the  aircraft  performance  on  the  ground  as  well  as  in  the  flare 
maneuver.  The  aircraft  mass,  wing  area,  and  the  maximum  landing  lift  coefficient  dictate  the 
approach  speed  and  hence  affect  the  ROT  and  exit  location.  It  is  also  used  to  estimate  the 
second  moment  of  inertia  of  the  aircraft  around  the  vertical  axis  (/**)  ultimately  influencing  the 
turning  aircraft  capabilities  through  an  exit.  Roskam  [Roskam,  1985)  suggests  a  logarithmic 
relationship  between  these  two  parameters  which  seems  to  correlate  very  well  for  all  aircraft 
TERP  categories.  The  regression  equation  in  metric  units  is  shown  in  Eq  3  29  where  the 
aircraft  mass  is  given  in  kilograms  and  the  moment  of  inertia  in  kg-m-m. 

The  aircraft  wheeltrack  (ACFWT)  is  used  to  estimate  the  maximum  track-in  distance 
present  during  the  turnoff  maneuver.  The  track-in  distance  is  defined  as  the  perpendicular 
distance  measured  from  the  geometric  center  of  the  aircraft  main  gears  to  the  imaginary  path 
followed  by  the  nose  gear.  Track-in  distances  are  used  to  assure  a  sound  geometric  design 
of  the  high-speed  turnoff.  It  should  be  pointed  out  that  in  general  track-in  distances  tend  to  be 
relatively  small  for  very  high  speeds  (i.e.,  >  30  m/sec.)  However,  for  large  aircraft  and 
medium  speeds  they  should  be  considered  in  the  geometric  design.  Fig.  3.9  shows  graphically 
the  nomenclature  used  to  model  the  aircraft  kinematic  behavior  including  the  estimation  of  the 
track-in  distance.  As  the  dynamic  simulation  executes  a  sample  record  of  the  main  gear 
position  {Xmi,Ymi )  is  kept  and  the  track-in  distance  is  evaluated.  A  simple  sorting  routine 
searches  for  the  largest  value  of  track-in  and  this  is  later  transfered  to  the  output  module  to 
calculate  the  corresponding  turnoff  geometry  that  satisfies  the  kinematic  constraints  of  the 
turnoff  track.  As  usual,  a  safety  distance  is  selectively  used  to  estimate  the  distance  from  the 
centerline  of  the  the  turnoff  track  to  the  edge  of  the  pavement.  No  judgemental  oversteering 
is  assumed  in  the  program  as  this  is  certainly  not  recommended  for  an  aircraft  traveling  at 
high-speed  on  the  ground. 
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Aircraft  Position  at  Time  : 


Aircraft  Position  at  Time  t 


(Xcg,  Ycg) 
(Xm,  Ym) 
(Xml,  Yml) 
(Xmr,  Ymr) 
Castor 
V 

T  rack-in 


Najngnc.lgt.urs 

Aircraft  Center  or  Gravity  Coordinates 
Coordinates  of  Main  Landing  Gear  Geometric  Center 
Left  Main  Landing  Gear  Position  Coordinates 
Right  Main  Landing  Gear  Position  Coordinates 
Castor  Angle 

Aircraft  Velocity  Vector  (at  Center  of  Gravity) 

Track-In  Distance  (from  Nose  Gear  to  Main  Gear 
Geometric  Center) 


Figure  3.9  Aircraft  Kinematic  Behavior  Nomenclature. 


•At 
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ACFWS  Aircraft  Wing  Span 

ACFWB  Aircraft  Wheel  Base 

NWTPLANE  Distance  from  Nose  Dear  to  Wing  Tip  Plane 

ACFWT  Aircraft  Wheel  Track 


Figure  3.10  Aircraft  Geometric  Characteristics. 
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Another  geometric  parameter  Included  here  is  the  distance  from  the  aircraft  nose  gear  to 
the  imaginary  plane  passing  through  the  airplane  wingtips.  This  distance  is  used  as  the 
controlling  point  to  ascertain  whether  or  not  the  aircraft  has  cleared  the  runway.  A  graphical 
description  of  some  of  these  parameters  is  seen  in  Fig.  3.10. 

/„  =  Antilogy  {1.7215  log10  (m)  -  1.6730}  {3.29} 


Table  3.2  Aircraft  Characteristics. 

Name 

Variable 

Remarks 

Aircraft  Mass 

ACFMASS 

Max.  Lending  Mass  (Kg.) 

Aircraft  Wheelbase 

ACFWB 

in  meters 

Aircraft  Wheeltrack 

ACFWT 

in  meters 

Aircraft  Load 
on  Main  Landing  Gear 

ACFLM 

At  aft  C.G.  (%) 

Aircraft  Wing  Area 

ACFWA 

Gross  wing  area  (sq.  m.) 

Aircraft  Maximum 

Lift  Coefficient 

ACFCL 

At  max.  flap  setting  (dim  ) 

Distance  from  Nose 
Gear  to  Wingtip 

NWTIP. 

in  meters 
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4.0  Optimization  Module 


The  dynamic  simulations  of  aircraft  landing  movements  calculate  the  best  turnoff 
locations  for  each  aircraft  in  both  dry  and  wet  runway  surface  conditions.  The  best  turnoff 
location  is  defined  as  the  nearest  location  from  the  runway  threshold  where  the  aircraft 
decelerates  to  the  pre-specified  desirable  exit  speed  with  the  pre-specified  reliability.  If  the 
aircraft  reduces  its  speed  to  the  pre-specified  exit  speed  before  reaching  the  assigned  turnoff 
location,  the  aircraft  will  be  considered  to  exit  the  runway  successfully.  Reliability  is  defined 
as  the  probability  that  the  aircraft  exits  the  runway  successfully.  For  example,  if  the  reliability 
is  specified  as  90%,  90  aircraft  out  of  100  landing  attempts  will  exit  the  runway  successfully. 

If  an  exit  is  constructed  at  the  best  turnoff  location  for  an  aircraft,  the  runway  occupancy 
time  (ROT)  of  the  aircraft  will  be  minimized  without  sacrificing  the  reliability  Though  some 
exits  constructed  ahead  of  the  best  turnoff  location  can  produce  less  ROT,  it  is  not  permissible 
to  assign  the  aircraft  to  these  exits,  since  reliability  must  be  sacrificed. 

Suppose  there  are  five  aircraft  in  consideration,  the  simulations  of  the  aircraft  landing 
movements  will  provide  ten  different  turnoff  locations  for  each  aircraft  and  two  runway  surface 
conditions.  The  goal  of  an  optimization  algorithm  is  to  find  a  few  locations  (e.g.  2  or  3)  at 
which  all  the  aircraft  in  consideration  can  exit  the  runway  with  the  minimum  weighted  sum  of 
ROT.  Since  each  aircraft  and  each  surface  condition  can  have  different  relative  frequency,  the 
weighted  sum  of  ROT  should  be  minimized  instead  of  total  ROT.  Figure  4.1  illustrates  the  best 
turnoff  locations  and  their  relative  frequencies.  Let  /,,  and  w*  represent  the  best  turnoff  location 
and  relative  frequency  for  aircraft  i  and  surface  condition  j. 

The  optimization  procedure  in  REDIM  conducted  with  the  following  steps: 

1.  Generate  the  complete  set  of  candidate  locations. 

2.  Calculate  the  ROT  of  every  aircraft  for  each  candidate  locations. 

3.  Find  the  optimal  location(s)  out  of  the  'mndidates. 

4.  Assign  aircraft  to  the  optimal  locations 
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Figure  4.1  The  Best  Turnoff  Locations  and  Their  Weights. 


Steps  1  and  2  are  the  data  preparation  for  step  3  which  is  the  mathematical  optimization.  The 
optimization  in  step  3  employs  a  dynamic  programming  technique.  Step  4  is  the  interpretation 
of  the  optimization  results  into  a  practical  solution.  A  flow  chart  of  the  optimization  is  depicted 
in  Figure  4.2.  The  notations  used  in  the  flow  chart  are  explained  in  the  following  sections. 

4.1  Generation  of  a  Complete  Set  of  Candidates 


Finding  optimal  turnoff  locations  is  a  continuous  optimization  problem.  That  is,  an  optimal 
turnoff  location  car.  be  at  any  place  on  the  runway.  Fortunately,  theorem  1  of  Appendix  B 
shows  that  the  optimal  solutions  can  be  found  by  searching  through  a  finite  set  of  potential 
turnoff  locations.  This  set  of  potential  solutions  consists  of  two  types  of  candidates  :  1)  primary 
and  2)  secondary  candidates.  Primary  candidates  are  the  best  exit  locations  for  each  aircraft. 
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which  are  found  during  the  simulations  of  individual  aircraft  landings.  Secondary  candidates 


are  exit  locations  for  each  aircraft  i  and  surface  condition  j  located  at  discrete  distance.  Dm,„. 
away  from  a  primary  candidate.  Secondary  candidates  are  generated  as  follows 

Ay  k  X  Dmln. 

where  Ay  is  the  primary  candidate  for  aircraft  i,  surface  condition  j. 
k  =  1,  2,  •  •  • 

Under  current  FAA  runway  turnoff  standards  [FAA,  1985)  two  adjacent  turnoffs  do  not  affect  the 
ru.,way  exit  index  unless  they  are  separated  by  at  least  750  ft.  (229  m.)  from  each  other,  in 
REDIM  we  have  added  more  flexibility  by  allowing  the  user  to  define  the  desired  minimum 
distance  Dmi„  between  adjacent  exits.  The  primary  and  secondary  candidates  comprise  the 
complete  candidate  set.  Suppose  a  runway  scenario  with  three  aircraft  in  consideration,  and 
the  best  locations  for  each  aircraft  and  two  pavement  conditions  are  (1000,  1100,  1300,1400, 
1600,  1700).  In  this  example,  the  complete  candidate  set  would  be  (1000,  1100,  1229,  1300, 
1329,  1400,  1458,  1529,  1558,  1600,  1629,  1687,  1700).  Theorem  1  of  Appendix  A  guarantees  that 
the  optimal  locations  should  be  some  of  the  candidate  set.  Thus  we  need  to  examine  only  the 
candidate  set  to  find  optimal  locations  instead  of  examining  infinite  points  on  the  runway. 


4.2  Estimation  of  Runway  Occupancy  Times 


The  simulations  of  landing  movements  provided  the  best  exit  locations  for  each  aircraft 
and  the  corresponding  times  required  to  reach  those  exit  locations.  The  complete  candidates 
were  enumerated.  The  next  step  is  to  find  out  the  time  required  to  clear  runway  for  every 
aircraft  for  every  candidate.  The  time  is  denoted  as  V  That  is,  T,lk  is  defined  as  the  runway 
occupancy  time  when  aircraft  i  takes  turnoff  candidate  k  on  surface  condition  j.  Every  T,(*  has 
three  components,  which  are  :  1)  time  to  reach  the  best  exit  location  (  Tb  ),  2)  time  to  travel 
from  the  best  location  to  the  candidate  (  7V, :  subscript  'fr'  stands  for  'f-ee  roll'.),  and  3)  time 
to  clear  the  runway  after  the  beginning  of  the  turnoff  (  Toff  )  Figure  4.3  illustrates  the 
components  of  the  Tijk- 
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Estimation  of  Runway  Occupancy  Tima  tor  Sacondaiy  Candldnios 


T„ 

-  Time  lo  reach  the  desired  exit  speed 

T,r 

-  Free  roll  time  to  reach  the  nearest  exit 

T.„ 

-  Time  in  iho  turnolt 

H*™. 

-  Threshold  crossing  alhludo 

Figure  4.3  The  Components  of  Tv*  Time. 


If  the  best  location  of  aircraft  i  is  farther  downrange  than  turnoff  candidate  k  from  the 
runway  threshold  ,  T,Jk  would  be  set  as  0,  which  means  aircraft  i  is  not  able  to  take  turnoff 
candidate  k.  Otherwise,  would  be  calculated  as  the  sum  of  T &  Tfr_  and  T0rr.  7^  is  calculated 
during  the  simulations.  Tlr  is  calculated  assuming  that  the  aircraft  coasts  on  the  runway 
'without  braking'  until  the  speed  of  the  aircraft  is  reduced  to  the  taxiing  speed  which  is 
specified  by  the  user.  'Without  braking'  implies  the  'oiling  friction  coefficient ,  f.  is  equal  to  0  03 
(  i.e.  deceleration  rate  is  9.81  m/s  2  (g)  *  0.03  (f)  =  0,2943  m/s  1  ).  Once  taxiing  speed  is 
reached,  the  aircraft  is  assumed  to  travel  on  the  runway  with  constant  speed  which  is  same 
as  taxiing  speed. 
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Suppose  Iho  bust  lurnoll  location  of  alia  oft  i  and  Th  an:  calculated  as  IlMK)  m  and  3'i 
seconds,  respectively  with  I  he  following  input  parameters  : 

•  Desired  exit  speed  =  30  m/s 

•  Taxiing  speed  -  7  m/s. 

The  distance  to  reach  taxiing  speed  would  be  : 

D  =  (  30  2-7J)/(2xgxf),  where  g  ^  9.81 ,  f=0.03 
=  1445  m. 

If  another  candidate  k  is  located  1300  m  downrange  from  the  runway  threshold,  ,  then  Tu,  the 
travel  time  to  reach  the  new  candidate  would  be  ; 

r  -  (2  y  D»> 

(Vi  +  Vt)  ' 

where  D,,  is  free  roll  distance  (  300  m  ), 

Vi  ic  initial  speed  (  30  m/s  ), 

V,  is  terminal  speed  {  Jv}  -  2  g  f  Dlr  =  26.3  m/s  ) 

=  10  5  seconds 

The  calculation  of  Toff  is  closely  related  to  the  turnoff  geometry  That  is,  Totl  is  the  travel 
time  along  the  turnoff  geometry  from  the  beginning  of  the  turn  to  the  clearance  of  the  runway. 
For  the  exact  calculation,  a  numerical  integration  requiring  large  computational  times  is 
needed.  Moreover,  this  integration  should  be  executed  for  every  V  unless  T,,k  is  set  as  0. 
T0„  is  therefore  approximated  by  the  method  described  in  Appendix  B,  to  reduce  the 
computation  time  T0„  usually  ranges  from  8  seconds  to  13  seconds  according  to  the  size  and 
exit  speed  of  the  aircraft  {for  moderate  exit  angles). 


4.3  Finding  Optimal  Locations 


In  this  section,  a  technique  to  find  optimal  turnoff  locations  is  described.  The  final  goal  is 
to  find  a  given  number  ol  turnoff  locations  which  minimize  the  total  weighted  sum  of  ROT  from 
the  set  of  candidates  The  number  of  turnoff  is  provided  by  the  user  The  optimization  task  can 
be  modeled  as  a  specific  linear  programming  model.  A  dynamic  programming  algorithm  is 
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applied  to  And  the  solution,  since  the  dynamic  programming  algorithm  is  more  efficient  than 
the  ordinary  linear  programming  algorithm  for  our  case. 


4.3.1  Mathematical  Model 


Suppose  M  different  types  of  aircraft  use  a  runway,  then  2  M  different  turnoff  locations 
would  be  calculated  for  every  aircraft  and  two  runway  surface  conditions  (  dry  and  wet )  during 
the  simulation  of  landing  movements.  A  complete  set  of  exit  candidates,  which  is  indexed  k 
=  1  to  K,  is  generated  based  on  the  2  M  initial  locations.  It  is  not  always  permissible  to  assign 
aircraft  i  on  surface  condition  j  to  candidate  k.  Let  us  define  A(i,j)  as  a  set  of  feasible 
candidates  for  aircraft  i  on  the  surface  condition  j,  for  i  =  1  to  M,  j  =  1  to  2  If  candidate  k  is 
nearer  from  the  threshold  than  the  primary  candidate  for  aircraft  i  and  surface  condition  j,  the 
candidate  k  does  not  belong  to  A(i,j). 

If  exit  candidate  k  is  selected  to  be  built,  the  candidates  which  are  within  „  1223  m  or 
750  ft  )  from  the  candidate  k  can  not  be  constructed.  Let  us  define  S(k)  as  the  it /.a  / 
exclusive  set  of  candidates  in  which  at  most  one  candidate  can  be  selected  to  be  built,  for  k  - 
1  to  K. 

In  order  to  the  minimize  the  weighted  sum  of  ROT,  information  about  weights  should  be 
provided  by  the  user.  Let  a,  be  the  proportion  of  aircraft  i  for  i  =  1  to  M,  and  let  p,  be  the 
probability  of  occurrence  of  the  surface  condition  for  j  =  1  to  2  (  if  j  =  1,  surface  condition  is 
dry.  otherwise,  surface  condition  is  wet  ). 

Suppose  the  number  of  exits  to  be  built  is  set  as  N.  The  binary  decision  variables  are 
defined  as  follows: 

1,  if  exit  candidate  k  is  selected 

x*  = 

0.  otherwise,  for  k  =  1  to  K 
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1.  if  aircraft  i  is  assigned  to  the  exit  candidates  k  on  surface  condition  j 


y*/*  - 

0,  otherwise,  for  i  =  1  to  M,  j  =  1  to  2,  k  eA(ij) 

Then,  the  model  which  attempts  to  design  a  feasible  runway  with  the  least  total  weighted 
runway  occupancy  time  may  be  formulated  as  follows: 

M2 

Minimize  Z  Z  Z  a'  pi  T'tk  '4  V 

l  =  1;=  U(eA(ij) 


subject  to  Z  Yijk  —  1 

k 

for  i  =  1,2 M ;  y  =  1,2 

{4.2} 

E «» s  1 

keS(k) 

for  /c  =  1,2,...,K 

{43} 

K 

Z*S" 

1 

{4  4} 

y,kj  ^  *k 

for  i  =  1,2 . M,j  =  1,2;  keA(ij) 

{4.5} 

*.y 

binary 

{46} 

The  objective  function  (Eq.  4.1)  represents  the  aggregate  expected 

runway  occupancy 

time.  Constraint  (Eq.  4.2)  requires  that  each  aircraft  type  should  be  assigned  to  one  (available) 
exit  under  each  surface  condition.  Constraint  (Eq.  4.3)  ensures  a  feasible  mix  of  exits,  while 
constraint  (Eq.  4.4)  enforces  a  maximum  limit  to  the  total  number  of  exits  constructed.  The 
fourth  constraint  (Eq.  4.5)  asserts  that  only  the  constructed  exits  must  be  used,  and  lastly,  Eq. 
4.6  enforces  the  logical  restrictions  on  the  .uriables. 
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The  same  formulation  given  above  may  be  used  to  model  the  problem  of  re-designing  or 
modifying  existing  runways,  by  simply  fixing  the  appropriate  variables  x*  to  be  one.  This 
option  can  also  be  adopted  for  a  priority  enforcing  choice  of  certain  exits. 


4.3.2  Dynamic  Programming  Formulation 


Suppose  the  number  of  exits  to  be  built  is  N,  the  number  of  candidates  is  K,  and  the 
candidates  are  sorted  based  on  the  distance  from  the  threshold.  For  the  dynamic 
programming  (DP)  formulation,  one  imaginary  candidate  need  to  be  introduced.  This 
imaginary  candidate  is  indexed  0,  and  is  located  229  m.  ahead  of  the  first  candidate.  The 
corresponding  TlJ0  is  set  as  0,  for  all  (i.j).  This  means  no  aircraft  can  take  exit  0  With  the 
imaginary  candidate,  we  can  observe  the  following  characteristics  of  T*: 

1)  There  exists  at  least  one  Tijk  =  0  for  all  (i.j). 

2)  If  Tl/ka  >  0,  then  T,lk  >  0,  for  A  £  A0.  for  alt  (i.j). 

3)  T,,„  <,  r,;(*  +  1),  fork  e  A{iJ),  for  all  (i.j). 

D(k)  is  defined  as  the  distance  from  candidate  k  to  candidate  1.  Thjs  D(0)  =  -229  m  ,  and 
0(1)  =  0  m.  Let  us  define  another  variable,  which  is  denoted  as  K0,  a  candidate  index  beyond 
which  at  least  one  exit  should  be  constructed.  K0  is  determined  by  : 

K0  =  Max  {  k  ;  TiJ{k ..  „  =  0  for  some  (i.j)  }  <,  K 

K0  ensures  that  each  aircraft  will  be  assigned  to  an  exit,  even  if  it  is  the  largest  aircraft.  With 
the  variables  defined  in  the  previous  section  and  above,  the  DP  formulation  is  as  follows  : 

Stages  ;  Stage  q  corresponds  to  a  situation  in  which  up  to  q  exits  can  be  located  to  the 
right  of  the  the  last  exit  already  located,  q  ranges  from  1  to  N.  For  1  <  q  <  N, 
(N-q)  exits  are  assumed  to  have  been  constructed. 

States  ;  The  state  s,  at  stage  q  is  a  candidate  index,  and  corresponds  to  the  right  most 
exit  currently  located.  For  1  <  q  <  N  -  1,  the  possible  values  of  s„  are  /,,••• 

,  K,  where  /,  is  the  smallest  exit  candidate  index  such  that  it  is  possible  to 
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construct  (N-q)  exits  in  candidate  1 subject  to  the  Dmin  separation 
restriction.  That  is,  /,  is  determined  by  : 

/,  =  Min  {  k  ;  D(k)  ^  ( N-q)x  Dmln  },  and 
If  q  =  N,  sN  =  0. 

Decisions  ;  Decision  dq  is  another  candidate  index.  Given  stage  q  and  state  s,,  the 
decision,  dq  ,  corresponds  to  the  next  exit  to  be  constructed  to  the  right  of  s„. 
Let 'dq  =  O'  mean  that  no  more  exits  are  constructed.  Then  the  possible  values 
of  d,  are  0,  and  L„,  •••,  K,  for  1  <,  q  ^  A/  —  1,  where  L,  is  the  smallest  exit 
index  such  that  D(L,)  -  D(/„)  >  Dm i„,  if  it  exists,  for  q  =  N,  LN  =  1. 

Given  any  stage  q  and  state  sq,  all  aircraft-surface  condition  combinations  (i,j)  for  which  T,ISq 
>  0  would  have  been  assigned  to  some  existing  exit,  due  to  the  characteristics  of  T,;*  Hence, 
the  problem  decomposes  into  locating  up  to  q  more  exits  to  Ihe  right  of  s,  with  the  minimum 
separation  constraint,  considering  only  (i,j)  combinations  whose  0-  which  implies  that 

(i,j)  is  not  yet  assigned.  Since  the  optimum  of  this  decomposed  problem  is  independent  of  the 
previous  decision,  and  depends  only  on  q  and  sq,  Bellman's  principle  of  optimality  holds,  and 
thus,  the  DP  application  is  valid. 

With  the  stage,  the  state,  and  the  decision  defined  above,  some  functions  need  to  be 
defined  for  the  complete  DP  formulation.  These  are  : 

Immediate  return  function 

The  return  function  cq(sq,dq)  is  the  'immediate'  stage  cost  incurred  by  making  decision, 
at  stage  q  in  state  sq.  This  cost  corresponds  to  the  additional  (i,j)  assignments  which  can  be 
made  with  a  given  dq.  Hence, 
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oo, 


if  D{dq)  -  D{sq)  <  Dmln,  and  dq  *  0 


Cg(Sq,  dq)  — 


<W):  T, 


I 


ai  Pi  TH<Jq 


if  D(dq)  -  O(Sq)  >  Dmln,  and  dq  *  0 


'l‘o 


:  0,  but  TIJd  >  0} 


{4.7) 


0,  if  dq  =  0 


Stage  transition  function 

q  -  1 ,  if  dqJ=  0 

*q(dq)  {4  8} 

0,  if  dq  =  0 

State  transition  function 

dq,  if  dq  #  0 

Tq(s<?'  <*q)  {4.9} 

Sq,  if  dq  =  0 

Recursive  formula 

Defining  /^(s,)  to  be  the  optimal  accumulated  return  function  with  given  input  state  s,  at  stage 
q.  the  recursive  formula  would  be: 


fq(sq)  =  minimum  {  cq{sq,dq)  +  Tq(Sq,  d„) )  },  {4.10} 

where  the  final  condition  is 

OO,  if  6q  Kq 

W«o)  =  (4in 

0,  otherwise 

By  iterating  the  recursive  formula  (  4,10  )  with  q  from  1  to  (N-1).  we  can  find  the  optimal 
accumulated  return  (  minimum  weighted  sum  of  ROT  )  for  all  possible  states  for  each  stage. 
At  the  final  iteration,  or  the  last  stage  (  q  =  N  ),  the  overall  weighted  sum  of  ROT  is  minimized, 
and  then  a  sequence  of  optimal  decisions,  cTq  (  q  =  1,  •  •  ,  N  ),  v.hich  minimizes  the  overall 
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weighted  sum  of  ROT  is  revealed.  These  c^'s  are  the  optimal  exit  candidate  indices  which  we 
are  looking  for. 


4.4  Aircraft  Exit  Assignment 


By  the  DP  technique,  the  optimal  exit  locations  are  found.  The  final  step  in  the 
optimization  sequence  is  to  assign  every  aircraft-surface  condition  combination,  (i.j),  to  an 
appropriate  exit.  This  step  is  performed  by  making  (i.j)  to  take  the  exit  which  is  permissible 
and  requires  minimum  ROT. 


4.5  A  Simple  Example 


In  this  section,  a  simple  example  is  discussed  to  illustrate  the  optimization  procedure 
developed  previously.  Suppose  three  aircraft  use  a  single  runway,  wnere  two  exits  will  be 
constructed.  The  aircraft  types  and  the  relative  frequncies  of  operation  are  :  1)  Learjet-31 
(30%),  2)  Airbus  A300-600  (30%),  and  3)  Boeing  B767-300  (40%).  The  desired  exit  speed  for 
all  the  three  aircraft  is  30  m/s  (  67  MPH  ).  The  exit  reliability  factor  is  90%.  The  chances  of 
dry  and  wet  conditions  occcuring  are  same  (i.e.,  50%  each).  With  these  data,  the  simulation 
of  landing  movement  calculates  six  primary  exit  candidates  for  three  aircraft  and  two  runway 
surface  conditions.  That  is,  the  best  exit  location  for  the  Learjet-31  on  dry  surface  is  906  m, 
and  on  wet  surface  is  968  m.  1546  m  is  the  best  location  for  the  A300-600  on  dry  surface,  and 
1711  m  is  the  best  location  under  wet  condition.  The  best  location  for  B767-300  is  1638  m  on 
dry  surface,  and  1816  m  on  wet  surface.  The  corresponding  ROT's  are  26.9,  28.5,  40.8,  44.6, 
42.3,  46.4  seconds,  respectively,  as  seen  on  Table  4.1. 

Based  on  the  primary  candidates  and  assuming  an  arbitrary  minimum  separation 
distance  (Dm)„)  of  213  m.  (700  ft  ),  eight  more  secondary  candidates  are  found  to  comprise  the 
complete  candidate  set.  These  fourteen  exit  candidates  are  are  906m,  968m,  1119m,  1181m, 


Optimization  Module 


75 


1332m,  1394m,  1545m,  1546m,  1607m,  1638m,  1711m,  1758m,  1795m,  and  1816m  (  STEP  1  ).  A 
TIJk  matrix  is  calculated  as  shown  in  Table  4.1  (  STEP  2  ).  The  optimization  is  performed  with 
the  7//*  data,  and  then  968m  and  1816m  are  selected  as  optimal  exit  locations  (  STEP  3  ). 
Finally,  the  aircraft  are  assigned  to  the  selected  exit  locations  as  shown  in  Table  4  2,  the 
weighted  average  ROT  is  calculated  as  43.2  seconds. 


4.6  Modified  Algorithm  for  ' Improvement '  Analysis 


The  optimization  algorithm  described  in  section  4.1  to  4.4  was  developed  for  design 
analysis  which  assumed  no  exits  were  available  on  the  runway.  With  some  modifications,  this 
algorithm  can  be  applied  to  an  improvement  analysis  scenario  in  which  some  exits  already 
exist  on  the  runway  and  a  few  more  exits  will  be  added  to  reduce  the  ROT. 

In  this  new  procedure,  the  existing  exit  locations  as  well  as  the  best  locations  are 
considered  as  primary  candidates.  The  complete  candidates  are  generated  with  the  same 
principles  used  in  design  analysis,  and  then  the  candidates  which  are  located  within  the 
±  Dmm  range  from  the  existing  exits  are  eliminated.  Stages,  states,  and  decisions  of  DP 
formulation  are  same  as  those  of  the  design  analysis.  The  immediate  return  function  should 
be  changed  to  consider  the  effect  of  the  existing  exits.  Suppose  the  ROT  of  aircraft  i  should 
be  accumulated  as  an  immediate  return  of  a  decision,  dq,  associated  with  a  state,  s,.  If  there 
are  some  existing  exits  in  the  region  of  {  sq,  dq  ),  and  a  existing  exit  requires  less  ROT  than 
dq  does,  then  the  less  ROT  required  by  the  existing  exit  is  considered  as  an  immediate  return. 


Optimization  Module 


76 


Optimization  Module 


77 


5.0  Output  Module 


The  Output  Module  plays  a  very  significant  role  in  the  program.  REDIM  is  structured  so 
that  at  each  stage  the  user  is  prompted  with  specific  questions  and  guidelines  that  are  to  be 
followed  in  order  to  view  the  appropriate  results.  A  significant  feature  in  this  module  is  that 
in  each  of  the  runmodes,  all  the  screens  clearly  display  information  regarding  the  aircraft  and 
the  airport  data.  This  information  provided  at  the  top  of  each  of  the  screens,  will  be  of  great 
use  to  the  user  as  (s)he  will  be  presented  with  the  general  information  pertaining  to  that  type 
of  analysis.  The  Output  Structure  is  shown  in  terms  of  a  Output  functional  flow  diagram  in  Fig 
5  1. 


After  defining  the  variables  in  the  Input  Module,  the  user  may  go  back  to  the  Main  Menu. 
From  this  menu,  the  user  has  the  option  to  go  to  the  Output  Menu  among  other  options.  When 
the  user  chooses  to  go  to  the  Output  Menu,  (s)he  is  provided  with  four  options.  The  options 
being:  (1)  View  the  Output  (2)  Print  the  Report  (3)  Help  and  (4)  Go  to  Main  Menu.  Fig  5  2 
clearly  depicts  the  Output  Menu  on  the  screen.  In  the  following  paragraphs,  all  the  options 
that  are  provided  to  the  user  in  the  Output  Menu  will  be  discussed  in  detail. 


5.1  View  the  Output 


When  the  user  selects  this  option,  the  program  automatically  goes  into  the  runmode 
variable  that  was  previously  specified  by  the  user  in  the  Input  Module  and  present  the  View 
Menu  screen  which  corresponds  to  the  aforesaid  option.  For  each  of  the  runmode  variable 
options,  the  program  presents  different  screens  which  prominently  display  the  results  in  color, 
while  at  the  same  time  guiding  the  user.  Each  screen  also  displays  the  general  input 
variables  that  were  initially  provided  by  the  user  in  the  Input  Module.  The  user  at  any  stage 
in  the  Output  Module  is  allowed  to  go  back  to  the  Main  Menu,  through  which  a  wide  range  of 
options  can  be  chosen.  Now  we  will  go  into  each  of  the  runmode  variable  option  and  see 
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Figure  5.1  Output  Functional  Flow  Diagram. 


OUTPUT  KflU 


Prei*  escape  (Esc)  key  to  return  to  WIN  ftrt). 


Use  arrou  (t  or  l)  key  to  change  your  choice. 
Press  enter  (  *J  1  key  if  the  choice  is  correct. 


DUieu  the  Output] 


2 (Print  The  Report 
3)Help 

f)Go  To  Hi  in  Menu 


Figure  5.2  Output  Menu. 


explicitly  as  to  what  is  presented,  and  also  the  extent  of  flexibility  that  is  made  available  to  the 
user. 

5.1.1  Evaluate  an  Existing  Runway 


This  runmode  is  specifically  designed  to  evaluate  and  existing  runway.  The  data  required 
to  evaluate  an  existing  runway  is  entered  by  the  user  in  the  Input  Module.  As  the  user  selects 
to  view  the  Output  from  the  Output  Menu  (Fig.  5.2),  (s)he  is  provided  with  a  different  screen 
which  displays  the  View  Menu  (Fig.  5.4).  This  View  Menu  corresponds  to  the  runmode  option 
which  is  used  to  evaluate  an  existing  runway.  Three  options  are  provided  in  this  menu  The 
options  being  (1)  ROT  /  REL  Table  of  Results  (2)  Turnoff  locations  and  their  Geometries  and 
(3)  Go  to  Output  Menu.  The  (low  pattern  for  this  runmode  option  is  shown  in  Fig.  5.3  and  the 
View  Menu  screen  is  illustrated  in  Fig.  5.4.  In  the  following  paragraphs,  we  will  discuss  in 
detail  these  options. 
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ROT  /  Reliability  Table  o(  Results:  This  screen,  shown  in  Fig.  5.5,  displays  the  number  of 
turnoffs,  their  location,  and  the  type  of  turnoffs.  The  FAA  standard  turnoffs,  30  degree,  45 
degree  and  90  degree  are  illustrated  in  Figs.  5.7,  5.8  and  5.9,  respectively.  For  each  of  the 
aircraft  selected  earlier  in  the  Input  Module,  the  corresponding  runway  occupancy  time  and 
the  reliability  associated  with  each  of  the  turnoffs  is  displayed  These  help  the  user  in  judging 
the  appropriate  turnoff  for  that  particular  aircraft.  One  more  significant  feature  in  this  table  is 
the  separation  of  the  values  for  wet  and  dry  airport  conditions.  The  user  is  also  in  a  position 
to  view  the  changes  in  runway  occupancy  time  and  the  reliability  associated  with  each  of  the 
turnoffs  when  two  runway  conditions  (wet  or  dry)  are  present.  In  addition  to  these,  the  relative 
frequency  of  occurrence  for  every  aircraft  under  every  runway  scenario  is  specified  as  a 
percentage  of  all  the  aircraft  occurrences.  This  table  is  shown  in  Fig.  5.5. 

In  each  screen  as  only  four  aircraft  are  shown,  it  is  necessary  to  press  'F  key  (for  forward) 
to  view  additional  aircraft  and  'B'  key  (for  backward)  to  view  the  previous  screen.  In  each  of 
the  screens  the  viewer  is  presented  with  the  average  runway  occupancy  time  which 
encompasses  the  whole  population  of  aircraft  selected  by  the  user.  These  include  dry  and 
wet  conditions  at  the  airport.  An  option  is  provided  to  print  the  table  by  just  pressing  the  'P' 
key  (for  print).  The  user  is  also  allowed  to  go  back  to  the  View  Menu  at  any  stage  by  entering 
V'  key  (for  View  Menu). 

Turnoff  Locations  and  their  Geometries:  In  this  screen,  shown  in  Fig.  5  6,  t  the  user  is 
presented  with  the  display  of  turnoff  locations  along  with  other  pertinent  airport  data.  The 
turnoff  locations  along  the  runway  downrange  are  shown  in  Fig.  5.6.  An  added  feature  is  the 
presentation  of  exit  numbers  and  their  location  in  a  tabular  form.  An  option  *o  view  the 
standard  FAA  geometry  for  each  of  the  turnoffs  is  also  provided.  The  user  is  prompted  to 
enter  the  exit  number  to  view  the  standard  FAA  geometry.  Each  of  the  turnoffs  may  represent 
any  of  the  standard  FAA  geometries  viz.,  30  degree  or  45  degree  or  90  degree.  Only  one  exit 
number  is  to  be  entered  to  view  the  complete  turnoff  geometry.  As  the  exit  number  is  entered, 
the  program  determines  the  type  of  turnoff  for  the  exit  number  from  the  Input  Module  Earlier, 
in  the  Input  Module,  the  user  selected  the  type  of  turnoff  geometry  for  a  particular  exit  number. 
This  turnoff  geometry  may  represent  any  of  the  three  standard  FAA  turnoff  geometries  shown 
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in  Figs.  5.7,  5.8  and  5.9  respectively.  After  viewing  the  standard  FAA  turnoff  geometry  for  the 
specified  turnoff,  the  user  is  returned  to  the  first  screen  where  a  choice  for  a  different  exit  can 
be  made.  This  screen  would  appear  as  shown  in  Fig.  5.5.  The  user  may  exit  from  this  screen 
by  entering  '0'  (zero). 

When  the  user  selects  to  view  the  standard  turnoff  geometry  for  a  particular  turnoff,  the 
program  automatically  takes  in  the  values  of  runway  width,  taxiway  width,  the  distance 
between  the  runway  and  taxiway  and  the  type  of  turnoff.  This  data  was  earlier  supplied  by  the 
user  in  the  Input  Module.  In  the  standard  FAA  turnoff  geometry  display,  the  specifications  are 
prominently  shown  for  the  benefit  of  the  user.  Although,  the  Metric  system  is  mainly  used  for 
computations  in  the  program,  the  units  for  specifications  are  also  displayed  in  the  English 
system.  This  helps  the  user  Who  might  be  still  using  the  FAA  standards  in  the  English  system. 

Go  to  Output  Menu:  This  option  is  provided  to  enable  the  user  to  go  back  to  the  Output  Menu 
from  the  View  Menu.  The  user  may  also  exercise  this  option  by  pressing  the  escape  'Esc'  ke / 


5.1.2  Improve  an  Existing  Runway 


This  runmode  option  is  made  available  to  the  user  to  improve  or  modify  an  existing 
runway.  The  required  data  of  the  airport  facility  that  needs  to  be  improved  is  entered  in  the 
model  through  the  Input  Module.  The  user  after  entering  the  data  may  go  to  the  Output  Menu 
through  the  Main  Menu.  The  Output  Menu  screen  is  as  shown  in  Fig.  5.2.  and  the  flow  pattern 
for  this  runmode  option  is  shown  in  Fig.  5.10.  As  the  user  selects  to  view  the  output  from  the 
Output  Menu,  he  or  she  is  presented  with  the  View  Menu  screen.  This  View  Menu  screen 
corresponds  to  the  runmode  option  which  is  to  improve  an  existing  facility.  The  View  Menu 
screen  is  provided  with  four  options:  (1)  ROT  Table  of  Results,  (2)  Turnoff  Locations  and  their 
Turnoff  Geometries,  (3)  Aircraft  Statistics  and  their  ROT's,  and  (4)  Go  to  Output  Menu.  The 
view  menu  screen  is  shown  in  Fig.  5.11.  In  the  following  paragraphs,  we  will  discuss  in  detail 
each  one  of  the  above  said  options. 
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Run  Mod*  »  2 

Redesign  of  an  Existing  Runwsy 
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igram  for  Runmoi 


ROT  Table  of  Results:  This  table  displays  the  number  of  turnoffs,  the  existing  and  the 
proposed  ones,  and  their  locations  in  a  format  similar  to  that  of  Fig.  5.12.  The  type  of  turnoff 
geometry  is  also  displayed  in  the  third  row  of  the  table  The  type  of  turnoff  geometry  either 
could  be  of  standard  FAA  turnoff  geometry  viz.,  30  degree  or  45  degree  or  90  degree,  for  the 
existing  ones  or  a  new  turnoff  geometry  developed  by  the  model  for  the  proposed  new  ones. 
The  new  turnoffs  are  differentiated  from  the  existing  ones  through  color  coding  This  helps  in 
knowing  the  location  of  new  turnoffs  at  a  simple  glance.  For  each  of  the  aircraft  selected 
earlier  in  the  Input  Module,  the  corresponding  runway  occupancy  time  associated  with  the 
appropriate  turnoff  is  displayed.  The  relevant  values  for  both  the  airport  conditions,  wet  and 
dry  are  displayed.  The  user  is  also  in  a  position  to  view  the  change  in  runway  occupancy  time 
associated  with  appropriate  turnoff  when  the  airport  conditions  (wet  or  dry)  are  changed.  In 
addition  to  these,  the  relative  frequency  of  occurrence  of  this  aircraft,  for  a  particular  scenario 
(wet  or  dry  conditions)  is  specified  as  a  percentage  of  all  the  aircraft  occurrences.  The  emp, 
boxes  in  the  table  imply  that  the  aircraft  in  question  cannot  negotiate  that  particular  turnoff. 
The  reliability  associated  with  the  aircraft  and  turnoff,  which  was  earlier  provided  by  the  user 
in  the  Input  Module,  is  displayed  at  the  bottom  of  the  table 

In  each  screen,  as  only  four  aircraft  are  shown  in  the  table,  it  is  necessary  to  press  'F' 
key  (for  forward)  to  view  additional  aircraft  and  'B'  key  (for  backward)  to  view  the  previous 
screen.  In  each  of  the  screens  the  user  is  presented  with  the  average  runway  occupancy  time 
which  encompasses  the  whole  population  of  aircraft  selected  by  the  user.  These  include  dry 
and  wet  conditions  at  the  airport.  The  user  is  provided  with  an  option  to  print  the  table  by  just 
pressing  the  'P'  key  (for  print).  The  user  is  also  allowed  to  go  back  to  the  View  Menu  at  any 
stage  by  entering  'V'  key  (for  View  Menu).  Figure  5.12  illustrates  and  ROT  table  with  four 
aircraft. 

Turnoff  Locations  and  their  Geometries:  In  this  screen  (Fig.  5.5),  the  user  is  presented  with  a 
graphical  display  of  turnoff  locations  along  the  runway  downrange.  In  the  display,  the  existing 
turnoffs  as  well  as  the  proposed  turnoffs  are  shown.  The  new  turnoffs  are  differentiated  from 
the  existing  ones  through  color  coding.  This  enables  the  user  in  noticing  the  location  of  new 


Output  Module 


88 


Output  Module 


89 
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Figure  5.13  Comparison  of  Centerline  Turnoff  Geometries. 


turnoffs  along  the  runway  downrange.  An  added  feature  in  this  screen  is  the  presentation  of 
exit  numbers  and  their  location  in  a  tabular  form. 

On  the  same  screen,  depicted  in  Fig  5.5,  a  small  menu  screen  is  provided  with  options 
to  view  the  complete  turnoff  geometry  or  to  compare  the  centerline  geometries.  For  the  latter 
option,  the  user  needs  to  input  the  exit  number(s).  The  user  has  the  flexibility  to  select  any 
combination  of  exit  numbers  for  comparing  the  geometry  of  different  turnoffs  or  may  even 
select  only  one  exit  number  to  view  the  centerline  geometry  of  a  particular  turnoff.  Fig.  5.13 
illustrates  the  output  screen  for  the  comparison  of  turnoff  geometries.  It  should,  however,  be 
noted  that  all  the  requested  turnoff  geometries  start  at  a  common  point  for  comparison 
purposes.  Also,  as  can  be  seen  from  Fig.  5.15,  the  user  has  the  benefit  to  view  the  centerline 
turnoff  geometry  coordinates  in  a  tabular  form  for  each  of  the  turnoffs.  The  user  needs  to 
press  'T  to  view  the  table  with  coordinates  of  centerline  turnoff  geometry  from  the  screen 
which  displays  the  centerline  turnoff  geometry.  If  the  user  does  not  intend  to  view  the 
coordinates,  (s)be  may  press  'E'  (for  exit  location)  to  return  to  the  turnoff  locations  screen. 
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The  user  after  returning  to  tho  previous  screen,  which  depicts  the  small  menu,  may  now 
view  the  centerline  turnoff  geometry  or  the  complete  turnolt  geometry  Tho  user  may  either 
select  to  view  the  complete  turnoff  geometry  or  may  even  again  choose  to  view  the  center  line 
turnoff  geometry  for  a  different  combination  of  turnoffs  When  the  user  select*  to  view  the 
complete  turnoff  geometry,  the  program  requests  the  user  to  enter  the  exit  number  to  be 
displayed.  Here,  the  user  needs  to  input  only  one  exit  number  and  not  any  combination  of 
exit  numbers.  The  program  displays  the  complete  turnoff  geometry  of  that  specific  exit 
number  entered  by  the  user.  For  the  exit  number  entered  by  the  user,  the  program  interna  / 
determines  the  type  of  turnoff  geometry.  The  type  of  turnoff  geometry  can  be  either  of  the  30 
degree  or  45  degree  or  90  degree  standard  FAA  turnoff  for  the  existing  turnoffs  or  can  be  a 
new  turnoff  geometry  developed  by  the  model  for  the  proposed  new  turnoff(s).  The  FAA 
standard  30  degree,  45  degree  and  90  degree  are  best  illustrated  in  Figs.  5.7,  5  8  and  5  9 
respectively.  The  new  turnoff  geometry  developed  by  the  model  for  the  proposed  new  turnoff 
is  shown  in  Fig.  5.16.  The  user  after  viewing  the  complete  turnoff  geometry  of  the  turnoff 
requested,  may  return  back  to  the  turnoff  location  and  their  geometries  screen  where  (s)he 
may  again  enter  a  different  turnoff  to  view  its  complete  turnoff  geometry. 

When  the  user  selects  to  view  the  standard  turnoff  geometry  for  a  particular  turnoff,  the 
program  automatically  takes  in  the  values  of  runway  width,  taxiway  width,  the  distance 
te'ween  the  runway  and  taxiway,  the  type  of  turnoff  and  other  data  supplied  in  the  Input 
Module.  In  the  standard  FAA  turnoff  geometry  display,  the  specifications  are  prominently 
shown,  and  the  units  are  displayed  in  both  English  and  Metric  systems  for  the  benefit  of  the 
user. 

Aircraft  Statistics  and  their  ROT's:  This  screen  displays  the  runway  occupancy  time  of  each 
of  the  aircraft  in  the  form  of  a  bar  chart.  At  the  top  edge  of  each  of  the  bars,  the  aircraft 
number  is  displayed.  For  each  of  the  aircraft,  selected  by  the  user  in  the  Input  Module,  the 
runway  occupancy  time  is  separately  displayed  for  wet  and  dry  runway  surface  conditions. 
For  the  benefit  of  user,  the  bars  representing  the  aircraft  and  its  name  are  of  the  same  color. 
On  the  bar  chart,  a  straight  line  is  drawn  across  the  bars  to  portray  the  average  ROT 
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In  each  screen  a  maximum  of  six  aircraft  are  shown,  each  with  two  bars,  one  tor  wet 
condition  and  the  other  for  dry  condition.  It  is  necessary  to  press  F'  key  (for  forward)  to  view 
additional  aircraft  and  'B'  key  (for  backward)  to  view  the  previous  screen  In  each  of  the 
screens,  the  average  runway  occupancy  time  which  encompasses  the  whole  population  of 
aircraft  selected  by  the  user  is  shown.  This  includes  dry  and  wei  conditions  at  iiie  airport 
An  option  is  provided  to  print  the  table  by  just  pressing  the  'P'  key  (for  print).  The  user  is  also 
allowed  to  go  back  to  the  View  Menu  at  any  stage  by  entering  'V'  key  (for  View  Menu)  The 
bar  chart  is  illustrated  in  Fig  5  17. 

Go  to  Output  Menu:  This  option  is  provided  to  enable  the  user  to  go  back  to  the  Output  Menu 
from  the  View  Menu.  The  user  may  also  exercise  this  option  by  pressing  the  escape  Esc'  key. 

5.1.3  Design  of  a  New  Runway 

i 

This  runmode  deals  with  the  design  of  a  new  runway  facility.  The  user  is  requested  to 

enter  the  relevant  data  in  the  Input  Module  for  the  design  of  a  new  runway  facility  The  user 

'  .  1 

after  entering  the  data  in  the  Input  Module,  is  returned  to  the  Main  Module  The  user  at  this 
point  may  go  to  the  Output  Menu  to  view  the  results  The  Output  Menu  is  as  shown  in  Fig 

i 

5.2.  and  the  flow  pattern  corresponding  to  this  runmode  option  is  shown  in  Fig.  5.15.  As  the 
user  selects  to  view  the  output  from  the  Output  Menu,  he  or  she  is  presented  with  a  new 
screen  which  displays  the  View  Menu  screen.  This  View  Menu  screen  corresponds  to  the 
runmode  option  which  is  to  design  a  new  runway  facility.  Four  options  are  provided  in  the 
View  Menu.  The  options  are  (1)  ROT  Table  of  Results  (2)  Turnoff  Location  and  their 
Geometries  (3)  Aircraft  Statistics  and  their  ROT's,  and  (4)  Go  to  Output  Menu.  The  View  Menu 
screen  is  best  illustrated  in  Fig.  5.11.  In  the  following  paragraphs,  we  will  discuss  in  detail 
each  of  these  options. 

ROT  Table  of  Results:  This  table  displays  the  number  of  new  turnoffs  for  the  proposed  new 
runway  facility  and  their  locations.  As  none  of  the  standard  FAA  turneffs  are  proposed  in  the 
new  runway,  the  type  of  turnoff  geometry  is  described  as  a  new  turnoff  geometry  developed 
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Run  node  *  3 

Design  of  a  New  Runway  facility 
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by  the  model.  For  each  of  the  aircraft  selected  earlier  in  the  Input  Module  by  the  user  the 
corresponding  runway  occupancy  time  associated  with  the  appropriate  turnoff  is  displayed 
These  values,  for  both  wet  and  dry  airport  conditions,  are  displayed  separately  in  the  table. 
In  addition  to  these,  the  relative  frequency  of  occurrence  of  this  aircraft,  for  a  particular 
scenario  (wet  or  dry  conditions)  is  specified  as  a  percentage  of  all  the  aircraft  occurrences 
The  empty  boxes  in  the  table  imply  that  the  aircraft  in  question  cannot  negotiate  that  particular 
turnoff.  The  reliability  associated  with  each  aircraft  and  turnoff(s),  which  was  earlier  provided 
by  the  user  in  the  Input  Module,  is  displayed  at  the  bottom  of  the  table. 


In  each  screen,  as  only  four  aircraft  are  shown  in  the  table,  it  is  necessary  to  press  F' 
key  (for  forward)  to  view  additional  aircraft  and  'B'  key  (for  backward)  to  view  the  previous 
screen  In  each  cf  the  screens  the  average  ROT  which  encompasses  the  whole  population 
of  aircraft  selected  by  the  user  is  shown.  This  includes  dry  and  wet  conditions  at  the  airport. 
An  option  to  print  the  table  by  just  pressing  the  'P'  key  (for  print)  is  also  provided  The  user 
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is  also  allowed  to  go  back  to  the  View  Menu  at  any  stage  by  entering  'V'  key  (for  View  Menu). 
This  table  is  illustrated  in  Fig.  5.12. 

Turnoff  Locations  and  their  Geometries:  This  screen  presents  the  user  with  a  graphical 
display  of  optimal  turnoff  locations  along  the  runway  downrange.  These  optimally  located 
turnoffs  are  all  proposed  by  the  model  for  the  new  runway  facility.  A  table  with  the  exit 
numbers  and  their  location  is  also  presented  on  the  same  screen  The  output  screen  is 
similar  to  that  shown  in  Fig.  5.5. 

On  the  same  screen,  depicted  in  Fig.  5.5,  a  small  menu  screen  is  provided  with  options 
to  view  the  complete  turnoff  geometry  or  to  compare  the  centerline  geometries.  For  the 
comparison  of  centerline  geometries,  the  user  needs  to  input  the  exit  number(s).  This 
provides  flexibility  to  select  any  combination  of  exit  numbers  for  comparing  the  geometry  of 
different  turnoffs  or  may  even  select  only  one  exit  number  to  view  the  centerline  geometry  of 
a  particular  turnoff.  Fig.  5.11  illustrates  the  output  screen  for  the  comparison  of  turnoff 
geometries.  It  should,  however,  be  noted  that  all  the  turnoff  geometries  start  at  a  common 
point  for  comparison  purposes.  Also,  as  can  be  seen  from  Fig.  5.15,  the  user  has  the  benefit 
to  view  the  centerline  turnoff  geometry  coordinates  in  a  tabular  form  for  each  of  the  turnoffs. 
The  user  needs  to  press  'T  to  view  the  table  with  coordinates  of  centerline  turnoff  geometry 
from  the  screen  which  displays  the  centerline  turnoff  geometry.  If  the  user  does  not  \r.'H~  ", 
view  the  coordinates,  (s)lie  may  press  'E'  (for  exit  location)  to  return  to  the  turnoff  locaf  or, s 
screen. 

The  user  after  returning  to  the  previous  screen,  which  depicts  the  small  menu,  may  now 
view  the  centerline  turnoff  geometry  or  the  complete  turnoff  geometry.  The  user  may  either 
select  to  view  the  complete  turnoff  geometry  or  may  even  again  choose  to  view  the  centerline 
turnoff  geometry  for  a  different  combination  of  turnoffs.  When  the  user  selects  to  view  the 
complete  turnoff  geometry,  the  program  requests  the  user  to  enter  the  exit  number  to  be 
displayed.  Here,  the  user  needs  to  input  only  one  exit  number  and  not  any  combination  of 
exit  numbers.  The  program  displays  the  complete  turnoff  geometry  of  that  specific  exit 
number  entered  by  the  user.  This  high  speed  turnoff  geometry  is  developed  by  the  model  for 
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the  proposed  new  turnoff(s)  and  is  shown  in  Fig.  5.16.  The  user  after  viewing  the  complete 
turnoff  geometry  of  the  turnoff  requested,  may  return  back  to  the  turnoff  location  and  their 
geometries  screen  where  (s)he  may  again  enter  a  different  turnoff  to  view  its  complete  turnoff 
geometry. 

When  the  user  selects  to  view  the  standard  turnoff  geometry  for  a  particular  turnoff,  the 
program  automatically  takes  in  the  values  of  runway  width,  taxiway  width,  the  distance 
between  the  runway  and  taxiway,  the  type  of  turnoff  and  other  data  supplied  in  the  Input 
Module. 

Aircraft  Statistics  and  their  ROT's:  This  screen  displays  the  runway  occupancy  time  for  each 
of  the  aircraft  in  the  form  of  a  bar  chart.  At  the  top  edge  of  each  of  the  bars,  the  aircraft 
number  is  displayed.  Each  aircraft  is  represented  by  two  bars.  The  first  one  represents  ROT 
for  dry  condition  whereas  the  second  one  represents  ROT  for  wet  runway  surface  condition 
For  the  benefit  of  user,  the  bars  representing  the  aircraft  and  its  name  are  of  the  same  color. 
On  the  bar  chart,  a  straight  line  is  drawn  across  the  bars  to  portray  the  weighted  average  ROT 

Each  screen  accommodates  six  aircraft,  each  with  two  bars,  one  for  wet  condition  and  the 
other  for  dry  condition.  The  user  needs  to  press  'F'  key  (for  forward)  to  view  additional  aircraft 
and  B'  key  (for  backward)  to  view  the  previous  screen.  In  each  of  the  screens  the  user  is 
presented  with  the  average  ROT,  which  encompasses  the  whole  population  of  aircraft  selected 
by  the  user.  This  includes  dry  and  wet  conditions  at  the  airport.  The  user  is  provided  with 
an  option  to  print  the  table  by  just  pressing  the  'P'  key  (for  print).  An  option  is  provided  to  go 
back  to  the  View  Menu  at  any  stage  by  entering  'V'  key  (for  View  Menu).  The  bar  chart  is 
illustrated  in  Fig.  5.15. 

Go  to  Output  Menu:  This  option  is  provided  to  enable  the  user  to  go  back  to  the  Output  Menu 
from  the  View  Menu.  The  user  may  also  exercise  this  option  by  pressing  the  escape  'Esc'  key. 
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5.2  Print  the  Report 


This  option  provides  a  hard-copy  report  which  is  a  complete  report  of  the  "Runway  Exit 
Design  Interactive  Model",  The  report  is  divided  into  two  sections.  The  first  one  deals  with  the 
summary  of  input  data  where  as  the  second  one  concentrates  on  the  results  of  the  analysis. 

The  input  data  summary  is  subdivided  into  five  categories,  as  done  earlier  in  the  Input 
Module.  Category  (1)  'Analysis  Type  and  Existing  Exits'  gives  the  type  of  analysis  selected 
and  the  number  of  exits,  their  type  and  their  location.  Category  (2)  'Aircraft  Mix  and 
Characteristics'  provides  a  table  with  the  names  of  the  aircraft  selected  and  their 
characteristics.  The  characteristics  of  the  aircraft  include:  wheelbase,  wheeltrack,  landing 
mass,  wing  area,  landing  run  distance,  load  on  main  gear,  distance  of  nose  gear  to  wingtip, 
and  the  maximum  clearance  distance.  Category  (3)  'Operational  Data'  provides  free  roll 
times,  taxiing  speed  and  their  standard  deviations,  and  the  safety  factor  for  skid.  Category  (4) 
'Environmental  Data'  provides  the  wind  speed,  wind  direction,  airport  elevation,  temperature, 
runway  visual  range,  runway  orientation,  runway  width  and  the  distance  to  the  runway. 
Category  (5)  'Runway  Gradients'  provides  the  runway  length  and  the  gradients  for  every 
one-tenth  of  the  runway  length  specified.  Category  (6)  'Weather  and  Exit  Speeds'  provides  the 
weather  conditions  (probability  of  dry  and  wet  condition)  in  percentage  and  the  speeds  for 
each  of  the  TERP  categories,  both  for  dry  and  wet  conditions.  A  sample  output  report  is  shown 
in  Figs.  6  24  and  6.25  of  this  report. 

In  the  results  of  the  analysis  section,  the  weighted  average  runway  occupancy  time  (ROT) 
and  a  table  with  the  number  of  exits,  their  location  and  their  type  is  provided.  For  runmode 
'1',  ie.,  for  the  analysis  of  an  existing  runway,  the  reliability  associated  with  each  of  the  aircraft 
and  exit,  for  both  dry  and  wet  runway  surface  conditions  is  provided  along  with  the 
corresponding  ROT's.  For  runmodes  '2'  and  '3',  ie.,  for  the  improvement  of  an  existing  runway 
or  for  designing  a  new  runway,  the  reliability  associated  with  all  the  turnoffs  is  given 
separately.  The  table  also  provides  the  turnoff  assignment  to  each  of  the  aircraft  Another 
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portion  of  the  results  pertaining  to  runmodes  '2'  and  '3'.  provides  the  coordinates  of  the 
centerline  for  each  of  the  turnoffs. 


5.3  Help 


The  main  prupose  of  this  option  is  to  help  the  user  by  explaining  each  of  the  options  that 
are  made  available  in  the  Output  Module.  The  help  screen  explains  briefly  the  first  option, 
'View  the  Output'  which  displays  the  output  on  the  screen  and  the  second  option,  'Print  the 
Report',  which  gives  the  hard-copy  of  the  complete  report.  The  user  may  press  any  key  to 
exit  from  this  screen  and  return  to  the  'Output  Menu'  screen. 

'Go  to  Main  Menu'  is  the  last  option  that  is  provided  in  the  Output  Menu.  This  allows  the 
user  to  go  back  to  the  Main  Menu  from  the  Output  Menu.  The  user  may  also  exercise  this 
option  by  pressing  the  escape  'Esc'  key. 
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6.0  Use  of  the  Model 


The  purpose  of  this  chapter  is  to  illustrate  the  use  of  the  REDIM  model  and  to  justify  the 
validity  of  the  results  obtained  throughout  the  report.  Perhaps  one  of  the  most  important 
issues  behind  this  model  is  the  capability  for  the  user  to  perform  sensitivity  analyses  for  a 
wide  range  of  airfield  environmental,  operational,  and  aircraft  dependent  variables.  As  it  was 
pointed  out  in  Section  2.1  the  model  incorporates  all  these  variables  in  a  very  interactive 
format  to  the  user  thus  minimizing  the  rerun  effort.  Fig.  6.1  illustrates  the  sensitivity  of  the 
model  to  airfield  elevation  and  exit  speed  (i.e.,  maintaining  a  constant  exit  probability).  The 
results  shown  apply  to  a  short-haul  transport  aircraft  (i.e.,  BAe  146-200)  and  depict  graphically 
the  increases  in  runway  occupancy  time  (ROT)  as  the  exit  speed  is  reduced  and  the  airfield 
elevation  is  increased.  Note  that  the  increase  in  ROT  with  decreasing  exit  speed  is  nearly 
linear  for  the  speed  ranges  tested  (10-40  m/sec.,  19.4-77.7  knots).  The  changes  due  to  airfield 
elevation  stem  from  the  larger  equivalent  airspeeds  (EAS)  during  landings  at  higher 
elevations.  The  magnitude  of  change  in  these  results  is  proportional  to  square  root  ire 
atmospheric  density  ratio.  Following  a  similar  treatment  Fig.  6.2  shows  the  variations  in  RC~ 
for  several  exit  speeds  and  airfield  temperatures  for  a  typical  short-haul  transport  aircraft  (i.e., 
Bae  146-200).  The  sensitivity  of  ROT  with  temperatures  is  again  deduced  from  the  changes 
to  the  aircraft  EAS  as  the  temperature  is  changed.  Computations  are  done  in  the  model  to 
estimate  an  equivalent  atmosphere  under  the  user-defined  conditions  and  then  estimate  the 
aircraft  equivalent  airspeed  during  the  landing  phase. 


Fig  6.3  illustrates  the  sensitivities  of  ROT  and  the  turnoff  location  parameters  with 
changing  aircraft  mass.  The  same  short-haul,  turbofan  engined  transport  aircraft  (BAe 
146-200)  operating  at  a  desired  exit  speed  of  15  m/sec.  is  used  for  illustrative  purposes.  The 
values  for  aircraft  landing  mass  cover  the  entire  allowable  landing  mass  envelope  for  this 
aircraft.  It  is  observed  that  the  variations  in  the  location  of  the  turnoff  could  be  significant  (230 
meters  between  end  points).  The  reader  should  realize,  however,  that  in  practice  a  large 
percentage  of  aircraft  are  operated  in  the  middle  of  the  region  shown  in  this  figure  (i.e., 
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80-85° b  of  the  maximum  allowable  landing  mass)  and  thus  the  changes  to  ROT  and  turnoff 
location  parameters  might  be  more  constrained  than  those  shown.  These  results  capture  the 
many  landing  performance  variations  observed  in  aircraft  flight  manuals  in  a  systematic  way. 


6.1  Example  1  (Evaluating  a  Runway  Facility) 


The  first  example  to  be  discussed  here  illustrates  the  use  of  REDIM  to  analyze  an  existing 
single-runway  airport  facility  serving  a  mix  of  general  aviation  (GA),  commuter  and  small 
transport  aircraft.  The  first  decision  faced  by  the  user  is  to  select  the  type  of  analysis  required 
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ROT  vs.  Airport  Temperature  Sensitivity 


Temperature 

♦  0  Deg.  C 
-e-  10  Deg.  C 
-o-  20  Deg  C 
30  Deg  C. 


Figure  6.2  ROT  Parameter  vs.  Exit  Speed  and  Airfield  Temperature. 


for  this  example.  Fig.  6.4  illustrates  the  main  menu  screen  where  the  'Start  a  New  Problem' 
option  is  selected  to  initiate  the  user  input  sequence.  Shortly  thereafter  REDIM  prompts  the 
user  to  name  his  working  file.  This  file  will  be  created  automatically  and  additions  will  be 
made  as  the  input  sequence  progresses.  The  model  will  ask  the  user  to  answer  yes  or  no  after 
every  input  screen  to  save  new  information.  Fig  6.5  selects  the  type  of  analysis  wanted  which 
in  this  case  corresponds  to  'Evaluation  of  an  Existing  Facility'.  Next  in  the  sequence  of  user 
inputs  is  the  definition  of  the  physical  characteristics  of  the  existing  scenario.  Lets  assume 
that  the  existing  runway  has  three  right-angled  (i.e.,  90-degree  angle  turnoffs)  located  at  both 
ends  of  the  runway  and  half  the  way  downrange.  For  a  2000  meter  long  runway  (another 
assumption  in  the  problem)  the  locations  will  be  at  0,  1000,  and  2000  meters  from  the  active 
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Landing  Mass  vs.  ROT  and  Location  Sensitivities 

BA*- 146-200  Aircraft  Data 


Maas (Tana) 


Figure  6.3  ROT  and  Exit  Location  Sensitivities  to  Aircraft  Landing  Mass  (BAe  146-200). 


threshold.  Fig.  6  6  illustrates  the  procedure  to  build  this  scenario  using  the  'Edit'  menu  screen 
containing  definitions  for  the  number  and  type  of  existing  runway  exits. 

Nine  representative  aircraft  spanning  three  different  TERP  categories,  A,  B  and  C  were 
selected  from  REDIM  aircraft  master  file  for  this  example.  Fig.  6.7  illustrates  the  aircraft  mix 
distribution  corresponding  to  an  equivalent  mix  index  of  30%.  After  the  user  has  decided  the 
complete  aircraft  population  operating  in  the  facility  pressing  the  'Esc'  key  reviews  the  aircraft 
characteristics  for  all  the  population  selected.  This  is  done  to  provide  the  user  with  some 
familiarization  of  the  aircraft  selected  and  to  allow  any  operational  changes  if  necessary.  The 
analyst  might,  for  example,  reduce  the  aircraft  landing  masses  by  a  specified  amount.  After 
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Figure  6.5  Evaluation  Mode  for  Example  1. 
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Use  arrou  <*  or  *)  ley  to  change  ijnur  choice. 
Press  enter  (  *J  )  key  If  the  choice  Is  correct. 


No.  of  existing  exits  .'  3 


<  LOCATION  (a)  > 
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30-dgr  15-djr  3£Hgr 


Figure  6.6  Selecting  Existing  Turnoffs  and  Their  Type  for  Example  1. 
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Figure  6.7  Aircraft  Mix  Editor  Screen  for  Lxample  1. 
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the  verification  of  the  aircraft  data  is  done  the  user  is  prompted  to  save  the  new  values  in  the 
user's  current  file. 

The  next  input  screen  deals  with  the  airport  operational  data  values  to  be  used  in  this 
example.  Fig.  6.8  displays  the  baseline  values  used  by  REDIM  for  the  analysis.  The  user  is  free 
to  change  any  one  of  these  values  by  typing  the  desired  values  over  the  existing  ones  A 
parameter  of  significant  importance  in  this  screen  is  the  safety  factor  for  skidding  coefficient 
to  be  used.  A  50%  safety  factor  has  proven  to  correlate  well  with  existing  empirical  turnoff 
data  (Horonjeff,  1959J  and  it  is  highly  recommended  for  a  broad  range  of  analyses  More 
conservative  designs  might  opt  for  higher  safety  factor  values  thus  increasing  the  radius  of 
curvature  of  the  proposed  optimal  turnoffs  and  also  increasing  the  resulting  weighted  average 
ROT  for  the  runway  For  this  example  the  use  of  the  baseline  values  seems  appropriate  and 
hence  Fig.  6.8  depicts  the  actual  values  used  for  this  analysis 

The  airport  environmental  conditions  used  for  this  example  are  shown  in  Fig  6  9  Salient 
features  include:  1)  calm  winds,  2)  sea  level  runway  location,  3)  25  C°  as  the  average 
temperature  of  the  hottest  month,  4)  0-18  runway  orientation  and  5)  a  distance  of  280  m  (918 
ft.)  from  runway  to  taxiway  centerlines.  This  latter  parameter  will  be  used  in  the  output 
module  to  construct  a  complete  high-speed  geometry  to  the  nearest  taxiway  using  a  terminal 
exit  angle  of  30  degrees.  However,  this  topic  is  currently  being  studied  to  investigate  various 
extended  turnoff  configurations  that  will  be  used  selectively  according  to  several 
runway-taxiway  configurations.  For  more  information  regarding  the  sensitivity  of  the  'terminal' 
turnoff  angle  us*d  refer  to  Section  6.3  of  this  document.  The  penultimate  input  screen  in  the 
'Evaluation  Mode'  prompts  the  user  to  select  the  local  runway  gradients  in  tenth's  of  the  total 
runway  length.  This  is  shown  in  Fig  6.10.  In  this  scenario,  a  2000  meter  runway  is  readily 
divided  into  ten  200  meter  segments  to  which  a  local  gradient  value  is  associated  For  the 
purpose  of  this  example  we  use  a  constant  -.5  %  gradient  (downslope)  throughout  the 
complete  runway  length.  Note  that  the  notation  used  in  the  program  is  consistent  with  that 
familiar  to  airport  and  highway  engineers.  Also,  REDIM  has  a  a  built-in  check  routine  to  verify 
that  local  gradients  will  not  exceed  the  maximum  allowable  by  FAA  standards.  This 
verification  is  accomplished  prior  to  the  actual  simulation  and  optimization  procedures  But 
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Figure  6.8  Airport  Operational  Values  for  Example  1. 
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Figure  6.9  Airport  Environmental  Values  for  Example  1. 
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Figure  6.10  Runway  Gradients  for  Example  1. 
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the  user  can  correct  the  error  Interactively  by  means  of  the  Main  Menu'  which  is  always 
within  reach  through  the  'Esc'  key. 

Finally,  in  the  Weather  and  Exit  Speeds'  screen  (Fig.  6.11)  the  user  enters  the  relative 
frequency  of  occurence  of  weather  conditions  to  be  factored  in  the  analysis  Two  weather 
conditions  are  modeled  in  the  current  program,  wet  and  dry  runways,  allowing  more  flexibility 
from  an  operational  point  of  view.  In  general  terms,  wet  scenarios  will  result  in  larger  values 
of  weighted  average  ROT  times.  However,  this  might  be  deemed  necessary  by  the  user  in 
order  to  account  for  airport  specific  conditions  at  the  location  being  analyzed.  Overall,  the 
runway  turnoff  designs  will  also  be  more  conservative  with  larger  radii  of  curvature  and 
further  downrange  turnoff  locations.  In  this  case  equal  weights,  50%  probabilities,  are  given 
to  both  runway  conditions. 

This  concludes  the  input  set  for  this  first  example.  At  this  point,  the  analyst  is  expected 
to  return  to  the  Main  Menu'  through  the  'Esc'  key  from  where  the  model  analysis  routines 
(i.e.,  simulation  and  optimization)  are  invoked  selecting  'Begin  Analysis'  from  this  menu  Th;s 
starts  the  landing  simulation  of  every  aircraft  subjected  to  the  operational  parameters  input 
by  the  user  to  find  candidate  exit  locations,  their  geometries,  and  finally  to  select  those 
considered  optimal  according  to  a  minimum  weighted  ROT  performance  index  criterion.  The 
execution  of  the  simulations  and  optimization  routines  can  take  anywhere  from  10  seconds  to 
a  few  minutes  depending  upon  the  number  of  aircraft  selected  In  order  to  provide  some  visual 
feedback  to  the  user  through  the  simulation  process  the  user  is  exposed  to  relevant  statistics 
for  each  aircraft  simulated.  The  statistics  include;  1)  the  aircraft  type  designator.  2)  the 
individual  runway  occupancy  time,  and  the  3)  exit  location  used.  For  this  nine  aircraft  example 
the  computation  time  is  about  10  seconds  for  the  dynamic  simulation  In  this  case  no 
optimization  is  necessary  as  only  the  'Evaluation  Mode'  subroutines  are  invoked. 

The  results  of  this  single-runway  scenario  are  shown  interactively  in  Figs.  6.12  to  6.15 
which  are  part  of  the  Output  Module  routines.  Fig.  6.12  illustrates  the  main  'Output  Menu' 
screen  shown  to  the  user  where  'View  the  Output'  and  'Print  the  Report'  constitute  the  two 
alternatives  to  obtain  screen  and  printed  output,  respectively.  Selecting  the  first  option,  'View 
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the  Output'  the  analyst  is  exposed  to  Fig.  6.13  where  a  triple  selection  screen  directs  the  user 
to  the  runway  occupancy  table  (ROT  Table),  the  turnoff  locations  and  their  geometries  or  back 
to  the  output  menu.  The  ROT  table  depicted  in  Fig.  6.14  shows  a  partial  view  of  the  ROT  Table 
where  individual  ROT  times  and  the  complete  aircraft  assignment  is  made  for  to  the  three 
existing  exits.  Notice  that  exit  number,  location  and  exit  type  are  clearly  identified  in  the  first 
row  of  the  table.  Thereafter,  every  aircraft  is  identified  by  its  program  designator  (i.e.,  see 
Appendix  A  for  more  details  on  designators  and  aircraft  representation).  As  detailed  in 
Chapter  5  of  this  document,  two  scenario  conditions  are  analyzed  by  REDIM,  wet  and  dry 
runway  conditions  with  relative  frequencies  of  occurence  specified  by  the  user  (see  Fig.  6.11). 
The  interpretation  of  the  ROT  table  results  is  as  follows:  every  aircraft  is  assigned  to  one  or 
several  turnoffs  where  potentially  a  successful  exit  maneuver  can  be  executed.  Taken  as 
example  the  Swedish  made  commuter  aircraft  SAAB  340  the  runway  occupancy  time  for  the 
dry  scenario  is  44.2  seconds  taking  the  second  turnoff  (Exit  #  2),  located  1000  meters  from  the 
active  runway  threshold.  The  probability  for  this  aircraft  taking  this  middle  exit  is  only  7.9  % 
suggesting  average  landing  rolls  greater  than  1000  meters.  Notice  that  if  the  third  exit  is  used 
to  clear  the  runway  the  remaining  percent  of  the  population,  92.1%,  is  able  to  exit  at  the 
expense  of  a  large  ROT  value  (138.6  sec.)  as  the  aircraft  is  required  to  travel  at  near  \a/-  'z 
speed  for  the  remaining  portion  of  runway  This  value  should  be  viewed  only  as  an  upper  limit 
since,  ground  operations  permitting,  the  SAAB-340  will  probably  execute  a  180-degree  turn 
and  still  take  the  second  exit.  Under  this  new  set  of  conjectures  an  estimated  ROT  time  closer 
to  95  seconds  is  more  realistic.  Following  the  same  aircraft  it  is  noted  that  wet  conditions 
lower  the  ROT  time  through  the  third  exit  since  the  aircraft  requires  longer  braking  distances 
with  the  corresponding  reduction  in  the  ROT  devoted  to  taxi  to  the  next  turnoff.  Note,  however, 
that  the  percent  of  the  SAAB-340  population  taking  the  second  exit  is  only  8%  implying  that 
a  small  percentage  of  the  operations  will  be  able  to  use  this  exit. 

Table  6.1  provides  the  baseline  exit  speed  values  used  in  REDIM  to  predict  the  turnoff 
location  and  reliability  parameter  for  individual  aircraft  using  three  FAA  standard  turnoff 
geometries.  Note  that  for  the  90-Degree  angle  turnoff  the  entry  speed  is  defined  by  the  user 
as  a  taxiing  speed.  This  speed  represents  the  safe  value  at  which  a  pilot  will  comfortably 
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maneuver  his  airplane  on  the  ground  to  reach  the  nearest  turnoff  location  once  a 
predetermined  exit  threshold  value  has  been  reached. 


Table  6.1 

REDIM  Baseline  Exit  Speeds  for  Standard  FAA  Turnoffs. 

Turnoff  Type 

Turnoff  Entry  Speed 

90-Degree 

Taxiing  Speed  (User  Defined) 

Typically  8  m./sec.  (18  MPH) 

45-Degree 

17.9  m./sec.  (40  MPH) 

30-Degree 

26.9  m./sec.  (60  MPH) 

In  Fig.  6.14  the  weighted  average  runway  occupancy  time  (WAROT)  is  also  indicated  for 
the  complete  population  analyzed.  In  this  case  102.94  seconds  represents  a  large  WAROT 
value  for  this  simple  example.  According  to  this  result  a  maximum  of  34  landings  per  hour 
would  be  the  upper  limit  for  this  single  runway  under  the  given  conditions.  For  mixed 
operations  this  value  could  increase  by  another  10%  or  so.  The  question  is  how  much  can  this 
facility  be  improved  by  adding  more  turnoffs?  The  answer  to  this  is  the  subject  of  Example  2 
to  be  discussed  in  the  next  section. 


6.2  Example  2  (Improving  a  Runway  Facility) 


This  problem  is  an  extension  to  the  previous  one  as  it  was  observed  that  the  existing 
turnoff  locations  were  'inefficient'  to  handle  the  hypothesized  aircraft  population  resulting  in 
large  ROT  times.  The  idea  behind  this  second  example  is  to  improve  the  existing  single 
runway  design  through  the  incorporation  of  additional  high-speed  turnoffs.  The  location  and 
geometry  of  these  will  be  found  by  REDiM's  dynamic-optimization  algorithms.  Since  it  is 
pressumed  that  the  analyst  has  created  a  file  with  the  airport  specifications  in  the  previous 
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Figure  6.16  Existing  90-Degree  Standard  Turnoff  Geometry  for  Example  1. 


example  and  only  those  screens  having  significant  changes  from  those  of  Example  1  will  be 
treated  in  detail. 

There  are  two  ways  to  approach  this  revised  problem.  The  first  one  is  to  use  the  existing 
data  file  from  Example  1  and  make  the  proper  modifications  through  the  complete  input 
sequence.  The  second  approach  is  to  start  a  new  problem  with  the  same  parameters  as 
Example  1  and  complete  the  pieces  of  information  left  out  in  the  previous  problem.  Taking  the 
second  approach  is  simpler  for  the  novice  but  the  first  one  saves  time. 

The  first  significant  change  is  made  to  the  'Analysis  Type'  menu  (see  Fig.  6.5)  where  now 
the  second  choice  is  selected,  'Improve  an  Existing  Runway'.  From  this  point  on  the  same 
screens  as  those  associated  with  Example  1  will  be  applicable.  A  reliability  parameter  needs 
to  be  specified  in  the  third  'Input'  screen  to  estimate  the  degree  of  exit  reliability  expected 
from  the  new  turnoffs.  For  this  example  we  use  90%  as  reliability  parameter  and  we  will 
specify  the  number  of  exits  to  be  built  at  the  end  of  the  simulation  stage.  Another  aspect  that 


Use  of  the  Model 


114 


we  need  to  specify  in  this  new  scenario  concerns  the  desired  entry  exit  speeds  for  the  new 
turnoffs.  This  is  crucial  as  part  of  the  input  as  REDIM  will  use  the  entry  turnoff  speeds  selected 
in  the  'Weather  and  Exit  Speeds'  menu  to  design  geometrically  the  optimal  turnoffs.  For  the 
improved  runway  the  exit  speeds  and  weather  conditions  are  shown  in  Fig  6.17.  Note  that  it 
is  possible  to  leave  voids  of  speed  data  as  far  as  these  are  not  used  in  the  program.  For 
example,  its  is  allowed  to  leave  blanks  for  the  speed  values  associated  with  TERP  categories 
D  and  E  as  there  are  no  aircraft  belonging  to  these  categories  in  this  example. 

After  the  pertinent  changes  have  been  made  to  the  Input  Module  screens  the  user  is 
expected  to  go  back  to  the  'Main  Menu'  screen  {pressing  the  'Esc'  key)  and  start  the  analytical 
procedures  of  REDIM.  Selection  of  item  3  in  Fig.  6.5  begins  the  dynamic  simulation  analysis. 
Just  as  for  the  previous  example,  there  is  some  feedback  information  displayed  on  the  screen 
in  terms  of  partial  ROT  values  for  every  candidate  solution  generated  (see  Fig.  6.18).  The 
dynamic  computations  take  on  the  average  9  seconds  per  aircraft-scenario  combination. 
Currently,  the  model  is  restricted  to  25  aircraft  per  run  (i.e.,  50  aircraft-scenario  combinations) 
due  to  internal  array  size  limitations.  This,  however,  seems  to  be  sufficient  for  most  of  the 
airport  scenarios  to  be  encountered  and  should  not  restrict  the  validity  of  the  results.  Once  the 
dynamic  simulation  is  done,  the  analyst  enters  the  number  of  new  turnoffs  to  be  constructed 
and  the  optimization  module  routines  are  executed.  Fig.  6.19  illustrates  a  partial  view  of  the 
optimization  results  using  5  exits.  That  is,  constructing  two  new  high-speed  turnoffs  to 
complement  the  three  existing  ones. 

The  new  average  ROT  value  being  45.91  seconds  represents  a  significant  improvement 
over  the  previous  example.  The  suggested  new  locations  are  727  and  1495  meters  from  the 
active  threshold.  Fig.  6.19  also  details  the  percentages  of  each  aircraft-scenario  combination 
exiting  through  each  turnoff  and  their  corresponding  individual  ROT  times.  It  should  be  noticed 
that  REDIM  constraints  adjacent  turnoff  locations  within  a  prescribed  distance  Dm(n  to  abide 
current  FAA  standards.  The  current  value  for  Dmm  is  considered  to  be  213  meters  (i.e.,  700  ft.) 
but  this  can  be  modified  by  the  user  in  the  'Airport  Environmental  Screen'  (see  Fig.  6.9). 
Running  REDIM  under  the  same  runway  environmental  conditions  but  increasing  the  number 
of  new  turnoffs  to  three  instead  of  two  as  before  a  15%  improvement  is  observed  in  the 
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Figure  6.18  Partial  Dynamic  Simulation  Results  for  Improved  Single  Runway. 
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average  ROT  time  (from  45.91  to  38.57  seconds).  Fig.  6.20  shows  a  partial  ROT  table  for  this 
new  scenario.  Note  that  the  new  turnoff  location  are:  727, 1279  and  1495  meters  from  the  active 
threshold.  As  one  might  suspect  an  increase  in  the  number  of  high-speed  turnoffs  yields  better 
runway  service  times  (i.e..  lower  values  of  WAROT)  at  the  expense  of  capital  cost.  It  can  also 
be  shown  with  several  consecutive  runs  of  the  model  that  the  gains  in  WAROT  are  small  for 
a  large  number  of  exits  (n  >  6). 

Figures  6.21-6.23  depict  graphically  the  location  and  geometries  generated  by  REDIM  in 
the  5-turnoff  runway  scenario.  Fig.  6.22  shows  a  characteristic  compound  plot  of  the  five  turnoff 
geometries,  three  already  available  (i.e.,  standard  FAA  90-Deg.  turnoffs)  and  two  more 
projected.  Fig.  6.23  represents  the  complete  turnoff  geometry  of  the  fourth  turnoff  located  1495 
m.  from  the  runway  threshold.  In  the  'Print  the  Report'  Option  the  user  receives  a  complete 
report  on  the  optimization  results  as  well  as  the  input  parameters  selected  for  that  particular 
run.  The  report  is  divided  into  two  sections:  I)  input  data  analysis  and  II)  analysis  results.  The 
former  is  in  term  subdivided  into  six  categories  corresponding  to  each  one  of  the  program 
input  screens:  1)  type  of  analysis,  2)  aircraft  mix,  3)  airport  operational  data  4 1  airport 
environmental  data,  5)  runway  gradients,  and  6)  weather  and  speed  characteristics  The 
analysis  results  section  of  the  report  contains  three  sub-sections;  1)  average  ROT,  2)  exit 
locations,  type  and  turnoff  assignment  table,  and  the  actual  centerline  turnoff  coordinates. 
Figs.  6.24-6.25  show  partial  listings  of  the  report  generated  for  Example  2  (5-turnoff  case) 
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1*1.  Analysis  type  and  Existing  Exits 

Analysis  type  ■  Improvement  of  an  Existing  Runway 


Exit  # 

Location  <■) 

Type 

1 

0 

90-dcg 

2 

1000 

90- deg 

3 

2000 

90- deg 

Figure  6.24  Input  Data  Summary  Report  (or  Example  2. 


Use  of  the  Model 


1-2.  Aircraft  Mix  ond  Characteristics 


acf  name  l 

r 

X 

UB 

(m) 

UT 

(m) 

LM 

(X) 

MASS 

(kg) 

L0 

(m) 

CL  MAX 

UA 

(m‘2) 

US 

(m) 

Nur 

(m) 

PA-38- 112  1 

10.0 

1.5 

3.0 

77.4 

757.0 

486.0 

1.575 

11.6 

10.4 

2.1 

PA-78-161  l 

10.0 

2.0 

3.0 

82.3 

1109.0 

416.0 

1.694 

15.8 

10.7 

2.4 

BE  -  58  1 

10,0 

2.7 

2.9 

84.6 

2500.0 

751.5 

1.486 

18.5 

11.5 

2.9 

Bl -  300  1 

10.0 

4.6 

5.2 

89.0 

6363.0 

857.2 

2.076 

28.2 

16.6 

5.3 

CE-407C  1 

10.0 

3.2 

5.5 

88.1 

3107.0 

655.8 

2.100 

21.0 

13.4 

3.8 

SAAB- 340  1 

10.0 

7.1 

6.7 

90.9 

12020.0 

1  J40.4 

2.578 

41.8 

21.4 

7.7 

FM8-120  1 

10.0 

7.0 

6.6 

90.5 

11250.0 

1269.5 

2.2/2 

39.4 

19.8 

7.6 

fOKKFR-1001 

15.0 

14.0 

5.0 

89.5 

39915.0 

1360.0 

2.531 

93.5 

28.1 

16.8 

BAo  146  I 

15.0 

11.2 

4.7 

92.3 

36  740.0 

1130.0 

3.305 

77. 3 

26.3 

12.6 

uo  -  Wheelbase 
MASS  *  landing  Moss 

UA  -  Uinq  Area 


WT  =  wheel  track  IM  *  lood  on  Mam  Gear 

LD  *  Landing  Run  Distance  CL  MAX  - 

US  *  Wing  Span  NUI  s  Dist.  Hose  Gear  to  Uingtip 


I -3.  Operational  Oata 


1st  free  rol l  t ime  : 

2.0 

(see) 

std.  dev. 

0.5 

2nd  free  rot l  time 

1.0 

(see) 

std.  dev. 

0.2 

taxing  speed  : 

safety  fac .  for  sk id  : 

8.0 

50.0 

(m/s) 

(X) 

std.  dev. 

1.0 

.  £nviror«ncntal  Oata 

wind  speed  : 

0.0  (m/s) 

wind  direction 

0.0 

airport  elevation 

0.0  (m) 

temperature 

25.0  (C-deg) 

min.  exit  interval  : 

213.0  (m) 

runway  orientation 

0.0 

runway  width 

45.0  On) 

distance  to  taxiwoy  : 

280.0  <m) 

1*5.  Runwny  Gradients 

runway  length  :  2000  (m) 

gradients  (X) 

-0.5  0.5  0.5  -0.5  0.5  -0.5  -0.5  -0.5  -0.5  -0.5 

I*  - - ♦ . ♦ . ♦ . ♦ . ♦ . * .  ♦ . ♦ .  I 

1000  7000 


1-6.  Weather  ond  Exit  Speeds 


weather  &  exit  speed  (m/s) 


1 

DRY 

1 

WET 

Probability  (X)  1 

50.0 

1 

50.0 

TERPS 

A  1 

20.0 

1 

20.0 

TERPS 

B  1 

25.0 

1 

25.0 

TERPS 

C  1 

30.0 

1 

30.0 

TERPS 

0  I 

30.0 

1 

30.0 

TERPS 

E  1 

0.0 

1 

0.0 

Figure  6.24  Input  Data  Summary  Report  for  Example  2  (Continuation), 


II.  ANALYSIS  RESULTS 


11-1.  Average  ROT 

average  ROT  •  AS. 91  (sec) 


-2.  Exit  Loca 

tions.  Types,  and  Turn-off  Assignment. 

Exit  t 

1  1 

2 

3 

A 

5 

Location  (m) 

1  0 

m 

1000 

1A95 

2000 

Type 

1  90-d 

new 

90-d 

new 

90-d 

PA-38- 1 12 
dry 

wet 

35.0 
34.  A 

PA-28- 161 
dry 

wet 

36.1 

35.7 

BE-58 

dry 

wet 

71. A 

68.1 

BE -300 

dry 

wet 

6A.0 

60.9 

CE-A02C 

dry 

wet 

25.0 

ZA.5 

SAAB-3A0 

dry 

wet 

55.5 

52.3 

EH8- 120 

dry 

wet 

50.1 

A7.0 

fOKKER-100 

dry 

wet 

41.1 

39.0 

BAC-1A6 

dry 

wet 

A6.9 

45. 1 

Figure  6.25  Analysis  Results  Format  for  Example  2. 


Um  of  the  Model 


1-3.  Turn-off  Centerline  Geometries 


1  Exit 

1 

1  Exit 

2  I 

Exit 

3  I 

Exit 

4  1 

1  X  1 

Y 

1  X  1 

Y  1 

X  I 

Y  1 

X  I 

Y  I 

1  5.1 

0.2 

I  19.8 

0.5  I 

5.1 

0.2  I 

29.9 

0.2  1 

1  10.1 

0.7 

1  39.3 

2.6  1 

10.1 

0.7  1 

59.4 

1.2  1 

1  15.2 

1.5 

1  58.1 

6.7  1 

15.2 

1.5  1 

88.6 

3.3  1 

1  20.3 

2.8 

1  76.1 

12.7  I 

20.3 

2.8  I 

117.4 

6.5  1 

1  25.4 

4.3 

I  93.1 

20.5  1 

25.4 

4.3  1 

145.7 

11.0  1 

1  30.4 

6.3 

1  109.0 

29.6  1 

30.4 

6.3  1 

173.5 

16.8  1 

1  35.5 

8.8 

I  124.7 

38.6  i 

35.5 

8.8  I 

200.7 

23.8  1 

I  40.6 

11.7 

1  140.1 

47.5  I 

40.6 

11.7  1 

227.2 

32.1  1 

1  45.6 

15.2 

1  155.3 

56.3  I 

45.6 

15.2  1 

253.0 

41.7  1 

1  50.7 

19.4 

1  170.1 

64.9  1 

50.7 

19.4  1 

278.0 

52.4  1 

55.8 

24.3 

1  184.8 

73.4  1 

55.8 

24.3  1 

302.1 

64.3  I 

I  60.8 

30.4 

1  199.2 

81.7  1 

60.8 

30.4  I 

325.4 

77.2  1 

1  65.9 

38.1 

1  213.3 

89.8  1 

65.9 

38.1  1 

348.2 

90.4  1 

1  71.0 

48.7 

I  227.2 

97.9  1 

71.0 

48.7  | 

370.7 

103.4  I 

I  76.1 

76.1 

1  240.8 

105.7  1 

76.1 

76.1  I 

393.0 

116.3  I 

1  76.1 

140.0 

1  254.2 

113.5  1 

76.1 

140.0  I 

415  0 

129.0  1 

1 

1  267.3 

121.0  1 

1 

434.2 

140.1  1 

] 

I  280.2 

128.5  1 

I 

1 

I 

1  292.8 

135.7  I 

1 

1 

1 

1  300.2 

140.0  1 

I 

1 

I  Exit 

5  1 

1  X  1 

Y  1 

1  5.1 

0.2  1 

I  10.1 

0.7  1 

1  15.2 

1.5  1 

I  20.3 

2.8  1 

1  25.4 

4.3  I 

1  30.4 

6.3  1 

1  35.5 

8.8  1 

1  40.6 

11.7  1 

1  45.6 

15.2  1 

1  50.7 

19.4  1 

I  55.8 

24.3  1 

1  60.8 

30.4  1 

I  65.9 

38.1  1 

1  71.0 

48.7  1 

1  76.1 

76.1  1 

1  76.1 

140.0  1 

Figure  6.25  Analysis  Results  Format  for  Example  2  (Continuation). 


7.0  Conclusions 


The  end  result  of  this  model  is  to  recommend  a  high-speed  geometry  that  will  minimize 
the  runway  occupancy  time  under  realistic  airport  scenarios.  As  it  was  explained  in  Chapter 
6  the  model  is  able  to  predict  turnoff  locations  and  geometries  that  optimize  the  weighted 
average  ROT  parameter  for  a  given  set  of  airport  conditions.  The  obvious  question  is  how 
sensitive  are  the  model  results  in  terms  of  some  of  the  input  parameters  such  as  aircraft  mix, 
wind  conditions,  airfield  parameters  and  so  on.  This  question  arises  naturally  since  these 
variables  are  highly  dynamic  and  fluctuate  during  the  day  and  from  season  to  season. 

Looking  at  existing  data  on  ru. iway  occupancy  time  [Koenig,  1978;  Ruhl,  1989)  it  is 
believed  that  REDIM  is  behaving  in  a  realistic  fashion  for  a  multitude  of  scenarios  tested.  Fig. 
7.1  illustrates  the  results  for  San  Francisco  International  Airport  runway  27R  where  two 
independent  sets  of  data  were  compared  with  the  predictions  made  by  REDIM  for  the  same 
scenario.  Fig.  7.2  depicts  the  results  of  REDIM  and  the  observations  reported  by  Koenig 
(Koenig,  1978).  In  all  cases  the  differences  are  below  5%  from  each  other.  A  generalized  trend 
on  ROT  times  versus  mix  index  for  various  numbers  of  exits  is  illustrated  in  Fig.  7.3.  In  this 
figure  it  is  seen  the  sensitivity  of  the  model  to  the  number  of  exits.  Without  any  doubt  one  of 
the  most  important  parameters  influencing  ROT  times. 


7.1  Suggested  High-speed  Standard  Geometry 


The  implementation  of  realistic  high-speed  turnoffs  seems  to  be  one  the  most  debatable 
issues  faced  by  airport  engineers.  On  one  hand  it  is  well  known  that  the  location  of  the  runway 
turnoffs  affects  significantly  geometry  of  every  turnoff;  However,  for  a  finite  aircraft  population 
a  single  turnoff  location  scheme  is  needed  to  minimize  the  desired  average  ROT  performance 
index.  The  problem  seems  then  to  be  that  each  planner  should  use  variable  geometry  turnoffs 
for  every  scenario.  This  makes  the  number  of  geometry  choices  almost  limitless  for  the 
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Figure  7.1  Comparison  of  REDIM  Results  at  San  Francisco  Inti.  Runway  28R. 


hundredths  of  airport  scenarios  around  the  nation.  This  issue  is  not  time  consuming  with  the 
use  of  REDIM,  but  even  a  complex  model  like  the  one  addressed  in  this  report  makes  use  of 
generalizations  in  order  to  reduce  the  magnitude  of  the  problem  within  a  confined  set  of 
choices. 


In  dealing  with  a  new  standard  geometry  our  approach  to  the  standardization  problem  is 
very  similar  to  that  confronted  by  the  Horonjeff  team  three  decades  ago.  If  a  standard  is  to  be 
accepted  by  the  aviation  community  it  not  only  needs  to  be  proven  in  simulators  and  in  fields 
demonstrations,  but  also  needs  to  address  the  needs  of  the  builder  in  terms  of  a  simple 
definition  of  the  geometry.  This  is  probably  the  most  difficult  task  to  address  since  a  fully 
variable  geometry  is  obtained  as  the  result  of  the  turning  equations  of  motion  of  aircraft 
negotiating  a  high-speed  turnoff.  The  specification  of  such  a  geometry  (i.e.,  fully  variable 
geometry)  is  difficult  to  justify  in  practice  since  every  position  coordinate  in  a  two-dimensional 
plane  needs  to  be  known.  From  an  operational  point  of  view  it  is  possible  to  approximate 
slow-varying  turnoff  geometries  (i.e.,  spirals  and  clotoids]  with  large  radius  of  curvature 
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Figure  7.2  Comparison  of  REDIM  Results  at  Denver  Inti.  Rurway  26R. 


entrance  curve  followed  by  a  reduced  radius  of  curvature  circular  segment.  This  approach 
was  suggested  by  Robert  Horonjeff  in  the  late  50's  (Horonjeff,  et  al.  1959)  but  interestingly 
enough  his  results  have  not  been  universally  accepted  by  all  (he  aviation  authorities  in  terms 
of  adopting  a  large  entrance  curve  as  geometric  design  standard.  The  current  FAA  practice 
uses  a  single  radius  of  curvature  to  define  the  geometry  of  a  high-speed  turnoff  (i.e.,  1800  ft 
for  30-Deg.  angled  exits).  In  our  findings  with  REDIM  we  have  to  acknowledge  that  HoronjeTs 
suggestions  were  justifiable  and  that  possibly  the  simplest  approach  to  define  a  new  standard 
is  to  consider  two  circular  arcs  with  a  common  tangency  point  as  a  viable  solution  to 
approximate  a  fully  variable  turnoff  geometry  (see  Fig.  7.4).  This  approach  is  revisited  in  this 
section  to  show  the  selection  process  behind  the  variable  geometry  standard. 


From  Fig.  7.4  it  is  seen  that  two  radii  of  curvature  defined  R ,,  Rz,  and  a  turnoff  exit  angle, 
^  form  the  basis  for  the  suggested  approximation.  The  first  radius  of  curvature  approximates 
the  jerk-limited  curve  corresponding  to  a  specified  entry  speed  ( V ,)  whereas  the  second  one, 
Rt,  models  the  aircraft  "steady-state  rotational"  inertia  characteristics  as  it  negotiates  the 
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Figure  7.4  Definition  of  the  High-Speed  Turnoff  Geometry. 


turnoff.  Through  many  simulations  using  REDIM  it  became  evident  that  extracting  two  specific 
values  of  R  an  excellent  approximation  to  this  fully  variable  turnoff  geometry  could  be 
obtained.  The  values  of  /?,  and  R2  then  were  obtained  as  a  function  of  turnoff  time  and  aircraft 
category. 

The  rationale  behind  the  time  factor  in  this  recommendation  is  to  account  for  the  aircraft 
inertia  resistance  motion  which  can  be  categorized  as  a  ’pseudo-first  order  model'  (see  Eqns. 
3.18-19  for  R)  where  the  radius  of  curvature  changes  slowly  as  a  function  of  time.  Looking  at 
Fig.  7.5  it  is  observed  that  an  equivalent  "time  constant’  characterizing  the  aircraft  rotational 
motion  about  the  z  axis  as  it  negotiates  a  high-speed  turnoff  is  proportional  to  the  aircraft 
mass  and  moment  of  inertia  about  this  axis  among  other  factors.  Knowing  this  fact  a  straight 
correlation  between  the  values  of  Rt  and  an  extraction  time  were  established.  Table  7.1 
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summarizes  the  nominal  extraction  times  used  in  REDIM  to  approximate  the  variable  turnoft 
trajectory. 

In  REDIM  nomenclature  these  times  are  labeled  as  easement  curvature  time.  TR1,  and 
steady-state  curvature,  time,  TF.2.  Note  that  for  heavy  transport-type  aircraft  (i.e.,  >  300,000 
lbs)  larger  time  lags  to  achieve  a  'steady-state'  radius  of  curvature  are  a  direct  result  of  larger 
time  constants  in  the  model. 

Fig.  7.4  also  illustrates  the  two  corresponding  encompassing  the  approximate  turnoff 
track.  Arcs  with  lengths  Li  and  Lt  are  defined  as  follows, 


L1  =  R,0,  {7.1} 

L2  =  R2  02  {7.2} 

where,  Li  and  L*  represent  the  turnoff  characteristic  lengths.  R,  and  R2  are  the  radii  of 
curvature  defining  the  turnoff,  and  0 ,  and  0t  are  the  arcs  defined  by  Ri  and  R2,  respectively 
measured  in  radians.  The  turnoff  arcs  are  characteristic  for  each  aircraft  since  the  transition 
and  runway  clearance  point  are  aircraft  speed  and  geometry  dependent.  It  should  be  kept  in 
mind  that  Li  is  a  linear  function  of  aircraft  speed  if  the  jerk-limited  equation  is  used  and  if  the 
values  of  a„  and  J„  are  substituted  in  Eqn.  7.3. 


,  amax 

Lf  —j -  -  —j -  {7.3} 

The  analyst,  however,  does  not  need  to  be  concerned  in  REDIM  since  the  actual  turnoff  track 
values  are  presented  in  tabular  form.  The  approximation  is  primarily  used  to  depict  the 
geometry  on  the  computer  screen. 
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Table  7.1  REDIM  Model  Extraction  Times  for  Estimating  Ri  and  R7. 

Category 

*1 

Ri 

A 

1  sec. 

4  sec. 

B 

1  sec. 

5  sec. 

C 

1  sec. 

6  sec. 

D 

1  sec. 

6  sec. 

E 

1  sec. 

6  sec. 

7.2  Comparison  of  REDIM  Geometries 


The  geometries  generated  by  REDIM  are  dictated  primarily  by  the  jerk  and  normal 
acceleration  in  the  first  few  seconds  of  the  trajectory  and  by  the  aircraft  rotational  inertia 
limitations  in  the  longer  term  <i.e.,  3  or  more  seconds  into  the  turn).  In  general,  the  geometries 
obtained  in  REDIM  differ  from  the  FAA  standard  acute  angle  exit  geometry  in  terms  of  their 
initial  and  steady  state  radii  of  curvature.  Fig.  7.6  depicts  two  exit  geometries  corresponding 
to  and  exit  speed  of  27  m.sec.  (60  MPH).  The  top  geometry  corresponds  to  the  standard  acute 
angle  exit  and  is  shown  for  the  sake  of  comparison.  The  bottom  geometry  was  generated  by 
REDIM  for  a  Boeing  727-200  operating  on  a  wet  runway.  Note  that  in  both  examples  the  finai 
exit  angle  has  been  maintained  at  30  degrees  and  as  can  be  seen  the  REDIM  geometry  is 
characterized  by  two  radii  of  curvature  (R1  =979.6  m.  and  R2  =  447.8  m.)  resulting  in  a  slightly 
larger  arc  length  to  reach  the  final  exit  angle.  Also  shown  in  Fig.  7.6  is  a  superposition  of  both 
geometries  revealing  in  greater  detail  their  differences.  Notice  that  the  width  of  the  turnoff  has 
also  been  maintained  at  30.5  m.  (i.e.,  100  ft.)  for  the  purpose  of  illustration.  It  is  important  to 
realize  that  currently  REDIM  evaluates  the  centerline  of  the  aircraft  trajectory  and  it  designs 
the  turnoff  edges  according  to  the  aircraft  design  group  classification.  It  seems  advisable. 
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Figure  7.5  Time  Variations  of  Radius  of  Curvature  for  Representative  Aircraft. 


however,  to  increase  the  existing  turnoff  widths  in  order  to  increase  the  pilot's  confidence 
while  negotiating  a  turn  at  high  speed. 

Fig.  7.7  illustrates  similar  results  for  a  Boeing  747-200.  In  this  particular  case  the  differences 
are  more  notorious  as  the  larger  aircraft  mass  and  inertia  penalize  its  maneuverability  on  the 
ground.  Again,  the  lateral  distance  and  the  turnoff  width  has  been  maintained  according  to 
existing  FAA  standards.  Fig.  7.8  illustrates  the  geometry  for  the  same  Boeing  747-200  when  the 
entry  speed  is  35  m./sec.  {78  MPH).  In  this  case  the  turnoff  width  has  been  increased  to  45.8 
mts.  in  order  provide  better  situational  awareness  to  the  pilot  It  is  believed  that  increases  in 
the  width  of  all  high  speed  exits  will,  in  general,  induce  pilots  to  maintain  faster  exit  speeds 
than  those  seen  today  at  major  airports.  When  one  considers  night  and  wet  pavement 
conditions  and  factors  a  reasonable  skidding  friction  parameter  the  resulting  geometries 
require  significant  longitudinal  and  lateral  distances  to  allow  sizeable  speed  reductions  on  the 
turnoff.  It  has  been  estimated  that  230  mts.  (750  ft.)  seems  to  be  the  minimum  lateral  distance 


Conclusion* 


131 


Conclusion* 


132 


allowing  safe  negotiation  of  a  35  m./sec.  entry  speed  turn.  This  distance  depends  upon  the 
terminal  speed  required  at  the  taxiway  junction  point  and  will  be  the  subject  of  further  analysis 
in  the  Phase  II  of  this  research. 


Conclusions 


133 


Run  i 


i»v 


r AA  Standard  Artitp  Anqlefxlt 


Rnnwny 


fttniri  ftpneratpd  High-Speed  Fxit 
(Hoeing  /4/  ?00  Parameters) 


Runway 


f  AA  Siftnrfnrd  Acute  Angle  t  *  i» 

RFDIM  Of* 1 1 or  alert  HS  F*rt  (  B  74?  ?00  Pirwmlw  *») 


Figure  7.7  Heavy  Transport  Turnoff  Geometries  Comparison  (27  m./sec.) 
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8.0  Model  Recommendations 


Although  the  REDIM  model  addresses  a  large  variety  of  parameters  affecting  the  aircraft 
landing  dynamics  and  the  airport  environmental  characteristics  several  features  of  the  model 
need  further  investigation  in  order  to  calibrate  and  verify  some  of  the  assumptions  made 
during  the  model  development.  Among  these  features  are:  1)  a  postoptimization  technique  to 
account  for  time  varying  aircraft  mixes  and  airport  environmental  parameters,  2)  human 
behavioral  factors  such  as  the  verification  of  the  lateral  acceleration  and  jerk  perception 
thresholds  used  in  the  model  and  the  incorporation  of  pilot  behavioral  factors  influencing  the 
selection  of  vehicle  deceleration  schedules  for  various  runway  lengths.  3)  added  flexibility  in 
the  turnoff  angle  parameters  (i.e.,  turnoff  angle  and  other  lateral  spacing  restrictions),  and 
We  now  try  to  address  each  one  of  these  topics  in  more  detail  pointing  out  some  of  the 
obstacles  and  methods  that  could  be  used  to  implement  these  recommendations  at  a  later 
research  stage. 


8.1  Postoptimization  Algorithm 


A  post-optmization  processor  that  could  factor  day-to-day  aircraft  traffic  mix  variability 
and  environmental  conditions  could  be  a  very  practical  addition  to  the  existing  model.  The 
major  constraint  to  this  seems  to  be  the  computer  storage  limitations  required  to  handle  the 
larger  size  matrices  generated  by  this  new  postoptimization  process.  Currently  REDIM  is 
limited  to  50  aircraft/runway  condition  pairs  in  a  single  run  to  limit  the  numerical  computations 
to  a  manageable  level.  If  a  giobal  optimization  scheme  is  be  implemented  under  varying  input 
parameters  throughout  the  life  cycle  of  the  facility  this  would  necessitate  complete  knowledge 
of  the  time  variations  of  some  parameters  accounted  for  in  REDIM  throughout  the  period  of 
interest  (i.e.,  the  airport  design  life-cycle).  This  of  course  could  only  be  done  for  a  few 
variables  such  as  aircraft  mix  and  airfield  environmental  conditions  (i.e.,  airfield  temperature, 
wind  conditions,  etc.)  in  order  to  maintain  a  reasonable  matrix  size  to  execute  the  problem 
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on  a  personal  computer.  Another  valuable  alternative  in  this  regard  could  also  be  the  addition 
of  an  iterative  procedure  that  would  search  for  a  user  defined  ROT  threshold  value  to  be  used 
as  runway  exit  location/geometry  design  parameter.  This  procedure  could  be  executed  in 
severa'  steps  allowing  at  least  two  parameters  to  be  varied  independently  to  achieve  the 
desired  ROT  value.  Parameters  of  great  influence  in  this  regard  are  the  number  of  exits  and 
the  exit  speeds  associated  with  each  one  of  them.  A  typical  searching  algorithm  to  achieve 
a  "goal"  ROT  value  would  cycle  these  two  parameters  sequentially  until  the  desired  ROT  value 
is  obtained. 


8.2  human  Behavioral  Factors 


Another  aspect  deserving  attention  in  this  section  is  that  dealing  with  some  of  the  safety 
margins  and  assumptions  made  in  the  present  modeling  effort.  In  the  overall 
conceptualization  of  REDIM  safety  margins  were  implemented  in  some  of  the  dynamic  module 
subroutines  to  account  for  the  usual  uncertainties  associated  with  manual  control  tasks,  such 
as  the  landing  of  an  aircraft,  the  activation  of  braking  devices,  etc.  However,  the  reduction  of 
these  uncertainties  could  significantly  reduce  the  runway  occupancy  time  (ROT)  by  reducing 
the  margins  of  safety  needed  to  cope  with  the  original  assumptions.  This  phenomena  is 
similar  to  the  anticipated  reductions  in  the  aircraft  interarrival  time  (IAT)  to  the  runway 
threshold  through  an  improvement  of  the  aircraft  delivery  accuracy  (e.g.,  by  reducing  the  final 
approach  IAT  separation  buffers).  The  underlying  assumptions  made  in  this  model  have  tried 
to  establish  a  good  balance  between  operational  safety  and  the  efficiency  of  the  runway 
subsystem.  This  compromise  was  necessary  because  the  model  is  expected  to  be  applied  in 
a  variety  of  scenarios  where  the  manual  control  uncertainties  could  be  quite  high.  That  is,  the 
model  could  be  either  applied  to  small  community  airports  where  the  proficiency  and  accuracy 
of  the  pilots  might  dictate  slightly  larger  safety  margins  or  to  large  transport-type  airports 
where  an  increased  number  of  automated  landing  rollout  operations  could  take  place  in  the 
future.  It  is  expected  that  REDIM  will  be  calibrated  with  the  help  of  simulation  and 
experimental  results  in  order  to  gain  more  confidence  in  the  output  results  of  the  model.  This 
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calibration  is,  in  fact,  one  of  the  most  important  steps  to  follow  the  development  of  REDIM.  It 
is  anticipated  that  the  second  phase  of  this  research  will  devote  time  to  validate  the  current 
model.  It  should  be  clearly  understood  by  the  analyst  that  scenario-specific  factors  such  as 
obstructions,  runway  length,  lighting  conditions,  etc.  could  affect  the  pilot's  behavior  to 
execute  manual  landings.  For  example,  it  is  well  known  that  the  runway  exit  location  and 
length  have  a  large  influence  in  ROT  as  pilots  adjust  their  piloting  behavior  under  scenario 
specific  circumstances  such  as  displaced  thresholds  and  short  runways.  Therefore  a  series 
of  empirical  observations  are  recommended  in  the  future  in  order  to  modify  REDIM  to  account 
for  some  of  these  human  operational  factors. 


8.3  Turnoff  Angle  Parameters 


The  turnoff  angle  plays  a  very  important  role  in  the  estimation  of  the  runway  occupancy 
times  (RCT's).  Through  simulations  it  can  be  shown  that  as  much  as  25%  of  the  runway 
occupancy  time  is  due  to  the  turnoff  for  high  speed  exits.  As  such  it  is  advisable  to  add  more 
fexibility  to  the  model  by  allowing  the  user  to  vary  the  turnoff  angle.  This  new  addition  will  be 
highly  beneficial  for  analysis  involving  airport  improvements  where  severe  lateral  separation 
restrictions  pose  a  problem.  In  those  cases  the  analyst  could  specify  small  turnoff  angles  to 
achieve  a  desired  taxiway-turnoff  intersection  speed.  In  practice  it  has  been  shown  (Fig  6  2; 
that  turnoff  exit  angles  lower  that  18  degrees  do  not  reduce  the  runway  occupancy  time  as  the 
aircraft  travels  for  large  periods  of  time  on  the  turnoff.  This  in  turn  reduces  the  ROT  times  to 
a  extend  provided  that  the  exit  angle  is  not  reduced  below  20  degrees  which  seems  to  be  a 
compromise  bethween  minimum  ROT  and  the  time  to  clear  the  runway.  It  is  suggested  that 
the  final  exit  angle  should  be  made  variable  in  order  to  provide  the  user  a  mechanism  to 
design  fast  turnoffs  under  drastic  lateral  restrictions  (i.e.,  the  presence  of  a  parallel  taxiway). 
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Appendix  B.  Optimality  Through  a  Discrete  Search 

Suppose  there  exist  feasible  ranges  r  (r=  1  to  R)  for  feach  aircraft-surface  condition 
combination,  where  ROT's  within  each  range  are  increasing  from  left  to  right.  Assuming  that 
N  exits  are  to  be  located  on  the  runway  at  any  points  such  trial  there  is  at  least  one  exit  for 
each  feasible  range,  and  that  the  exits  are  separated  by  at  least  a  distance  of  Dm,„.  the  optimal 
exit  locations,  which  minimize  the  weighted  sum  of  ROT,  can  be  found  from  a  finite  collection 
of  points.  Let  L,  and  R,  be  the  left  hand  and  the  right  hand  interval  end  points  for  the  range 
r,  respectively.  L,  and  R,  are  actually  distances  measured  from  the  start  of  the  active  runway 
threshold.  Define  a  set  of  breakpoints  as  points  on  the  runway  which  are  of  the  type 
Lr  +  q  D min  for  q  S  0  and  integer  valued,  for  r  =  1  to  R.  Then  the  optimal  locations  are  found 
from  the  set  of  breakpoints  by  the  following  theorem. 

THEOREM  1  Assume  that  N  is  large  enough  so  that  the  above  problem  has  a  feasible 
solution.  Then  at  optimality,  each  location  will  coincide  with  some  breakpoint. 

PROOF  We  will  prove  this  by  induction  on  the  exit  index.  Consider  the  leftmost  exit 

location.  This  exit  must  coincide  with  L,  for  some  r  e  {1 . R}  because  if  not 

by  sliding  its  location  leftwards  until  it  coincides  with  such  a  location,  v/e  wi!i 
maintain  feasibility  (since  all  aircraft  which  could  take  this  exit  can  continue 
to  do  so),  and  the  objective  value  will  strictly  improve.  Inductively,  suppose 
that  the  result  is  true  for  the  location  of  exit  1, 

. .  .  ,  t,  and  consider  exit  t  + 1,  where  t  {1, . . . ,  N-1).  If  exit  t  + 1  coincides 
with  some  L,  for  R  e  {1,  ...  ,  R},  then  the  result  is  true.  If  exit  t  +  1  is  at  a 
distance  Dmi„  from  exit  t  to  its  left,  then  by  the  induction  hypothesis,  and  the 
construction  of  breakpoints,  the  result  is  again  true.  If  neither  of  these  cases 
holds,  then  we  can  slide  the  location  of  exit  t  +  1  leftwards  until  one  of  these 
conditions  holds,  thereby  maintaining  feasibility  and  improving  the  objective 
value.  Hence,  the  result  must  be  true  for  the  location  of  exit  t  +  1,  and  this 
completes  the  proof. 
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COROLLARY  1  For  any  pair  of  exits  t  and  1 4  1  separated  by  a  distance  greater  than  Dm,„,  the 
location  of  exit  t  + 1  must  lie  in  {  Li .  L*} 

PROOF  Evident  from  the  proof  of  Theorem  1. 

COROLLARY  2  Given  that  the  ROT'S  are  nondecreasing,  rather  than  strictly  increasing  within 
each  feasible  range,  there  exists  an  optimal  solution  in  which  the  exit  locations 
coincide  with  the  defined  breakpoints. 

PROOF  Evident  from  the  proof  of  Theorem  1. 

COROLLARY  3  (Improvement  problem)  Given  existing  exit  locations  at  points  D, . D, ,  define 

additional  breakpoints  as  the  points  D>  4-  q  Dmi«forq  s  1  and  integer,  for  i  =  1, 
...  ,  e.  Furthermore,  delete  from  the  set  of  breakpoints  thus  defined,  those 
which  lie  at  a  distance  less  than  D from  an  existing  exit  location  (  on  either 
side  of  it).  Then  again,  any  optimal  solution  will  have  the  new  exit  locations 
coinciding  with  these  defined  breakpoints. 

PROOF  Can  be  constructed  similar  to  that  of  Theorem  1. 

REMARK  By  Corollary  1,  for  N  and  Dm )n  small  enough,  optimal  locations  of  exits  will 

coincide  with  the  points  L,,  r=1 . R.  For  larger  values  of  these 

parameters,  the  other  breakpoints  will  begin  to  play  a  role.  This  is  of 
consequence  since  the  points  L,  represent  the  critical  locations  given  by  the 
simulations  of  aircraft  landing  movement.  Also,  given  our  emphasis,  Corollary 
3  is  of  most  importance. 
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Appendix  C.  Approximation  of  Turnoff  Times 

The  estimation  of  the  turnoff  time  plays  a  very  important  role  in  the  mathematical 
optimization  module  of  REDIM  as  the  dynamic  programming  technique  used  tries  to  minimize 
a  time  related  performance  index.  It  was  said  in  Chapter  4  of  this  report  that  in  order  to  save 
valuable  computational  time  it  was  necessary  to  approximate  the  time  spent  by  aircraft  in  the 
turnoff  maneuver  under  two  scenario  conditions  (dry  and  wet).  Furthermore,  every  secondary 
candidate  solution  (i.e.,  those  generated  from  the  actual  aircraft  landing  simulations  to  comply 
with  Bellman's  principle  of  optimality  as  explained  in  Appendix  B)  has  an  associated  turnoff 
time  (TOT)  for  every  aircraft  and  scenario  condition  and  thus  the  estimation  of  these  times 
would  consume  large  amounts  of  time  if  performed  through  the  complete  simulation  scheme 
used  to  estimate  primary  candidates  and  described  in  Section  3.3  of  this  report. 

Since  the  geometry  for  every  primary  candidate  is  completely  known  from  the  simulation 
results  it  is  possible  to  extract  two  representative  values  of  the  radius  of  curvature,  R1  and 
R2.  to  approximate  the  turnoff  geometry  until  the  aircraft  has  cleared  the  runway  as  depicted 
in  Fig.  C.l.  It  should  be  emphasized  that  although  this  is  an  approximation  the  results  are 
usually  accurate  if  R1  and  R2  are  selected  appropriately.  In  the  late  fifties  Horonjeff  (Horonjeff, 
1959]  used  this  scheme  to  approximate  high-speed  turnoff  tracks  with  satisfactory  results. 

Through  hundredths  of  simulations  of  the  REDIM  model  its  was  observed  that  the  values 
of  R1  and  R2  could  be  extracted  from  the  turnoff  simulation  as  a  function  of  time  and  aircraft 
category.  This  segmentation  per  category  was  somewhat  expected  from  equations  3.13-3.18 
in  Section  3.4  if  one  realizes  that  the  aircraft  turning  capability  is  related  to  the  inertia, 
centripetal  and  scrubbing  forces  resisting  the  aircraft  turning  motion.  Results  depicting  the 
time  rate  of  change  variations  of  the  radius  of  curvature  for  representative  aircraft  using 
REDIM  are  shown  in  Fig.  C.2.  It  was  then  decided  through  examination  of  all  the  data  to 
estimate  R1  as  the  instantaneous  radius  of  the  curvature  occuring  one  second  after  the  turning 
maneuver  started  whereas  R2  was  varied  selectively  between  four  and  six  seconds  depending 
upon  the  aircraft  category.  The  four-second  R2  is  used  with  category  A  aircraft  which  display 
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Figure  C.1  Turnoff  Time  Approximation  Nomenclature. 


very  fast  behavior  in  the  turnoff  dynamics  whereas  the  six-second  R2  is  used  to  predict  heavy 


transport  aircraft  turnoff  dynamics  having  larger  time  lags  to  achieve  a  "steady-state"  radius 


of  curvature,  in  REDIM  nomenclature  these  times  are  labeled  as  easement  curvature  time, 


TR1,  and  steady-state  curvature  time.  TR2. 


Once  the  exact  turnoff  path  is  known  the  next  step  is  to  estimate  the  time  required  to  clear 


the  runway.  This  is  done  under  the  assumption  that  a  turning  aircraft  decelerates  due  to 


rolling  friction  alone.  Actual  aircraft  speed  measurements  performed  by  by  Horonjeff 


(Horonjeff  et  al,  1959,  1960]  and  Hosang  (Hosang,  1978]  in  high-speed  taxiways  show  nearly 


constant  deceleration  rates  similar  to  those  associated  with  a  moderate  value  of  rolling 


friction  alone.  This  can  be  attributed  to  the  small  aircraft  castor  angles  present  while 


negotiating  a  high-speed  turnoff.  A  conservative  value  of  Fm,  of  .03  has  been  used  throughout 


the  program  to  model  the  rolling  friction  deceleration  rates  experienced  by  every  aircraft. 
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Radius  of  Curvature  Time  Histories 


Figure  C.2  Time  Variations  of  Radius  of  Curvature  for  Representative  Aircraft. 


Fig.  C.1  illustrates  how  the  turnoff  time  is  estimated  using  two  simple  radii  of  curvatire  to 
approximate  the  actual  turnoff  track.  Two  turnoff  arcs  with  lengths  Li  and  L2  are  defined  as 
follows, 


Li  =  *1  0y 

{Cl} 

L 2  —  /?2  $2 

{C.2} 

where,  Li  and  lt  represent  the  turnoff  characteristic  lengths.  R ,  and  /?*  are  the  radii  of 
curvature  defining  the  turnoff,  and  0,  and  02  are  the  arcs  defined  by  /?i  and  Rz,  respectively 
measured  in  radians.  The  turnoff  arcs  are  characteristic  for  each  aircraft  since  the  transition 
and  runway  clearance  point  are  aircraft  speed  and  geometry  dependent.  A  further 
simplification  regarding  the  easement  length,  Li  can  be  introduced  using  results  derived  from 
highway  geometric  design  principles  where  the  length  of  a  spiral  transition  curve  Li  is  made 
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a  function  of  exit  speed.  Horonjeff  later  on  showed  that  a  short  transition  spiral  could  well  be 
approximated  with  a  large  radius  of  curvature  segment  and  this  approximation  is  easily 
implemented  in  the  model  [Horonjeff,  1959). 


L1  =  2.914  +  3.1701  V, 


exit 


for  8  m/sec.  <,  "exit  ^  45  ml  sec. 


[C  3) 


The  aircraft  speed  at  the  transition  point  between  the  two  radii  of  curvature  is  obtained  from 
Eqn.  C.4  whereas  the  speed  at  the  rurlway  clearance  point  is  shown  in  Eqn.  C.5 


/ tran~  exit  2  fron 


{C.4} 


V final  ~  (V Iran  2  fr0//  Lj} 


{C.5} 


where.  V,,*  is  the  desired  aircraft  exit  speed  (m./sec.),  V,,*,  is  the  transition  speed  (m./sec), 
Vfmii  is  the  final  speed  at  the  runway  clearance  point,  (Xe  ,YC  ),  g  is  the  gravity  constant 
(m./sec.-sec.)  and  f„„  is  the  rolling  friction  coefficient  (dimensionless).  The  travel  time  across 
each  of  the  turnoff  segments  is  estimated  as  shown  in  Eqns.  C.8  and  C.7. 


Ti  = 


' exit 


+  V, 


tran 


{C.6} 


'tran 


+  v< 


final 


{C- 7} 


where,  Ti  and  Tt  are  the  travel  times  from  the  start  of  the  turnoff  to  the  transition  point 
(X, ,  V,)  and  from  transition  point  to  runway  clearing  point  (  Xe ,  Yc  ),  respectively.  The  total 
turnoff  time  is  the  summation  of  these  two  previous  contributions. 

T,0i  +  T2  {C.8} 


where,  7,„  is  the  turnoff  time  until  clearing  the  runway.  This  procedure  to  estimate  the 
turnoff  time  is  implemented  for  the  secondary  candidates  whose  locations  are  q(Dmn)  meters 
away  from  primary  candidate  solutions  (for  q  =  1,2, ....  N)  as  explained  in  Chapter  4  of  this 
report.  It  should  be  noticed  that  the  secondary  candidate  solutions  obtained  for  small  aircraft 
far  downrange  from  an  active  threshold  will  usually  be  unfeasible  for  large  aircraft  since  these 
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will  nol  be  able  to  negotiate  the  turnoff  with  the  desired  margin  ot  safety 
even  more  the  candidate  set  to  be  used  in  the  optimization  module. 


This  process  reduces 
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